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Abstract  of  Dissertation  presented  to  the  Graduate  School 

of  the  University  of  Florida  in  Partial  Fulfillment  of  the 

Requirements  for  the  Degree  of  Doctor  of  Philosophy 

TURBIDITY  CURRENTS  AND  SEDIMENTATION  IN 
CLOSED-END  CHANNELS 

By 

CHUNG-PO  LIN 

May,  1987 

Chairman:   Dr.  Ashish  J.  Mehta 
Major  Deparment:   Civil  Engineering 

As  a  result  of  concern  over  sedimentation  and  water  quality  in 
closed^end  channels  such  as  docks,  pier  slips  and  residential  canals, 
an  investigation  on  the  front  behavior  of  turbidity  currents  and 
associated  sedimentation  in  a  closed»end  channel  was  carried  out.  The 
foeus  of  the  present  study  was  on  the  case  of  a  closed-end  channel  of 
rectangular  cross  section  and  horizontal  bed  connected  orthogonally  to 
a  main  channel,  where  turbulent  flow  and  high  suspension  concentration 
existed.  Experiments  were  conducted  in  two  specially  designed  flume 
systems  to  investigate  front  characteristics  of  turbidity  currents, 
flow  regime,  and  sediment  deposition. 

A  group  of  dimensionless  parameters  for  each  physical  aspect  of 

concern  derived  from  dimensional  analysis  provided  the  basis  for 

selecting  the  neaaurement  items  in  laboratory  experiments  and  for 

presenting  the  results.  A  two-dimensional,  explicit,  coupled  finite 

difference  numerical  model  for  simulating  vertical  flow  circulation  and 

sediment  ion  in  a  closedfend  channel  was  made.  In  addition,  analytic 

xix 


cal  developments  for  the  longitudinal  distribution  of  mean  concentra*- 
tion  below  interface,  flow  field,  the  maximum  rise  of  water  surface, 
sediment  flux  through  the  entrance)  rind  tide-induced  deposition  were 
also  attempted. 

No  significant  differences  were  found  between  the  front  shapes  of 
turbidity  currents  and  non*settiing  gravity  currents.  The  rate  of  the 
decrease  of  front  speed  of  turbidity  currents  with  distance  was  faster 
than  that  of  nonHsettling  gravity  currents.   Characteristics  which 
showed  an  exponential-type  decrease  with  distance  in  the  closedfend 
channel  include  front  speed,  mean  concentration  below  interface,  front 
concentration,  settling  velocity  of  suspension,  dispersed  particle  size 
of  deposit,  and  deposition  rate.  The  relative  settling  velocity, 
i.e.  the  ratio  of  particle  settling  velocity  to  the  densimetric 
velocity,  was  found  to  be  the  best  parameter  to  interpret  the  different 
ces  between  tests.  Durations  of  front  propagation  in  the  initial 
adjustment  phase,  inertial  selffisimilar  phase,  and  viscous  self.-simiiar 
phase  were  found  to  be  mainly  dependent  on  sediment  type  and  the 
relative  settling  velocity.  A  relationship  for  predicting  the  sediment 
flux  into  a  closed-end  channel  through  entrance  was  found.   According 
to  the  relationship,  the  sediment  flux  is  proportional  to  H1/2  and 
Cf  ,  where  H  is  the  total  water  depth  and  C,  is  the  depth-mean 
concentration  at  the  entrance  of  the  closed'end  channel.  It  was  shewn 
that  the  numerical  model  could  simulate  the  suspension  concentration 
distribution  in  the  channel  satisfactorily. 
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CHAPTER  1 
INTRODUCTION 

1 . 1  Turbidity  Currents  and  Sedimentation  in  a  Closed-End  Channel 
Sedimentation  problems  in  pier  slips,  tidal  docks,  and  closed-end 
channels  have  been  mentioned  by  many  investigators  (e.g.,  Simmons, 
1966;  O'Dell,  1969;  Kanless,  1975;  Hoffman,  1982).   In  order  to 
maintain  the  required  depth  of  navigation  in  these  shoaled  channels, 
frequent  costly  dredging  becomes  inevitable.   Recently,  increases  in 
dredging  costs  and  the  problem  of  dredged  spoil  disposal  have  motivated 
the  search  for  other  methods  of  sediment  control  (Hoffman,  1982; 
Ballard  et  al.,  1985),   Because  of  the  complexity  and  variety  of  the 
factors  involved  in  the  shoaling  processes  at  each  shoaled  area,  thus 
far  no  efficient  and  standard  method  of  sedimentation  control  is  avail- 
able.  In  fact,  the  prevention  of  sedimentation  at  a  particular  site 
must  be  based  on  a  knowledge  of  the  local  dominant  forces  and  associat- 
ed mechanisms. 

The  main  causative  factors  of  the  sedimentation  problem,  for 
example,  in  ?oastal  closed-end  residential  canals  of  southwest  Florida, 
suggested  by  van  de  Kreeke  et  al.  (1977),  are  density  and  wind-induced 
flows,   in  tidal  docks  on  the  Mersey  Estuary,  England,  investigated  by 
Haliiwell  and  O'Dell  (1970),  tida]  flows  and  density  currents  are  the 
primary  causes  of  sedimentation.   Accordingly,  tide,  wind,  and  density 
currents  are  the  most  likely  physical  factors  involved  in  shoaling 
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processes  in  areas  such  as  pier  slips,  tidal  docks,  and  clcsed^end 
channel.;.   Density  currents  include  currents  Induced  by  salinity, 
temperature,  and  turbidity  gradients.  There  is  a  large  body  of  papers 
describing  mechanisms  related  to  tide-  and  wind-driven  circulation  as 
well  as  salinity-  and  thermally-induced  density  currents  in  bays  and 
estuaries.   However,  less  consideration  seems  to  be  given  the  mechanism 
of  turbidity-induced  density  currents  (turbidity  currents)  in  shallow 
water  environments.   In  some  areas  with  high  suspended  sediment 
concentrations,  turbidity  currents  may  play  an  important  role  in  the 
sedimentation  problems  of  channels. 

Inside  a  closed-end  channel,  suspended  sediments  tend  to  settle 
onto  the  channel  bottom  because  of  the  quiescent  conditions.  Subsequ- 
ently, the  suspension  concentration  inside  the  channel  will  become 
lower  than  that  of  the  outside  water  body.   Consequently,  a  turbidity 
gradient  between  outside  waters  and  inside  waters  appears,  resulting  in 
the  intrusion  of  a  turbidity  current  into  the  channel  along  the  bottom. 
During  the  penetration  of  the  turbidity  current,  the  clear  water 
overlying  the  layer  of  turbid  water  is  forced  to  flow  out  of  the 
channel  due  to  the  water  surface  gradient  generated  by  the  intrusion  of 
denser  fluid.   In  essence,  therefore,  a  vertical  circulation  of  fluid 
is  set  up.   Such  a  circulation  pattern  can  provide  exchange  between 
outside  and  inside  waters,  and  may  affect  the  water  quality  in  the 
channel . 

One  significant  feature  of  the  turbidity  current  is  that  the 
density  current  caused  by  sediment  in  suspension  itself  tends  to  keep 
sediment  in  suspension  (Bagnold,  1962).   The  type  of  turbidity  current 
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•--depositing,  eroding,  or  lutosuspmalon  (neither  deposting  nor 
eroding)  —  depends  on  thi  p  .  iertl%a  of  suspended  sediment,  the  front 
speed  of  the  turbidity  current,  and  the  period  of  consolidation  of  the 
bed  material  (Pallesen.  1983).  In  the  present  study,  the  turbidity 
currents  considered  are  depositing  currents,  which  carry  suspended 
sediments  from  the  outside  area  into  the  closed-end  channel  and  deposit 
them  on  the  channel  bed. 

In  this  study,  turbidity  currents  and  associated  sedimentation  in 
closed-end  laboratory  flumes  were  investigated.   Prior  to  the  start  of 
a  test,  steady  flow  and  sediment  concentration  were  established  in  the 
main  channel  with  the  entrance  gate  closed,  separating  the  quiescent 
clear  water  in  the  closed-end  channel  from  turbid  water  in  the  main 
channel.  The  gate  was  lifted;  then  measurements  were  made  as  the 
intrusion  turbidity  front  propagated.  Characteristics  of  resulting 
deposits  were  determined  subsequently.  A  schematic  diagram  showing  the 
geometry  of  the  laboratory  flume  and  an  intruding  turbidity  front  in 
the  closed-end  channel  is  shown  in  Fig.  1.1(a)  and  (b).   The  gyre  zone 
shown  in  the  figure  indicates  a  circulation  cell  developed  just  inside 
the  entrance.  This  circulation  was  induced  by  the  laterai  shear  force 
exerted  by  the  main  flow  at  the  flow  boundary  near  the  entrance.   This 
gyre  zone  was  relatively  short  in  longitudinal  extent.   The  suspended 
sediment  was  fairly  well  mixed  vertically  in  the  gyre.   Beyond  the 
gyre,  the  turbidity  front  was  stratified.   Definitions  of  the  other 
te'-ms  indicated  in  the  figures  dre  given  in  section  ^  1,   it  suffices 
to  note  here  that,  the  major  influencing  parameters  are  the  sediment 
properties,  channel  geometry,  and  suspension  concentration  at  the 
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entrance  of  the  olosed^end  channel,  which  is  connected  to  a  nair:  f  o« 
channel . 

1.2  Approach  to  the  Problem 

In  order  to  study  the  front  behavior  of  the  turbidity  current  and 
the  associated  sedimentation  in  a  closed»-end  channel,  dimensional 
analysis  was  first  carried  out  to  relate  the  phenomena  of  Interest  to 
their  influencing  dimensionless  parameters.  Next,  physical  and 
mathematical  models  were  developeded  to  Investigate  the  general 
characteristics  of  the  problems  of  concern.  An  attempt  was  made  to 
solve  the  equations  of  eonvective«dif fusion  and  momentum  for  concen*- 
tration  distribution  and  flow  regime  in  the  channel  under  a  steady 
state  condition. 

Dimensional  analysis  of  several  important  aspects  of  this  investi- 
gation was  carried  out  using  the  Buckingham  <-  it  theorem  which  provided 
a  group  of  dimensionless  parameters.  These  dimensionless  parameters 
can  aid  understanding  of  the  phenomena  of  concern  in  a  physical  sense, 
to  reduce  the  number  of  measurement  items  in  laboratory  experiments, 
and  to  express  the  investigation  results  in  a  compact  format. 

The  study  of  the  problem  by  physical  model  includes  two  sets  of 
laboratory  experiments,  which  were  conducted  in  two  differentt-sized 
flume  systems.   Each  measurement  taken  in  the  flume  tests  and  the 
presentation  of  the  results  were  based  on  the  dimensionless  parameters 
obtained  from  the  dimensional  analyses  mentioned  above.  Two  types  of 
f  ine**grained  sediments,  cohesive  and  cohesionless,  were  used.  The 
settling  velocities  and  median  particle  sizes  of  these  sediments  were 
determined  by  settling  column  analysis  (McLaughlin,  1959)  and  hydro- 
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meter  or  pipette  analysis  fASTM,  1983;  uuy,  1975).  Owing  to  the 
complexity  Involved  in  3caltng  the  properties  of  fine»-grained  sediment 
and  internal  density  oonditions,  no  presently«known  raouel  technique  can 
reproduce  tne  flelc:  phenomena  in  the  small-scale  models  (Tesaker, 
1975).  Therefore,  the  results  of  this  flume  study  are  regarded  as 
indicative  of  a  qualitative  approach  to  the  prototype  scale  problem. 

A  two«dimensional  (x  and  2  directions),  explicit,  finite  differ- 
ence mathematical  numerical  model  was  developed.  In  the  model,  the 
coupled  momentum  and  convective«diffusion  equations  govern  the  movement 
of  water  and  suspended  sediments.  A  common  procedure  used  for  simulate 
ing  tide*-  and  wind^lnduced  sediment  transport  in  estuaries  by  a 
numerical,  model  is  to  solve  the  flow  field  directly  from  the  momentum 
equation,  then  the  transport  of  suspended  sediment  can  be  obtained  by 
solving  the  convective'-diffusion  equation  with  the  known  flow  field. 
In  the  present  study,  since  the  motivating  force  in  the  momentum 
equation  is  the  density  gradient  resulting  from  suspended  sediment 
concentration,  the  flow  regime  and  the  distribution  of  suspended 
sediment  closely  affect  each  other.  Hence  the  two  governing  equations 
must  be  coupled  and  solved  simultaneously.  The  four  coefficients  M 
longitudinal  and  vertical  momentum  diffusion  coefficients,  and  longitu- 
dinal and  vertical  mass  diffusion  coefficients  N<  were  calibrated  until 
the  model's  predictions  showed  a  reasonable  agreement  with  the  measured 
data  from  a  selected  laboratory  test.  With  the  chosen  coefficients, 
thr  model's  predictions  of  sediment  distribution  and  deposition  rate 
were  compared  with  the  corresponding  experimental  results  for  other 
te3ts. 
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Five  physical  aspects  of  concern  in  the  ciosed"end  channel  at  r.oe 
steady  state  condition  are  discussed  in  the  analytical  development 
section.  They  include  1 )  the  longitudinal  variation  of  mean  suspension 
concentration  below  the  interface;  2)  the  flow  regime  in  the  channel; 
3)  sediment  flux  into  the  channel;  4)  the  maximum  rise  of  water  surface 
in  the  channel;  and  5)  the  tide»-induced  to  turbidity  currentMnduced 
deposition  ratio.  The  longitudinal  variation  of  mean  concentration 
below  the  interface  is  obtained  by  solving  the  verticallyrtintegrated 
form  of  the  convectiver-dif  fusion  equation  at  steady  state.   Based  on 
this  result,  the  steady  state  equation  of  motion  is  solved  for  the  flow 
field  inside  the  closed<-end  channel.  The  formula  for  predicting  the 
sediment  flux  into  the  channel  through  the  entrance  is  obtained  as  a 
product  of  the  mean  velocity  and  suspension  concentration  below  the 
interface  at  the  channel  entrance;  furthermore,  the  mean  velocity  is 
related  to  the  densimetric  velocity,  which  is  the  characteristic 
velocity  induced  by  density  differences  (Keulegan,  1957). 

A  stationary  sediment  wedge  can  be  established  by  a  force  balance 
between  the  force  induced  by  the  dense  fluid  at  the  entrance  of  the 
closed-end  channel  and  the  hydrostatic  pressure  of  clear  water  with  an 
elevated  water  surface  in  the  region  beyond  the  toe  of  the  stationary 
sediment  wedge.  A  formula  for  predicting  the  magnitude  of  the  rise  of 
water  surface  needed  for  the  establishment  of  a  stationary  sediment 
wedge  in  the  closed-end  channel  is  derived  based  on  the  conservation  of 
momentum. 

The  ratio  of  tide-induced  deposition  to  turbidity  current-induced 
deposition  in  the  closed«end  channel  can  be  computed  based  on  the 
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concept  of  tidal  prism  and  the  sediment  flux  of  turbidity  current 
mentioned  above. 

1.3  Scope  of  Investigation 

Three  major  foci  of  this  study  were  1)  front  behavior  of  the 
turbidity  current,  2)  mechanics  of  a  stationary  sediment  wedge,  3) 
turbidity-induced  sedimentation. 

Front  behavior  of  turbidity  currents  (which  play  a  significant 
role  in  sediment  transport)  was  examined  primarily  by  laboratory  flume 
tests.  The  investigation  included  initial  front  speed,  characteristics 
of  front  propagation,  and  front  shape.  Results  are  compared  with  the 
corresponding  front  behavior  of  salinity-induced  density  currents 
reported  in  previous  studies. 

For  the  stationary  sediment  wedge,  the  longitudinal  variation 
of  mean  concentration  (mean  sediment  concentration  in  the  layer  below 
the  interface)  along  the  closed-end  channel  was  Investigated  experi- 
mentally and  analytically.  The  flow  pattern  in  the  channel  was 
analytically  derived  from  the  equation  of  motion.  In  addition, 
formulas  for  predicting  sediment  flux  into  the  closed^end  channel  and 
for  computing  the  magnitude  of  the  water  surface  rise  in  the  region 
beyond  the  toe  of  sediment  wedge  were  derived,  and  the  predictions  were 
compared  with  experimental  data. 

Characteristics  of  turbidity-induced  sedimentation  in  a  closed-end 
channel  were  investigated  experimentally  and  by  means  of  the  numerical 
model.  Herein,  the  longitudinal  variations  of  the  deposition  rate, 
median  dispersed  particle  size  and  settling  velocity  of  deposits  along 
the  flume  are  presented  and  discussed  in  Chapter  5. 
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t  RUfiber  of  publications  are  reviewed  In  Chapter  2.     These  papers 
describe  the  front  behavior  of  density  currents  'most  are  salinlty 
induced  density  currents),  the  steady-state  sediment  wedge  in  a 
closed-end  channel,  mathematical  models  for  estuarial  sediment  trans- 
port, and  the  depositional  behavior  of  cohesive  sediments. 

Chapter  3  presents  and  discusses  a  dimensional  analysis  for 
several  subjects  of  concern  in  the  present  study.  Finite  difference 
formulations  of  the  coupled  governing  equation  in  the  mathematical 
model  are  developed  and  the  stability  conditions  of  the  numerical 
scheme  are  examined.  Analytical  solutions  for  mean  sediment  concentra- 
tion below  interface  and  water  velocity  at  steady  state  are  obtained 
from  the  equations  of  convective"diffusion  and  momentum.  In  addition, 
the  solutions  for  sediment  flux  into  the  closed-end  channel,  the  rice 
of  water  surface  in  the  closed'-end  channel  at  steady  state,  and 
tide-induced  to  turbidity  current-induced  deposition  ratio  are  present- 
ed. 

In  Chapter  H,    the  set-up,  equipment,  test  materials,  and  procedure 
of  two  sets  of  laboratory  experiments  are  described.   The  first  set  was 
performed  in  a  9.1  m  long  flume  at  the  U.  S.  Army  Corps  of  Engineers 
Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi;  the  second 
set  was  conducted  in  a  14.7  m  long  flume  at  the  Coastal  Engineering 
Laboratory  (COEL)  of  the  University  of  Florida. 

Chapter  5  shows  typical  experimental  observations  of  a  selected 
laboratory  test.  It  also  includes  the  comparison  between  experimental 
results  and  the  corresponding  predictions  obtained  by  the  calibrated 
numerical  model  and/or  analytical  approach.   Subsequently,  characteris- 
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tic  remits,  which  induce  the  data  of  ail  experimental  tests,  are 
preser'tea  mostly  in  a  dimenaionless  format.   Some  of  the  characteristic 
results  are  compared  with  the  corresponding  results  reported  by 
previous  studies  and/or  from  the  present,  analytical  approach.   In 
addition,  the  results  of  two  tidal  effect  tests,  (one  for  "flood  tide," 
one  for  "ebb  tide"),  and  a  reference  test  are  discussed. 

In  Chapter  6,  a  summary  of  this  study,  conclusions,  and  recommen- 
dations for  future  work  are  presented. 

In  Appendix  A,  the  spatial  grid  system,  flow  chart,  and  stability 
condition  for  the  mathematical  modeling  are  presented.   Appendix  B 
includes  the  derivations  for  the  relationship  between  suspension 
concentration  and  density  of  water-sediment  mixture,  and  the  flow 
regime  in  the  clcsed>-end  channel  at  steady  state.  In  Appendix  C,  a 
computation  of  the  maximum  water  surface  rise  in  the  closed-end  channel 
at  steady  state  for  a  select  laboratory  test  is  given.  The  characteri- 
stics of  closed-end  channel  deposits  are  presented  in  Appendix  D. 


CHAPTER  2 
LITERATURE  REVIEW 


2.1  Front  Behavior  of  Turbidity  Current 
2.1.1  Classification 

Gravity  currents,  sometimes  called  density  currents,  are  formed  by 
fluid  flowing  under  the  influence  of  gravity  into  another  fluid  of 
different  density.  Due  to  the  greater  weight  of  the  denser  fluid,  a 
larger  hydrostatic  pressure  exists  inside  the  current  than  in  the  fluid 
ahead,  which  provides  the  motive  force  to  drive  the  gravity  current. 
Gravity  currents  occur  in  many  natural  and  man-made  situations,  for 
example,  a  turbidity  current  resulting  from  a  landslide  at  the  sea 
bottom;  oil  spillage  on  the  sea  surface;  spreading  of  warm  water 
discharged  from  power  plants  into  rivers;  a  cold  air  front  confronted 
by  a  warm  air  front  in  the  atmosphere;  and  accidental  release  of  dense 
industrial  gases.  Because  of  the  wide  and  frequent  occurrence  of 
gravity  currents,  investigations  of  the  front  behavior  of  such  current 
have  continued  to  the  present  time  (O'Brien  &  Cherno,  1 931* ;  Von  Karman, 
19^0;  Ippen  &  Harleman,  1952;  Barr,  1963a,  1963b,  1967;  Middleton, 
1966a,  1966b,  1967;  Benjamin,  1968;  Simpson,  1982;  Akiyama  t   Stefan, 
1985).  In  these  studies,  the  front  shape,  which  indicates  the  coverage 
range  of  the  current,  and  the  front  speed,  which  shows  the  spreading 
rate  of  the  current,  were  of  primary  concern. 
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The  gravity  current  is  characterized  by  a  "head  wave"  where  the 
head  height,  h2,  Is  highei  thai,  neck  height,  h3  (see  Fig. 1.1(b)).  On 
the  rearward  side  of  the  nead  wave  there  is  a  highly  turbulent  zone 
suggestive  of  wave-breaKing  and  mixing.  These  features  appesr  to  be 
commor  Lo  rcany  physical  phenomena  that  may  be  clas3fied  as  gravity 
currents. 
2.1 .2  Front  Shape 

Von  Karman  (1910)  investigated  the  front  shape  of  the  gravity 
current  based  on  Stokes'  deduction  of  the  extreme  sharp*crested  shape 
of  water  waves,  and  concluded  that  the  stagnation  point,  denoted  as  S' 
in  Fig. 2.1,  lies  on  the  no^slip  bottom  boundary,  and  a  tangent  line 
drawn  from  the  point  S'  on  the  interface  of  the  front  has  a  slope  of 
it/ 3.  The  same  result  was  derived  in  Benjamin's  (1968)  extensive 
treatment  which  was  based  on  inviscid  fluid  theory.  However,  from 
laboratory  experiments,  the  stagnation  point  was  observed  to  be 
elevated  (Fig. 2. 2),  unlike  the  analytical  predictions  by  Von  Karman 
(1910)  and  Benjamin  (1968). 

From  their  laboratory  experiments,  Ippen  and  Harleman  (1952), 
Keulegan  (1958),  and  Middleton  (1966)  have  suggested  a  universal 
dimensionless  front  shape  profile  which  will  approximately  fit  all 
gravity  current  heads  of  different  sizes  and  velocities.  In  these 
profiles,  estimates  of  the  ratio  of  nose  height  to  head  height  h-|/h2, 
where  h-|  Is  the  nose  height  or  the  distance  from  the  bottom  to  the 
point  S',  and  h2  is  the  head  height,  shown  in  Fig. 2. 2,  vary  from  0.1  to 
0.27.   This  result  reveals  that  in  real  flows,  because  of  friction 
resistance  at  the  stationary  boundary,  the  lowest  streamline  in  the 
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Figure  2.1.  Front  Shape  of  Gravity  Current  From  Theoretical  Model 
(after  Von  Karraan,  19H0). 


•'lgure  2.2  Schematic  of  Mean  Flow  Relative  to  a  Gravity  Current 
Head  Subject  to  No«Slip  Conditions  at  Its  Lower 
Boundary  (after  Simpson,  1972). 
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flow  relative  to  the  head  must  be  towards  the  rear,  and  the  stagnation 
point  ?'  then  is  raised  some  small  distance  above  the  fioor  (Mg..'<.2). 

Simpson  (1972)  further  found  that  the  ratio,  h,/h?,  where  h,  is 
the  nose  height,  is  not  a  universal  constant,  but  is  dependent  on  a 
local  head  Reyno.as  number,  Reh  =--  ufh2/v,  where  Uf  is  the  front  speed 
of  the  gravity  current  and  u  is  the  kinematic  viscosity  of  water.  Bas- 
ed on  his  experiments,  he  proposed  a  relationship  between  tn/h7   and 
Reh,  according  to  which  h1/h2  decreases  as  Reh  increases.  In  Chapter 
5,  the  trend  between  ni/h2  and  Reh,  observed  in  the  laboratory  measure- 
ments of  the  present  study,  is  compared  with  the  results  of  previous 
investigations. 

Another  important  aspect  of  front  shape  is  the  ratio  of  h2/h3, 
where,  as  noted.  h2  ana  h3  are  the  heights  of  the  head  and  the  neck, 
respectively.  An  approximately  constant  value  of  2.08  was  obtained  by 
Keulegan  ftfjt'j  in  the  range  of  Ren  between  1  o3  and  10?,  where  Ren  - 
ufh3/v  designates  the  local  neck  Reynolds  number,  using  the  characte- 
ristic length.  h3.  Simpson  and  Britter  (1979)  found  empirically  that 
the  ratio  h2/h3  is  not  a  universal  constant,  but  depends  upon  the 
fraction  depth  fej/B,  where  H  is  the  total  water  depth.   A  comparison 
between  the  variation  of  h2/h3  with  Ren,  obtained  from  the  present  and 
previous  investigations,  is  given  in  Chapter  5. 

The  existence  of  a  "head  wave"  in  the  gravity  current  id  due  to 
the  excess  mass  flux  of  the  lower  fluid  brougnt  forward  to  the  head 
from  m   neck  -egion  and  lifted  upward  within  the  head,  tnen  returnea 
hack  toward  the  neck  region  by  means  of  interface  waves  breaking  just 
behind  the  head.   In  order  to  explore  the  intensity  of  mixing  oetween 
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the  two  fluids  at  the  hppr  I nterface,  which  is  related  tc  the  front 
head  height,  the  ratio  iif.u,,  where  u3  is  the  depth-averaged  velocity 
in  the  neck,  was  examined  by  Komar  (1977).  Baaed  on  his  experimental 
results,  Komar  pointed  out  that  uf/u3  is  dependent  on  the  densimetric 
Froude  Number,  FrA  ■  uf/[ApgH/pw]1/2>  where  Pw  ls  the  density  of  fresh 
water  and  g  is  the  acceleration  of  gravity.  Note  that  u3  is  greater 
than  Uf.  As  the  ratio  Uf/uj  becomes  smaller  with  increasing  FrA,  the 
ratio  n2/h3  will  become  larger  (i.e.  more  intensive  mixing  between 
two  fluids  will  take  place).  Since  the  rate  of  mixing  at  the  head 
between  the  two  fluids  results  from  the  excess  mass  flux  into  the  head; 
therefore  this  rate  can  be  computed  as  (ua  -  uf)ho. 
2.1 .3  Front  Speed 

The  initial  front  speed,  uf ( ,  in  lock  exchange  flow  of  saline 
water  has  been  investigated  by,  among  others,  O'Brien  and  Cherno 
(1931),  Yih  (1965)  and  Keulegan  (1957).   These  investigations  found 
that  in  a  mutual  intrusion  between  two  liquids  of  a  lock  exchange  flow, 
a  stagnation  point  appears  at  the  elevation  of  one-half  of  the  total 
water  depth  at  the  cross  section  of  the  barrier  (used  to  separate  two 
fluids  of  different  densities  before  the  test).   The  initial  front 
speed,  uf 1 ,  can  be  expressed  as 


Ap     1/2 

(  2.1  ) 


where  Cp  is  a  "universal"  constant,  and  Ap  and  pm  are  the  density 
difference  and  mean  density  between  the  two  liquids,  respectively.  In 
eq.2.1  the  square-root  ten;  was  denoted  as  the  densimetric  velocity, 


16 

uA,  by  Keulegan  (1957): 

Ap     1/2 

"A  -  £  —  8"  ]  (  2.2  ) 

Pm 

The  coefficient,  Cp,  in  eq.2.1  was  found  to  be  0.50  by  Yin  (1965)  using 

an  approach  based  on  the  conservation  of  energy,  and  0.16  from  the 

laboratory  experiments  conducted  by  Yin  (1965)  and  Keulegan  (1957). 

Initial  front  speeds  were  examined  in  laboratory  experiments  of  the 

present  study,  and  the  results  are  compared  in  Chapter  5  with  those 

obtained  by  Keulegan  (1957). 

Benjamin  (1968)  pointed  out  that  the  initial  front  speed  which 

occurs  after  removing  the  barrier  in  lock  exchange  flow  will  remain 

unchanged  over  a  distance,  during  which  the  gravity  current  propagates 

without  energy  losses.  After  the  front  propagates  beyond  a  certain 

distance,  energy  conservation  no  longer  holds,  because  of  interfacial 

as  well  as  bottom  friction.  During  energy  dissipation,  a  "head  wave" 

at  the  front  is  formed,  and  both  front  speed  and  head  height  decrease. 

A  dimensionless  factor,  FrA  -  ur/[Ap' gh3/Pw] 1 '2,  termed  the  loca] 

densimetric  Froude  number,  is  used  to  express  the  relationship  between 

the  front  speed,  Uf,  and  local  densimetric  velocity,  where  Ap'  is  the 

local  density  difference  between  the  liquid  in  the  front  head  and  tne 

liquid  ahead.   Keulegan  (1958)  found  that  Fr&  is  independent  of  the 

local  neck  Reynolds  number  ufh3/v  when  the  number  is  greater  than  100, 

but  is  a  function  of  urh3/v  when  Ufhy'v  <  100.   A  similar  conclusion 

that  the  dimensionless  front  speed  is  strongly  dependent  or,  the 

densimetric  Reynolds  number,  ReA  =  uAH/v,  as  ReA  <  order(103),  and 

Reynolds  number-independent  as  ReA  >  order(1o3)  was  made  by  Barr  (1967) 
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based  or,   a  large-scale  flume  study.  In  the  present  Investigation, 
values  of  FrA  were  computed  from  measured  data  in  the  range  of  3214  < 
ReA  <  3,830  and  are  compared  in  Chapter  5  with  those  obtained  by 
Keulegan  (1958). 

A  series  of  recent  studies  on  the  front  propagation  of  salinity- 
induced  gravity  currents  in  lock  exchange  flow  have  been  carried  out 
by  Huppert  and  Simpson  (1980),  Didden  and  Maxworthy  (1982),  Huppert 
(1982),  and  Rottman  and  Simpson  (1983).  According  to  their  studies, 
there  are  three  distinct  phases  which  may  exist  during  front  propagat- 
ion along  the  flume.  In  each  phase  the  front  movement  is  controlled 
by  different  combinations  of  dominant  forces  and  can  be  characterized 
by  a  specific  function  of  the  elapsed  time  t.  In  the  first  phase 
(initial  adjustment  phase),  during  which  the  initial  conditions  are 
important,  the  front  position,  xf,  changes  in  linear  proportion  to  the 
elaped  time,  t.  This  was  empirically  and- numerically  verified  by 
Rottman  and  Simpson  (1983).  The  second  phase  is  the  inertial  self-sim- 
ilar phase,  during  which  the  gravity  force  (or  buoyancy  force)  is 
balanced  by  the  inertia  force.  That  xf  moves  as  t2/3  was  theoretically 
derived  by  Huppert  and  Simpson  (1980),  and  was  experimentally  and 
numerically  confirmed  by  Rottman  and  Simpson  (1983).   At  the  stage  when 
the  viscous  force  becomes  more  important  than  the  inertia  force,  a 
viscous  self-similar  phase  is  reached.   In  this  third  phase  the  force 
balance  is  governed  by  gravity  and  viscous  forces.  That  xf  moves  as 
t1/5  was  obtained  by  Didden  and  Maxworthy  (1982)  by  estimating  the 
order  of  magnitude  of  forces  involved.   Also,  Huppert  (1982)  found  the 
same  characteristic  relationship  between  xf  and  t  for  the  viscous 
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aelf'-almilar  phase  after  staying  the  viscous  gravity  current  using  a 
lubrication  theory  approximation.  Within  the  first  two  phases  (when 
the  Reynolds  number  of  the  gravity  current  is  high),  the  front  movement 
is  locally  controlled  by  the  conditions  at  the  front.   In  tne  viscous 
phase  (]  >w  Reynolds  number  gravity  current),  the  front  shape  and  front 
speed  are  independent  of  the  conditions  at  the  front.  A  figure 
(Fig. 2. 3)  given  by  Rottman  and  Simpson  (1983)  clearly  illustrats  the 
existence  of  these  three  distinct  phases  as  the  salinity  front  propaga-- 
ted  along  a  flume.  In  the  figure,  x0  is  the  lock  length,  and  tg  is  a 
characteristic  time  parameter.  Evidence  of  these  three  phases  was 
found  from  the  laboratory  measurements  in  the  present  study,  and  a 
discussion  in  detail  is  given  in  Chapter  5. 

The  effect  of  the  internal  stratification  in  the  gravity  current 
front  on  the  front  movement  was  of  particular  interest  for  this  study 
because  a  strong  stratification  in  the  turbidity  front  can  exist  due  to 
the  settling  of  suspended  sediment.  Stratification  considerably 
reduces  front  speed,  according  to  the  findings  of  Kao  (1977). 

2.2  Stationary  Sediment  Wedge 
If  the  leading  wedge  (toe)  of  the  turbidity  current  has  propagated 
to  the  point  where  all  suspended  sediment  at  the  front  head  has 
deposited,  s  steady  or  quasi-steady  state  sediment  wedge  can  be 
established  in  a  closed-end  channel  where  the  force  balance  is  achieved 
by  a  raising  of  the  water  level  in  the  region  beyond  the  toe.  McDowell 
(1971)  mentioned  that  such  a  stationary  sediment  wedge  may  exist  in  an 
enclosed,  quiescent  basin  connected  to  open  water.  He  proposed  a 
simple  means  of  determining  the  sediment  fiux  from  outside  waters  into 
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Figure  2.3.  The  Front  Position  as  a  Function  of  Time  for  Different 
Releases  (after  Rottman  &  Simpson,  1983) 


the  basin,  as  the  product  of  ttie  mean  settling  velocity  of  parUoXtt 
and  the  arrested  wedge  length.   In  addition,  based  on  analytical 
considerations,  he  found  that  for  the  case  of  particles  with  small 
settling  velocity  and  high  suspension  concentration,  inflow  velocity  in 
the  lower  layer  would  be  approximately  equal  to  0.71uA,  where  u4  is  the 
densimetric  velocity.  Yet,  no  evidence  was  provided  to  verify  his 
findings. 

Gole  et  al.  (1973)  conducted  an  experimental  investigation  of  a 
stationary  sediment  wedge  in  a  two-dimensional  (i.e  relatively  narrow) 
closed'-end  channel  by  permitting  the  denser  silt-laden  saline  water  to 
intrude  into  silt-free  saline  water.  Under  the  equilibrium  condition, 
the  surface  outflow  velocity  of  silWree  water  at  the  entrance  of  the 
channel  was  measured  and  found  to  be  approximately  0.35uA,  which  is 
only  one-half  of  McDowell's  (1971)  finding.   Based  on  the  volumetric 
conservation  of  water,  the  magnitudes  of  volume  flux  in  lower  and  upper 
layers  should  be  identical.  Accordingly,  if  the  zero«velocity  interface 
(stagnation  point)  occurs  at  the  midfdepth.  it  can  be  assumed  that  the 
depth*averaged  velocities  in  both  layers  at  the  channel  entrance  are 
equal.  Onoe  the  inflow  velocity  is  known,  then  the  sediment  flux  into 
the  channel  can  easily  be  computed  based  on  the  inflow  velocity,  inflow 
depth,  and  sediment  concentration  at  the  entrance. 

One  of  the  objectives  In  the  present  study  was  to  find  out  the 
longitudinal  suspended  sediment  distribution,  sediment  flux  into  the 
closedt-end  channel,  and  the-  flow  regime  in  a  stationary  sediment  wedge 
by  means  of  experimental,  analytical,  and/or  mathematical  model 
approaches.  The  results  are  discussed  in  Chapter  5. 
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Z.j   Mathematical  Mode.; 
Fischer  v 1976)  attempted  to  simulate  density  currents  in  estuaries 
by  a  numeri  i]  nodel  and  found  that  numerical  instability  or  large 
numerical  errors  which  occur  in  solving  the  convective^dif fusion 
equation  will  be  encountered  if  a  common  numerical  method  (e.g. 
leapfrog  method)  is  used  to  solve  the  eqution  for  the  density  current 
with  a  fairly  sharp  density  discontlnity  (i.e.  at  the  density  front). 
A  higher-order  difference  approximation,  Hermite  interpolation  func- 
tion, was  suggested  by  Fischer  to  be  used  to  formulate  the  convec- 
tive*diffusion  equation.  A  conclusion  resulting  from  comparing  the 
solutions  obtained  by  a  semi-analytical  method,  the  higher'-order 
Hermite  approximation,  and  other  methods,  is  that  the  numerical  errors 
were  considerably  reduced  by  using  the  higher-order  Hermite  function 
approximation. 

Miles  (1977)  developed  a  multi-layered,  one-dimensional  model  to 
study  the  detailed  vertical  structure  of  the  flow  and  salinity  distri- 
butions in  estuaries.  Non-uniform  layer  thicknesses  can  be  deployed 
over  the  water  depth  to  give  thinner  layers  at  the  locations  where 
vertical  velocity  or  salinity  varies  rapidly.  In  the  equations  of 
motion  and  convective-diffusion  equation,  a  mixing  length  turbulence 
model  was  adopted  to  express  the  vertical  momentum  and  mass  transfer  in 
the  stratified  flow  with  stratified  buoyancy  effects.  A  semi-implicit 
numerical  scheme  was  used  to  formulate  the  governing  equations  in  which 
explicit  finite  differences  for  the  convective  terms  and  Implicit 
finite  differences  for  the  remaining  terms  were  used  in  the  equations. 
Subsequently,  the  double  sweep  method  (Abbott,  1979)  was  applied  to 
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3olve  the  finite  difference  formulations  with  known  boundary  cor.-ii- 
tions. 

Perrels  and  Karelse  (19H1   developed  a  two-dimensional,  laterally 
averaged  mode]  to  study  salinity  intrusion  in  estuaries.  A  finite 
difference  method  was  selected,  in  combination  with  a  coordinate 
transformation,  for  the  numerical  integration  of  the  system  of  differ- 
ential equations.  An  explicit  technique  was  used  in  the  longitudinal 
direction  and  an  implicit  technique  in  the  vertical  direction  for  both 
equations  of  momentum  and  convective-dif fusion.  For  the  continuity 
equation  a  central  difference  scheme  was  used.  The  calibration  and 
verification  of  the  model  were  carried  out  by  comparing  model  results 
with  experimental  data  obtained  from  laboratory  flume  tests. 

One  of  the  first  mud  transport  models  was  developed  by  Odd  and 
Owen  (1972).  It  is  a  twoi-layered,  oneMimensional  coupled  model 
which  simulates  both  the  tidal  flow  and  mud  transport  in  a  well-mixed 
estuary.   The  two  layers  can  be  of  unequal  thickness,  with  uniform 
properties  (e.g.  flow  velocity,  suspension  concentration)  assumed  for 
each  layer.  The  equations  of  motion  and  continuity  for  each  layer  are 
solved  using  a  finite-difference  formulation,  while  the  convective- 
diffusion  equation  governing  the  transport  of  suspended  sediments  in 
two  layers  is  solved  using  the  method  of  characteristics.  Erosion 
and  deposition  are  included  in  this  model. 

Ariath;irai  <19T>)   and  Ariathurai  and  Krone  (1976)  developed  an 
uncoupled  two-dimensional,  depth-averaged  cohesive  sediment  transport 
model  which  used  the  finite  element  method  to  solve  the  convective-dif- 
fusion  equ&tion.  The  model  simulates  erosion,  transport,  and  deposi- 
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tion  of  suspended  cohesive  sediments.   Aggregation  of  cotMSivt  '-lertlment 
was  accounted  for  by  determining  the  sediment  settling  velocity  as  a 
function  of  the  suspension  concentration.   Ariathurai  et  al.  (1977) 
modified  the  model  to  solve  the  two-dimensional,  laterally-averaged 
suspended  sediment  transport  problem.  The  model  was  verified  jsmg 
field  observations  from  the  Savannah  River  Estuary.  Required  data  for 
the  model  include  the  two-dimensional,  lateral ly«averaged  velocity 
field,  diffusion  coefficients,  and  sediment  deposition  and  erosion 
properties. 

Kuo  et  al.  (1978)  devloped  a  two-dimensional,  lateral] y-averaged, 
coupled  model  which  simulates  the  motion  of  water  and  suspended 
sediment  near  the  turbidity  maximum  of  an  estuary.  The  vertical  dimen- 
sion was  divided  into  a  number  of  layers,  and  a  finite  difference 
method  was  used  to  solve  the  equations  of  motion,  continuity  and 
convective^dif fusion  for  each  layer.  Stratification  effects  on  the 
vertical  momentum  and  mass  diffusion  coefficients  were  considered  in 
this  model.  Furthermore,  an  empirical  formula  was  used  for  the  rela*< 
tionship  between  the  vertical  momentum  or  mass  diffusion  coefficient 
and  the  local  Richardson  number,  suggested  by  Pritchard  (1960).  The 
longitudinal  momentum  and  mass  diffusion  coefficients  were  obtained  by 
multiplying  a  constant  ( =  1 0^ )  times  the  vertical  momentum  and  mass 
diffusion  coefficients,  respectively.  Erosion  and  deposition  were 
accounted  for  in  the  convective*-diffusion  equation  for  the  bottom 
layer. 
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2.4  uepositional  Iroperties  of  fine-Grained  Sediments 

Two  different  types  of  fine*gralnea  sediments,  cohesionless  and 
cohesive  sediments,  were  used  in  the  present  investigation  of  sedimen- 
tation in  closedfcend  channels.  Understanding  the  sedimentation 
characteristics  in  the  channel  requires  knowledge  pertaining  to  the 
settling  characteristics  of  suspended  sediment  and  the  depositional 
properties  of  sediment.  Several  important  depositional  properties  of 
fine-grained  sediments  are  reviewed  as  follows. 
2.4.1  Settling  Velocity 

For  cohesionless  fine  sediments,  the  settling  velocity  can 
practically  be  assumed  to  be  concentration* independent  and  can  easily 
be  determined  by  Stokes'  formula  (Daily  &  Harleman,  1966): 

2 

dm  (V13  8 

ws  ■  (  2.3  ) 

18  v 

where  ws  is  particle  settling  velocity,  dm  is  median  grain  size,  and  Gs 

is  the  specific  gravity  of  sediment. 

However,  for  cohesive  sediments  (particle  sizes  less  than  about  20 

microns),  the  settling  velocity  is  dependent  on  floe  or  particle  size, 
suspension  concentration,  local  physico-chemical  conditions,  and 
microbiological  activity  in  the  water  or  at  the  particle  surface.  Of 
these  factors,  tne  effect  of  suspension  concentration  on  the  settling 
velocity  war,  found  to  be  very  significant  due  to  the  change  of  the 
frequency  of  inte—particie  collision  in  suspension  with  changes  in 
concentration.   Collision,  which  Influences  the  rate  and  degree  of  floe 
aggregation,  is  caused  by  Brownian  motion,  the  presence  of  a  velocity 
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gradient,  and  differential  particle  nettling  velocities  (Hunt,  1980i. 
Krone  (1962)  ha3  discussed  the  effect  on  inter-partlele  collision  in 
suspension  by  each  mechanism.  It  wis  noted  that  when  settling  occurs 
under  static  or  quiescent  flow  conditions,  i.e.  the  conditions  of  the 
present  study,  Brownian  motion  and  differential  settling  velocities 
mainly  determine  the  frequency  of  collision  and,  consequently,  the 
degree  ot  aggregation. 

Mehta  (1986)  summarized  the  relationship  between  settling  velocity 
of  cohesive  sediment  and  Initial  suspension  concentration  in  three 
different  concentration  ranges.  At  very  low  concentrations,  the  rate 
of  aggregation  is  negligible,  and  the  settling  velocity,  ws,  does  not 
depend  on  the  suspension  concentration.  At  moderate  concentrations, 
aggregation  causes  ws  to  increase  with  concentration,  and  a  relationr 
ship  was  found,  from  laboratory  settling  column  analyses,  of  the  form 

"1 
ws  -  ki  C  (  2.n  • 

where  k-  depends  on  the  sediment  composition,  while  n^  was  theoretical** 
ly  and  experimentally  found  to  be  equal  to  1.33  by  Krone  (1962)  for  San 
Francisco  Bay  sediment,  and  to  be  less  than  unity  by  Teeter  (1983)  for 
Atchafalaya  Bay  sediment. 

At  high  concentrations,  w3  decreases  with  increasing  concentration 
because  ws  is  hindered  by  the  mutual  interference  of  particles  and  by 
the  upward  flux  of  fluid  escaping  through  the  small  spaces  among  the 
network  of  aggregates.  Brownian  motion  becomes  important  under  these 
conditions.  The  following  relationship  has  been  suggested  by  the  work 
of  Richardson  and  Zaki  (1951)  to  be  applicable  in  this  range: 
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n? 

wa  -  wsO  (I'-KoC)  (  2.5  ) 


Where  wa0  is  a  reference  settling  velocity,  kj  is  a  coefficient  which 
depends  on  the  sediment  composition,  and  n2  is  a  coefficient  which  has 
been  analytically  derived  to  be  equal  to  1.65  by  Richardson  and  Zaki 
(1951)  and  was  empirically  found  to  be  equal  to  5.0  by  Teeter  (1983). 

An  example  of  the  variation  of  the  settling  velocity,  ws,  with 
concentration  is  shown  in  Fig. 2.1,  which  is  based  on  measurements  in  a 
settling  column  using  mud  in  salt  water  from  the  Severn  Estuary, 
England  (Thorn,  1981). 

Settling  column  tests  for  the  settling  velocities  of  the  fine' 
grained  sediments  used  in  the  present  flume  study  were  carried  out  at 
the  Coastal  Engineering  Laboratory  of  the  University  of  Florida.  In 
these  tests,  the  settling  velocities  of  kaolinite  were  examined  over 
the  concentration  range  from  0.15  g/1  to  10.0  g/1,  and  the  results  are 
presented  in  Chapter  1. 
2.1.2  Deposition  Rate 

Krone  (1962)  and  Mehta  (1973)  conducted  laboratory  deposition 
experiments  in  which  they  monitored  the  variation  of  suspension 
concentration  with  time  under  a  given  applied  bed  shear  stress, 
T5.  In  several  tests  1-5  was  selected  to  be  smaller  than  the  critical 
deposition  shear  stress,  i0(j,  i.e.  the  shear  stress  below  which  all 
initially  suspended  sediment  deposits  eventually. 

With  reference  to  the  timenconcentration  curves  in  Fig. 2. 5  (Mehta, 
1973).  the  linear  portions  corresponding  to  the  range  of  low  concentre- 
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Figure  2.t.  Median  Settling  Velocity  vs.  Concentration  for 
Severn  Estuary  Mud  (after  Thorn,  1981). 


TIME  (mens) 


ipth   7b     l  3 


Figure  2.5.  Suspended  Sediment  Concentration  vs.  Time  for 
Three  Cohesive  Sediments  (after  Mehta,  1973). 
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tiorcs,  where  aggregation  is  negligibly  slow,  leads  to  the  relation 

dm  Tj 

-  -  WSC  (1 I         ,  Tb  <  T0d     (  2.6  ) 

dt  tej 

where  ra  is  the  mass  of  suspended  sediment  per  unit  bed  area  over  the 
depth  of  flow.  Eq.2.6  is  a  mass  balance  equation  which  essentially 
represents  dilution  of  the  suspension  with  time.  Eq.2.6  can  also  be 
used  to  compute  the  deposition  rate  of  suspended  particles. 

When  aggregation  is  proceeding  at  a  significant  rate,  as  in  the 
initial  phase  of  the  curves  of  Fig. 2. 5,  the  dilution  rate  of  suspension 
has  beer,  described  by  Krone  (1962)  as 


dm  it,      1 

■  *  [k3C  (1  N  )]  ,     xh   <  Tcd    (  2.7  ) 

dt  tc(j     t 


where  ko  is  a  factor  that  must  be  determined  empirically,  and  t  is  the 
time  since  the  beginning  of  aggregation. 

Mehta  (1973)  also  carried  out  extensive  deposition  tests  under 
steady  flows  using  kaolinite  as  well  as  muds  from  San  Francisco  Bay  and 
Maracaibo  Estuary,  Venezuela,  and  in  addition  reanalyzed  some  previous 
flume  data.  The  critical  deposition  shear  stress,  Tcd,  was  found  to 
depend  on  the  sediment  composition,  and  varied  from  O.Ci  to  0.15  N/m2. 
Note  that  t0(j  of  kaolinite  was  equal  to  0.15  N/m2  which  will  be 
referred  to  in  comparison  with  tc(j  of  kaolinite  obtained  from  the 
laboratory  measurements  of  the  present  study. 
2.1.3  Sorting  of  Deposited  Sediments 

Dixit  et  il.  (1982)  analyzed  a  series  of  data  from  deposition 
experiments  conducted  by  McAnally  (Mehta  et  al .  1982)  in  a  100  m  long 
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flume  at  the  U.S.  Army  Corp:-  of  Engineers  Waterway  Experiment  Station, 
Vicksburg,  Mississippi.  In  one  test,  a  kaolin! te-water  mixture  was 
introduced  from  the  headbay  of  the  flume  with  an  initial  concentration 
on  the  order  of  10  g/1.  The  depth  of  flow  was  0.17  m,  the  bed  shear 
stress  was  0.11  N/m2  ana  Manning's  roughness  coefficient  was  approxi- 
mately 0.01.  The  test  was  run  for  one  hour,  and  a  smooth  surface 
without  any  ripple-scaled  bed  features  was  found  at  the  end  of  the 
test.  The  variation  of  meaian  grain  size,  dm  of  the  dispersed 
sediment  from  the  deposit  with  distance  from  the  flume  headbay  is  shown 
as  Fig. 2. 6(a).  Comparing  these  values  with  dra  =  0.001  mm  of  the 
injected  sediment,  a  sorting  effect  is  evident.  In  Fig. 2. 6(b)  the 
settling  velocity  of  the  same  deposit  computed  from  the  mass  balance 
equation  is  plotted  against  the  distance  from  the  headbay.  Similar 
sorting  trends  of  dm  and  ws  with  distance,  x,  were  found  from  the 
deposits  of  the  present  laboratory  experiments,  and  are  presented  in 
Chapter  5. 

Furthermore,  the  flocculation  factor,  F,  of  cohesive  sediment  is 
defined  as  the  ratio  of  the  median  settling  velocity  of  the  floes  to 
the  settling  velocity  of  the  primary  constitutive  particles.  According 
to  the  relationship  between  F  and  the  median  primary  particle  diameter, 
dm,  found  by  Migniot  (1968)  and  confirmed  by  Dixit  et  al.  (1982),  the 
dependence  of  cohesive  sediment  aggregation  on  primary  particle  size 
can  be  demonstrated  in  a  quantitative  manner.   A  similar  relationship 
between  F  and  dm  was  found  in  the  present  study  and  is  discussed  in 
Chapter  5. 
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Figure  2.6.  (a)  Median  Dispersed  Grain  Size  vs.  Distance; 
(b)  Suspended  Sediment  Settling  Velocity  vs. 
Distance  (after  Dixit  et  al.  198?) 


CHAPTER  3 
METHODOLOGY 


3.1  Pimenslonal  Analysis 
3.1.1  Front  Shape 

The  height  of  front  nose.  The  geometry  of  the  experimental 
channel  system  is  given  in  Fig. 1.1  (a)  and  (b).  The  channel  has  the 
same  cross  section  from  the  entrance  to  the  closed  end.  The  bottom 
of  the  channel  system  is  in  the  same  horizontal  plane,  and  the  water 
depths  at  main  and  closed"end  channels  are  the  same  before  removing  the 
gate  located  at  the  connection  between  both  channels.  The  suspended 
sediment  source  at  the  main  channel  remains  constant.  Here  L,  B,  and  H 
denote  channel  length,  channel  width,  and  water  depth,  respectively. 
Let  pw  and  v  be  respectively  the  density  and  kinematic  viscosity  of 
fresh  water,  and  pw  ♦  Ap  and  v\    be  the  corresponding  quantities  for  the 
waterf-sediment  mixture.  Because  the  density  and  temperature  variations 
occurring  throughout  a  test  were  relatively  small,  the  average  value  of 
the  densities  can  be  expressed  as  pu,  and  of  the  kinematic  viscosities 
as  v. 

As  mentioned  in  Chapter  2,  the  front  shape  of  gravity  currents 
has  interested  many  investigators.  During  the  present  laboratory 
experiments,  a  typical  front  characteristic  of  turbidity  current  was 
observed  and  is  shown  In  Fig. J.I.  Let  h-|  and  h2  be  the  heights  of 
front  nose  and  front  head,  respectively  (Fig. 2. 2).  Then  an  expression 

31 


^^^^^»    r 

Figure  3.1.  F>ont;  Shape  cf  a  Propagating  Turbidity  Current. 
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for  h^  la  given  as 

hj  -  hi  (h2,  K,  Uf,  v)  (  "?.1  ) 

where  Uj  is  the  front  speed.   An  operation  using  Buckingham  -  it  theorem 
(Streeter  &  Wylie,  1975)  leads  to  the  desired  form: 


h^       Ufh2    h2 

—  -  f  (  ,  —  )  (  3.2  ) 

h2       v     H 


which  shows  that  the  ratio  of  nose  height  to  head  height  may  be  depend^ 
ent  on  the  local  head  Reynolds  number,  Ufh2/v,  and  fraction  aepoh, 
h2/H. 

The  height  of  front  head.  One  significant  feature  of  the  front 
shape  is  that  the  height  of  the  front  head  is  larger  than  the  height  of 
front  neck,  and  an  intensive  mixing  takes  place  at  the  interfa;e  right 
behind  the  front  head  (Fig. 2. 2). 

A  series  of  photographs  of  the  front  shape  taken  in  the  present 
experiments  indicated  that  the  size  of  the  advancing  front  and  the 
degree  of  front  mixing  decrease  as  the  front  speed  decreases.  The 
front  head  completely  disappears  after  the  front  travels  for  a  relati* 
vely  long  distance.  Let  tta  denote  the  height  of  front  neck,  then  the 
expression  of  the  height  of  the  front  head  is 

h2  -  n2  (H,  hj,  ur,  v,  ws,  B)  (  3.3  '; 

where  w3  is  the  particle  settling  velocity  at  the  front.  It  leads  to 
the  following  expression: 


h2      Ufh^   03   ws    B 
113        \i     H    Uf    H 


_ 
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"hf  quantity  t^/h-^  w^s  found  to  be  independent  of  B/H  By  Keulegan 
(1958),  but  iependenl  upon  the  '"faction  depth,  ha/H,  by  Simpsor  and 
Britter  (1979).  A  fina]  form  of  hj/h-^  can  be  written 


h2      ufho  ho       w3 

—  -  f  < ,  —  ,  —  )  (  3.5  ) 

hj        v     H    Uf 


which  indicates  that  besides  the  fraction  depth  l^/H,  h2/hj  may  be 
dependent  on  the  local  neck  Reynolds  number,  Ufho/v  ,  and  the  relative 
settling  velocity  of  suspended  sediment,  ws/Uf. 
3.1.2  Front  Speed 

Initial  front  speed.  Let  Ufi  denote  the  initial  front  speed  of 
turbidity  front,  that  is  the  velocity  immediately  after  the  gate  is 
opened.  JdeaJly  the  opening  is  very  smooth  and  sudden.  Within  such  a 
short  time,  the  effect  of  the  sediment  settling  velocity  on  the  front 
movement  is  negligible.  The  expression  of  Uf|  with  the  related 
geometrical  and  physical  quantities  is  given  as 

Ufi  -  Ufi  (  g,  Api  ,  pw,  v,  H,  B  )       (  3.6  ) 

where  Ap-|  is  the  density  difference  between  the  main  channel  and 
closed^end  channel  waters.  One  important  velocity  characteristic, 
densimetric  velocity  u^  ,  was  introduced  in  eq.2.2.   Another  signifi- 
cant dimensior.less  parameter  wa3  noted  as  the  densimetric  Reynolas 
number,  ReA  ,  by  Keulegan  (1957): 


uA  H 
Re.  ■ (  3.7 


3b 
A  final  form  for  eq.  (3.6)  Is 


uf1       uAH   B 

-  f  (  ,  -  )  (  3.8  ) 


which  suggests  that  the  dimenslonless  initial  front  velocity,  Ufj/u^, 
may  depend  upon  the  densimetric  Reynolds  number,  u4H/«,  and  the 
wldth'-to-depth  ratio,  B/K. 

Instantaneous  front  speed.   A  major  concern  with  front  behavior 
of  advancing  turbidity  current  in  the  present  study  was  how  these 
physical  parameters  (e.g.  channel  geometry,  boundary  suspension 
concentration,  and  sediment  properties)  influence  the  front  speed.  It 
will  be  discussed  in  this  section  from  two  viewpoints;  one  is  by 
relating  the  front  speed  to  the  channel  geometry  and  suspension 
concentration  at  the  entrance,  and  the  other  is  by  relating  the  front 
speed  to  the  local  front  shape  and  suspension  concentration  at  the 
front. 

The  front  speed,  Uf,  and  these  possibly  influencing  parameters  are 
grouped  as  follows: 

Uf  ■  Uf  (x,  H,  B,  g,  Api ,  pw,  v,  ws  )    (  3.9  ) 

This  equation  can  be  reduced  to  be  in  a  final  form 


Uf      x  B  uAH  w3 

—  -  f  (  -,  -,  .  —  )  (  3. 10 

uA       H  H   v   uA 


which  shows  that  the  dimensionless  front  speed  may  be  a  function  of 
front  position,  width^tc-depth  ratio,  densimetric  Reynolds  number,  and 
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relative  settling  velocity.     Eq.3.1G   la  similar  to  the  dimensional 
analysis  result  of  front  speed  of  a  sa.' ine  wedge,   which   is  advancing 
against   the  river  flow,    obtained  by   Keuiegan   (1971): 


Uf  x     B     u^H     ur 

—  -   f   (  -,  -,  ,   —  )  (    3.11    ) 

uA  H     H        v       uA 


except  for  the  term  ur/uA,  where  ur  is  the  velocity  of  the  river 
opposing  the  advancing  wedge. 

Barr  (1967)  investigated  lock  exchange  flow  of  saline  water  in  a 
large-scale  flume  (i.e.  flume  width,  B,  much  greater  than  water  depth, 
H),  and  found  that  the  dimensionless  front  speed,  Uf/u&,  is  a  function 
of  the  dimensionless  distance,  x/H,  and  the  densimetric  Reynolds 
number,  Re*: 


Uf       x  uAH 

—  -  f  (  -,  )  (  3.12  ) 

uA       H   v 


For  a  small-scale  flume,  it  would  be  expected  that  Uf/uA  would  also  be 
a  function  of  B/H.  Therefore,  the  dimensionless  grouping  of  instanta- 
neous front  speed  for  the  lock  exchange  flow  of  saline  water  would  be 
the  same  as  eq.3.10  except  for  the  relative  settling  velocity  term, 
ws/uA. 

It  was  noted  in  the  experiments  during  the  front  propagating  along 
the  channel  that  whenever  a  decrease  in  the  front  speed  was  observed,  a 
decrease  in  the  front  height  was  also  observed.  This  suggested  that 
possibly  a  relation  between  the  front  speed  and  the.  front  height 
existed.  The  local  density  of  the  front,  denoted  as  pw+Ap' ,  decreases 
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by  dropping  a  portion  of  suspended  sediments  to  the  bottom  as  it 
propagates  downstream.  Thus  the  effect  of  density  decrease  of  the 
front  on  the  instantaneous  front  speed  must  be  taken  into  account.  The 
expression  of  the  front  speed  can  be  written  as 

Uf  -  Uf  (  113,  g,  Ap',  pw,  v,  ws,  H,  B  )       (  3.13  ) 

here  the  neck  height,  I13,  is  used  as  a  height  characteristic  of  the 
front.  Equation  3.13  can  be  reduced  to 


Uf  Ufh^   ws  ho  B 

f  (  -  .—,—,-)     (  3.11  ) 


"  A   T72 

Up'gh3/pw]  v     uf  H   H 

where  Fr&  designates  local  densimetric  Froude  number.  Here,  the 
parameter  of  ws/uf  will  Introduce  a  stratification  effect  of  suspension 
concentration  at  the  front  on  the  dimensionless  front  speed,  as  found 
by  Kao  (1977).  Besides  the  internal  stratification  in  the  front,  the 
local  densimetric  Froude  number,  FrA,  may  also  be  influenced  by  the 
local  neck  Reynolds  number,  Ufhj/v,  the  fraction  depth,  hj/H,  and  the 
width-to-depth  ratio,  B/H. 
3.1 .3  Front  Position 

The  identification  of  three  phases,  each  controlled  by  specific 
dominant  forces,  in  the  front  propagation  along  a  flume  were  mentioned 
in  section  2.1.3.  In  each  phase  the  front  position,  Xf,  was  found  to 
be  a  specific  function  of  the  elapsed  time,  t.  In  the  present  study  an 
attempt  using  dimensional  analysis  was  made  to  provide  the  proper 
dimensionless  forms  of  Xf  and  t  for  analyzing  measured  data  and 
presenting  the  results.  An  expression  of  the  front  position  may  be 
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stated  as 

*f  -  *t   <  B.  H,  B,  t,  Api,  pw,  v,  ws       (  3.1b  ) 

An  operation  using  Buckingham  h  i  theorem  leads  eq.3.1">  to  a  final  form 
with  relevant  dimensionless  parameters: 


Xf       B   lut   ws  uAH 

—  -  f  (  -  ,  ,  — ,  )  (  3-16  ) 

H        H    H    u&   \) 


which  indicates  that  the  dimensionless  front  position  is  dependent  on 
the  width-tcdepth  ratio,  B/H  ,  the  relative  settling  velocity,  Wg/u^, 
the  densimetric  Reynolds  number,  u^H/v,  and  a  dimensionless  elapsed 
time,  uAt/H. 

Barr's  (1967)  study  on  the  lock  exchange  flow  of  saline  water  in  a 
large-scale  flume  provided  a  group  of  dimensionless  influencing 
parameters  for  the  dimensionless  front  position,  Xf/H,  i.e. 

Xf       u&t   u^H 

—  =  f  (  ,  )  I   3.17  ) 

K        H    v 

Comparing  eq.3.16  with  eq.3.17,  it  can  be  seen  that  two  additional 

dimensionless  factors  are  included  in  eq.3.16.  They  are:  1)  the 

widtlr-to-dept>i  ratio,  B/H,  which  will  introduce  an  effect  on  the  front 

position  Xf/I!  when  the  gravity  current  propagates  in  a  small-scale 

flume,  and  ?)  the  relative  settling  velocity,  ws/uA,  which  is  a 

particular  characteristic  of  turbidity  currents. 

3.1,1  Deposition  Rate 

One  of  the  major  concerns  in  the  present  study  was  the  deposition 

of  fine-grained  sediments  in  a  closed^end  channel.  The  sediment 
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deposition  rate  in  each  subsection  of  the  channel  bed  can  be  computed 
using  the  deposited  sediment  mass  divided  by  the  effective  duration 
(i.e.,  the  time  elapsed  from  ore  passing  of  the  turbidity  front  to  the 
end  of  the  test)  and  depositt-d  area.  The  local  deposition  rate  would 
be  expected  to  vary  with  distance  due  to  the  distance  variations  of 
mean  suspension  concentration  and  of  settling  velocity  of  suspended 
sediment.  The  variation  of  deposition  rate  with  distance  and  other 
relevant  parameters  are  expressed  as 

6  -  6    (  x,  H,  ws,  g,  B,  pw,  4P1)        (  3.18  ) 

where  6  is  the  local  deposition  rate.  Eq.(3.18)  can  be  further  reduced 
to  a  desired  form  as 


5         x    B    ws    Api 

—  -  f  ( . , ,  )        (  3.19  ) 

5,        H    H    uA     pw 


where  Si  is  the  deposition  rate  at  the  entrance  of  the  closed»-end 
channel.  The  normalized  deposition  rate,  6/6],  is  thus  observed  to  be 
a  function  of  x/H  and  B/H.  Also,  that  6/ 6-]  depends  upon  the  relative 
settling  velocity,  ws/u&,  is  obvious.  The  influence  of  Api/pw  on  5/6-j 
can  be  made  apparent  by  introducing  different  sizes  of  the  flocculated 
particles  through  the  entrance. 

The  dimensionles3  parameters  obtained  in  eq.3.19  have  also  been 
used  to  examine  the  experimental  results  on  the  local  sediment  settling 
velocity,  ws,  and  the  median  dispersed  particle  size,  dm,  since  all 
data  on  6,  ws,  and  dm  were  obtained  by  analyzing  deposited  samples  from 
the  flume  tests. 


3.2  Devel opement  of  Mathematical  Model 
3.2.1  Moriel  Description 

Basic  governing  equations.  The  mathematical  model  developed  Is  a 
time  varying,  two-dimensional ,  coupled  finite  difference  model  that  is 
capable  of  predicting  the  vertical  and  temporal  variations  in  the 
suspension  concentration  of  fine-grained  sediments  and  in  flow  veloci- 
ties in  a  coastal  closed'-enj  waterway  with  a  horizontal  bottom.  In 
addition,  It  can  be  used  to  predict  the  steady-state  or  unsteady 
transport  of  any  conservative  substance  or  nonf-conservative  consti^ 
tuent,  if  the  reaction  rates  are  known. 

The  governing  equations,  which  control  the  flow  regime  and 
sediment  transport,  are  expressed  in  conservation  forms  (Roaehe,  1972) 
and  shown  a3 


3u   3w 

—  4  —  =  0  (  3.20  ) 

3x   3z 


3u   3u2   3(uw)     1  3p   3      3u 
3t    3x     3z       p  3x   3x     3x 


3      3u 

—  (ez  — ) 
3z     3z 


(  3.21  ) 


1  3p 

0  -  -g (  3.22  ) 

p  3z 


3C   3(uC)    3  3      3C    3      3C 

—  ♦ +  —  f(w-ws)C]  =  —  (ex  — )  +  —  (ez  — ) 

It    3x     3z  3x     3x    3z     3z 

-  d  +  r  (  3-23  ) 


m 

u,  M  -  velocity  components  in  the  x  and  I   directions, 
respectively 

p  =  pressure  force 

ex,  ez  ■  momentum  diffusion  coefficients  in  the  x  and  z  directions, 
respectively 

ex,  ez  -  nass  diffusion  coefficients  in  the  x  and  z  directions, 
respectively 

d  ■  sediment  deposition  rate 

r  =  sediment  erosion  rate 

In  fact,  the  physical  quantities  involved  in  these  equations,  such 
as  u,  w,  p,  and  C,  are  time-averaged  quantities  over  a  period  which  is 
greater  than  the  time  scale  of  turbulence  fluctuation  when  the  flow 
is  turbulent. 

Equation  3.20  is  the  continuity  equation  for  an  incompressible 
fluid.  Equation  3.21  is  the  equation  of  motion  for  an  incompressible 
fluid,  and  represents  the  longitudinal  momentum  conservation  of  the 
flow.  Equation  3.22  is  the  hydrostatic  equation  which  results  when  the 
vertical  components  of  the  flow  velocity  and  acceleration  are  smaller 
relative  to  the  horizontal  flow  velocity  and  acceleration.  Equation 
3.23  is  the  convective^diffusion  equation  for  suspended  sediment,  with 
resuspension  and  deposition  as  source  and  sink,  respectively. 

To  obtain  the  time-varying  solutions  of  the  longitudinal  and 
vertical  velocity  field,  eq.3.21  must  be  solved  with  the  continuity 
equation  (eq.3.20).  With  the  velocity  field  solved,  it  may  be  substi- 
tuted into  eq.3.23  to  solve  for  the  time^varying  concentration  field  of 
suspended  sediment. 
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Vertical  integration  of  basin  equations.  Since  the  physical 
quantities  in  the  channel  can  change  rapidly  over  a  short  vertical 
distance,  they  require  a  grid  size  that  is  much  smaller  in  the  vertical 
direction  than  in  the  longitudinal  direction.  The  fluid  motion  will  be 
considered  in  horizontal  slice3  with  an  exchange  of  mass  and  momentum 
between  these  slices.  A  sof-called  sparse  grid  system  (Miles,  1977)  in 
space  used  in  the  model  and  the  location  of  physical  quantities  within 
the  grid  are  shown  in  Fig.A.1  of  Appendix  A.  At  the  top  layer,  the 
free  surface  is  allowed  to  occupy  any  position  within  the  layer,  but 
the  thicknesses  of  the  others  are  fixed. 

Integration  over  the  height  of  the  j  th  layer  can  be  performed  by 
assuming  that  all  variables  are  practically  constant  through  the  depth 
of  any  layer,  and  that  the  fluxes  of  momentum  and  mass  normal  to  the 
bottom  of  the  channel  and  to  the  surface  are  zero. 

Employing  Leibnitz's  rule  for  the  vertical  integration  of  equa- 
tions 3.20,  3.21,  and  3.23  over  the  j  th  layer  and/or  the  surface 
layer,  results  in  the  following  equations: 

3h       3 

—  -  wN (uNAzN)  (  3.214  ) 

3t        3x 


3 
WT  "  wb  *  —  (UjAZj)  (  3.25 


3x 


3uj    3    2  1                  1    3p 

♦  —  (uj)  i  (wTuT  -  wbub) (— )j 

«    3x  Azj                 p,   3x 
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where: 

Uj ,  Azj,  pj  =  longitudinal  velocity,  height,  and  density  for  the 
j  th  layer,  respectively 

uj,  wT,  if  ■  longitudinal  velocity,  vertical  velocity,  and  shear 
stress  at  the  top  of  a  layer,  respectively 

ub>  wt>'  Tb  *  longitudinal  velocity,  vertical  velocity,  and  shear 
stress  at  the  bottom  of  a  layer,  respectively 

Ujj,  wk,  Az(j  =  longitudinal  velocity,  vertical  velocity,  and  height 
for  the  top  layer,  respectively 

h  =  water  surface  elevation. 

Equations. 3. 24  and  3.25  are  the  continuity  equations  for  the  top 
layer  and  all  other  layers,  respectively.  Equation  3.26  is  the 
longitudinal  equation  of  momentum,  and  eq.3.27  is  the  conveetive'-diffuf- 
sion  equation  for  suspended  sediment. 

Further  derivation  of  the  pressure  term,  ( 1 / p j ) ( 3P/3x) i  in  eq.3.25 
is  necessary.   The  hydrostatic  pressure  at  an  elevation  z  in  the  fluid 
can  be  obtained  from  eq.3.22 


p(z)  -g  J  p  dz  (  3.28 


and  its  longitudinal  gradient  is 
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where  pn  is  the  aensity  at  water  surface. 

The  mean  pressure  gradient  in  the  top  layer  can  be  derived  thus 
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where  pN  is  the  average  density  in  the  top  layer.  For  two  adjacent 
layers,  j  and  j+1,  the  relationship  between  the  layeri-averaged  pressure 
gradients  can  be  derived  from  eq.3.29,  and  the  final  expression  is 

3p        3P        g        3pj+i       3pj 

(  —  )j  -  {  —  )J  +  i  ♦  -  (Azj  +  1  ♦  Azj  )   (  3.31  ) 

3x        3x        2         3x         3x 

where  all  of  PJ  ,  (3p/3x)j,  Pj+1,  and  (3p/3x)j+1  are  layer-averaged 
quantities.  From  eq.3.30  the  mean  pressure  gradient  at  the  top  layer 
can  be  computed  with  a  known  mean  density  gradient  at  the  same  layer 
and  the  slope  of  water  surface.   Once  the  mean  pressure  gradient  at  the 
top  layer  is  solved,  the  mean  pressure  gradient  at  the  other  layers  can 
be  calculated,  based  on  eq.3.31. 

Boundary  conditions.  At  free  surface  boundary,  no  flow  will  occur 
across  the  free  surface,  and  the  vertical  velocity  on  tne  water  surface 
must  be  equal  to  zero.   The  shear  stress  induced  by  wind  acting  on  the 


water  surface  is  not  considered  in  this  study;  thus  the  shear  stress  at 
surface  boundary  can  be  regarded  as  zero.  No  mass  flux  of  suspended 
sediment  can  flow  through  the  free  surface.  This  condition  is  express*- 
ed  as 


3C 

—  =0  at  z  =  h        (  3.32  1 

3z 


At  the  bottom  boundary,  no  flow  car.  pass  the  bottom;  therefore, 
the  vertical  velocity  at  the  channel  bed  must  be  zero.  The  relation- 
ship between  the  friction  factor,  t0,    and  Manning's  roughness  coeffi- 
cient, n,  can  be  obtained  from  the  formulas  of  Manning  and  of  Darcy- 
Weisbaeh  (Daily  &  Harleman,  1966)  as 


Rl/'3 


f  3.33  ) 


where  R  is  hydraulic  radius.  Accordingly,  the  bottom  shear  stress,  tb, 
can  be  expressed  as 


Tb  '  Pb  ublubl  (  3.34  ) 

Azb 

where  pD,  ub,  and  Azb  are  the  density,  horizontal  velocity,  and 
thickness  of  the  bottom  layer,  respectively.   Like  the  free  surface 
boundary,  no  mass  flux  of  suspended  sediment  is  allowed  to  pass  through 
the  bottom  boundary,  that  is 
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in  addition,  iwc  boundary  conditions  at  the  entrance  of  the 

closed-end  channel  are  necessary;  I)  concentration  profile  of  suspended 
sediment,  and  2)  the  water  surface  elevation.  The  data  from  either 
condition  can  be  time-varying  quantities.  To  investigate  the  problem 
under  tide  effects,  the  surface  elevation  at  the  entrance  of  tne 
closed-end  channel  can  be  specified  with  a  periodic  function  including 
prototype  tidal  amplitude  and  period.  The  concentration  profile  of 
suspended  sediment  and  water  surface  elevation  involved  in  the  present 
study  did  not  vary  with  time.  Also,  the  boundary  conditions  at  the  end 
of  the  closed-end  channel  include  zero  horizontal  velocity  and  zero 
horizontal  concentration  gradient  at  the  vertical  wall. 

Momentum  and  mass  diffusion  coefficients.   In  turbulent  flow  of  a 
homogeneous  fluid  there  are  no  buoyancy  effects  and  lumps  of  fluid 
moved  by  turbulence  fluctuations  have  no  restoring  force  to  return  to 
their  original  position.  When  the  flow  is  stratified,  the  buoyancy 
effect  tends  to  restore  the  moved  lumps  back  to  their  original  posi- 
tion, and  results  in  the  reduction  of  the  turbulent  transfer  of 
momentum  and  mass.  Bowden  and  Hamilton  (1975),  among  others,  have 
pointed  out  the  necessity  of  considering  the  buoyancy  effects  on  the 
momentum  and  mass  transfer  in  stratified  flows.  In  order  to  consider 
the  buoyancy  effect  of  stratification,  the  momentum  and  mass  diffusion 
coefficients  were  related  to  a  dimensionless  parameter,  which  is  a 
function  of  vertical  density  and  velocity  gradients.  Furthermore, 
according  to  many  investigations  based  on  measurements,  a  typical 
relationship  between  either  vertical  momentum  diffusion  coefficient, 
ez,  or  mass  diffusion  coefficient,  cz,   and  the  dimensionless  parameter, 
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Rj,  is  given  as 


11 
sz  =  e0  (  1  ♦  aiRi  )  (  3-36  ) 


12 
e0  (  1  -  a2Ri  )  (  3.37  ) 


where  e0  and  e0  are  the  vertical  momentum  and  mass  diffusion  coeffici- 
ents in  a  homogeneous  flow,  respectively;  c^  ,  Qi  ,  u£,  and  q^  are 
constants  needed  to  be  calibrated  using  data;  and  Bj  is  the  gradient 
Richardson  number  defined  as  follows: 


3p 

g   3z 

Ri  -  '■ (  3.38  ) 

p   3u  2 
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which  is  generally  used  as  a  stability  index  in  stratified  flows 
(Turner,  1979  K  For  the  vertical  diffusion  coefficients,  e0  and  aj  are 
assumed  identical  to  e0  and  02,  respectively,  while  q^  is  taken  to  be 
-1/2. 

One  well-known  formula  for  the  vertical  mass  diffusion  coefficient 
has  been  suggested  by  Pritchard  (I960).  The  formula  was  developed  from 
a  study  of  velocity  and  salinity  distributions  in  the  James  River 
estuary  by  fitting  observational,  results.  The  formulas  of  vertical 
mass  diffusion  coefficient  in  the  stratified  and  homogeneous  flow 
fields  are  given  respectively  as  follows: 


ne 

- 

ez  -  e0  <  ■  ♦  C.276  Rj  )  (  3.39  ) 


-3    iujz?(H-z)2 

t0  -  8.59  x  10 (  i.W    ) 

H3 


where  H  is  the  total  water  depth  and  z  is  the  elevation  at  which  >7  is 
being  calculated.   In  eq.3.39,  ap  and  q2  are  equal  to  0.276  and  -?, 
respectively. 

For  the  longitudinal  momentum  diffusion  coefficient,  ex,  Festa  and 
Hansen  0976)  changed  the  value  of  ex  from  ex  =  ez  to  ex  =  Mr  ez  with 
negligible  effects  on  the  results  of  their  tidal  model.  This  indicates 
that  the  results  are  not  sensitive  to  the  actual  value  of  ex.  Yet  they 
found  that  varying  the  longitudinal  mass  diffusion  coefficient,  ex, 
from  ex  =  ez  to  ex  =  10?  ez  did  produce  significant  changes  in  their 
model  results.  Dyer  (1973)  and  Kuo  et  al.  (1978)  suggested  that 
longitudinal  momentum  and  mass  diffusion  coefficients  are  on  the  order 
10J  of  the  corresponding  vertical  diffusion  coefficients. 

Sediment  settling  velocity  and  deposition  rate.   In  a  mathematical 
model  of  sediment  transport,  model  prediction  is  very  sensitive  to  the 
settling  velocity  of  suspended  sediment  given  as  input  data.  Thus  the 
determination  of  sediment  settling  velocities  for  a  mathematical  model 
requires  much  care. 

There  are  essentially  four  types  of  field  or  laboratory  measure- 
ment of  the  suspension  settling  velocity  mentioned  by  Mehla  (1986):  1) 
the  use  of  an  in  situ  fube  for  prototype  measurement,  2)    fitting 
analytical  solution  to  measured  suspended  sediment-depth  profile  from 
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the  prototype  measurement,  3)  use  of  laboratory  settling  column,  and  '•>) 
measurement  of  rates  of  deposition  in  the  flume. 

The  measurement  by  means  of  a  laboratory  settling  column  for  the 
settling  velocities  of  fine-grained  sediments  used  in  the  tests  was 
made  in  the  present  study.  The  principle  of  this  test  is  essentially 
based  on  the  known  relationship  between  the  downward  settling  flux  of 
suspended  sediment  and  the  dilution  rate  of  suspension  concentration  in 
the  settling  column.  During  the  test,  the  suspension  samples  are 
withdrawn  at  designated  times  from  several  vertically  aligned  taps, 
which  are  attached  tc  the  side  wall  of  the  settling  column.  The 
experimental  facility  and  procedure  of  the  settling  column  test 
performed  in  the  present  study  were  designed  with  reference  to  the  work 
by  Interagency  Committee  (1943),  Mclaughlin  (1959),  Owen  (1976),  and 
Christodoulou  et  ai.  (1971).     The  results  of  the  median  settling 
velocity  versus  initial  suspension  concentration  for  each  raw  sediment 
are  presented  in  Chapter  H .      However,  these  settling  velocities  cannot 
appropriately  represent  the  local  suspension  settling  velocities  in  the 
closed-end  channel  due  to  the  losses  of  a  part  of  raw  sediments  before 
they  reach  the  closeCend  channel  entrance  and  the  non*-uniform  flow 
regime  in  the  closed-end  channel. 

The  suspension  settling  velocity  in  the  closed-end  channel,  which 
is  subject  to  local  flow  and  concentration  conditions,  was  computed 
based  on  the  measurement  of  deposition  rate  from  each  test.   These 
computed  settling  velocities  are  expected  to  be  more  realistic  than 
those  obtained  using  settling  column  analysis  for  representing  the 
local  settling  velocities  of  suspension  in  the  closed-end  channel. 
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Therefore,  the  computed  settling  velocity  was  adopted  as  Input  data  in 
the  mathematical  model. 

In  the  present  study,  the  bed  shear  stress  generated  by  turbidity 
current  was  mucn  smaller  than  the  critical  erosion  shear  stress  which 
is  on  the  order  of  0.2  -  0.6  N/m?  (Parchure  and  Mehta,  1985).  Thus, 
the  term  for  resuspension  rate  of  bed  material  in  the  convective-diffu- 
sion  equation  was  omitted. 

The  mass  balance  equation  (eq.2.6)  was  used  for  computing  the 
deposition  rate  of  suspended  sediment  in  the  mathematical  model.  In 
this  equation,  the  critical  deposition  shear  stress,  tc(j,  must  be 
known  before  any  computation.  Therefore,  an  attempt  at  finding  tccj  for 
kaolinite  was  carried  out  based  on  the  measured  data  in  the  main 
channel  and  is  discussed  in  detail  in  Chapter  5.  Note  that  the  term 
for  sediment  deposition  rate  in  the  convective<-diffusion  equation  is 
only  considered  at  the  bottom  layer  of  the  numerical  grid  system. 
3.2.2  Finite  Difference  Formulation 

There  are  many  ways  to  present  the  derivative  terms  of  the 
governing  equations  in  the  finite  difference  formulation  for  solving 
the  flow  velocity  and  sediment  concentration  numerically.  An  explicit 
numerical  scheme  was  used  to  develop  the  mathematical  model.   In  order 
to  gain  better  stability  and  accuracy  of  the  solutions  of  the  finite 
difference  formulations,  the  temporal,  convective,  and  diffusive  terms 
in  the  equations  of  motion  and  convective'-diffusion  were  time^  and 
space-centered. 

With  reference  to  the  grid  system  of  Fig.A.1  of  the  Appendix, 
three  independent  variables,  x,  ?.,   and  t  denote  longitudinal  and 
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vertioa]  space  coord ina let.  Mid  time,  respectively.  Using  i  and 
represent  the  number  of  apace  Intervals  in  the  x  and  z  directions, 
respectively,  and  using  n  to  denote  the  number  of  time  intervals  that 
have  elapsed,  variables  will  be  represented  using  i,  j,  n  subscripts, 
where  i,  j,  or  n  ■  0,  1,  2,  3,  etc.  The  surface  elevation,  h,  is  a 
function  of  x  and  t  only,  while  the  layer  thickness,  Az,  is  a  function 
of  t  for  the  surface  layer  and  a  function  of  z  only  for  the  other 
layers. 

The  finite  difference  approximations  of  eqs.3.24,  3.25,  3-26,  and 
3.27  are: 
the  free  surface  elevation 

n+1    n-1  n        n  n       n 

ni   *  hj      n  1  Azi  +  itN  t  AZj  [j   n    Azi.N  +  izi-1 ,N  B 

■  w  *  —  [ u   «  u     ] 
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the  vertical  velocity 
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equation  of  motion 
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n  n  n         n 


C(ui,j+ui,j  +  V'W,,j*i+wiJ.1  ,j4l  )    -   (ulf  j.-i-Ui^J'wi^+Win  _j)] 


n 


n-1        H*1  n-1  n-1  n-1  n-1        n"1      n"1 

Hex        +ex  )(u  -u        )    -    (ex  +ex        )(u        -u  )] 

i.J        1  +  1.-J  I+1.J      i,J  l»1,J        i.j        i,j      i  - 1 .  j 

2Ax2 


8*1             n-1                  n-1           n-1                  n-1         n-1             n-1      n-1 
(ez            +ez                 )(u            -u        )         (ez        *ez             )(u        -u            ) 
1              i,j  +  1        i+l,j  +  1        i,j  +  1      i,j                 i,j        i  +  l,j        i,j      i.jH 
-[ ] 


n-1  n-1  n-1  n-1  n-1 

Azj  (4zj        ♦    AzJ  +  1)  (aZj        +   Azj-,) 

3P      2 

2   (  -  h,j 
3x 


n  n 

(Pi.j    +   Pi  +  1 ,  i > 


(    3.13    ) 


and  the  convective-dif fusion  equation  for  sedierant 
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In  addition,  the  finite  difference  formulation  of  the  gradient 
Richardson  number,  eq.3.38,  can  be  expressed  as 
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(  3.15  ) 
The  procedure  of  numerical  calculation  starts  with  all  the 
variables  assigned  at  their  initial  values,  and  moves  from  the  channel 
entrance  to  the  closed  end  in  x  direction,  from  the  tup  layer  to  bottom 
layer  in  z  direction.  With  all  variables  known  at  the  n  th  time  step, 
the  continuity  equation  for  the  top  layer  (eq.3.41)  is  used  to  calcula- 
te hj   .  With  hj   known,  i)j  <  and  Cji  can  be  calculated  for  all  i's 
and  j's  using  equations  3-^3  and  3-tf.  Knowing  all  the  u's  for  the 
(n+1)st  time  step  allows  all  the  w's  for  the  (n+1)st  time  step  to 
be  calculated  using  eq.3.^1. 

Next,  the  density,  pressure  gradient,  and  longitudinal  and 
vertical  mass  diffusion  coefficients  are  calculated  for  the  (n+1':st 
time  step.   ~ubst-qu<-ntly ,  the  entire  procedure  described  is  repeated  to 
calculate  the  values  of  the  variables  at  the  (n+?)nd  time  step.  A 
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brief  flow  chart  which  summaries  the  sequence  of  numerical  calculations 
is  shown  In  Appendix  A.?. 
3.2.3  Stability  and  Accuracy 

As  the  derivatives  in  the  governing  equations  are  centered  In  time 
and  apace,  the  finite  difference  scheme  has  second  order  accuracy. 
Unconditional  stability  does  not  exist  because  explicit  finite  differ- 
ence formulations  are  used.  A  stability  analysis  of  the  full  set  of 
governing  equations  is  impossible  due  to  the  involvements  of  nonli- 
near terms  and  of  the  coupling  characteristic.  Instead,  parts  of  the 
finite  difference  formulation  representing  different  physical  process- 
es have  been  analyzed  separately.  The  stability  condition  of  each  part 
is  essential  for  the  overall  stabitity  of  the  model,  but  is  not 
sufficient  to  ensure  that  it  is  the  stability  criterion  for  the  entire 
equation  system.  The  stability  condition  (Courant'-Friedrirhs-Lewy 
condition)  of  the  numerical  solution  when  considering  only  temporal  and 
linear  convective  terms  in  the  finite  difference  equations  of  motion  or 
convectivet-diffusion  is 


<Ax)mln 
At  <  (  3-t6  ) 


which  is  valid  for  the  explicit  discretization  of  a  hyperbolic  equati- 
on. Tn  eq.3.^6,  4t  is  the  time  interval  of  each  time  step,  (Ax),,,^ 
is  the  minimum  longitudinal  space  interval  used  in  the  grid  system,  and 
"max  's  tne  maximum  velocity  occurring  in  the  study  area. 

A  linear  parabolic  equation  including  the  temporal  and  diffusive 
terms,  tfhieh  are  time-  ana  space-centered  with  three  successive  time 


, 
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levels,  is  expressed  as 

n+1   n-1         n*»1     n<-l   n-1 


(  3.17  ) 


2At  Ax^ 


where  f  denote  velocity  or  sediment  concentration,  and  A  can  be 
momentum  or  mass  diffusion  coefficient.  The  stability  condition  for 
eq.3.17,  is  derived  using  complex  Fourier  series  analysis 
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The  derivation  of  the  inequality  condition  (eq.3.18)  is  given  in 
Appendix  A. 3 

The  stability  conditions  found  above  provide  guidance  for  select- 
ing Ax,  At,  and  momentum  or  mass  diffusion  coefficients. 
3.2.1  Calibration  and  Verification 

In  the  present  study,  calibration  is  defined  as  the  adjustment  of 
coefficients  based  on  an  assumption  that  the  relations  usea  in  the 
model,  in  which  the  coefficients  appear,  adequately  describe  the 
physical  phenomena  under  consideration.  The  unknown  coefficients  in 
the  model  are:  1)  longitudinal  and  vertical  momentum  diffusion  coeffi- 
cients, ex,  ez;  2)  longitudinal  and  vertical  mass  diffusion  coeffi- 
cients, £x,  tz;  and  3)  Manning  roughness  coefficient,  n.   Besides,  some 
physical  parameters  need  to  be  determined  empirically;  these  include 
the  settling  velocity,  ws,  and  the  critical  deposition  shear  stress, 
tC(j.  The  test  COEL-I  was  used  as  a  calibration  test,  such  that  its 
conditions,  i.e.  water  depth,  water  density,  sediment  specific  gravity, 
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and  concentration  profile  at  the  entidr""  t^f   the  closed'-end  channel, 
were  related  as  the  input  data  in  the  mathematical  model  simulation. 
In  addition,  for  grid  size  (which  is  based  on  how  detailed  predictions 
in  space  one  likes  to  obtain)  and  time  step  (which  is  referred  to  the 
stability  conditions  aforementioned),  values  of  Ax=0.5  m,  Az-1 .0 
cm,  and  4t=0.15  sec  were  adopted  for  simulating  all  tests  of  WES  and 
COEL.  Also,  ex,  ez  are  usually  assumed  to  be  identical  to  the  longitu*- 
dinal  and  vertical  momentum  diffusion  coefficients,  respectively,  by 
virtue  of  Reynolds  analogy.  Since  the  turbidity  currents  were  laminar 
in  all  the  laboratory  tests  in  the  present  study,  the  kinematic 
viscosity  of  water  was  used  for  ez  (i.e.  ez  ■  10*6  m2/see).  The 
longitudinal  momentum  diffusion  coefficient,  ex,  was  determined  with 
reference  to  the  stability  condition  given  as  eq.3.t8  and  was  found  to 

be  5.0  x  IO^j  m^/sec.  Manning's  roughness  coefficient,  n,  was  0.03 

1/3 
sec/m   .  This  n  value  selection  can  be  justified  by  comparing  with  an 

1  /  3 
illustrative  value  n  •  0.022  sec/m   ,  which  is  obtained  by  utilizing 

the  Moody  diagram  (Daily  &  Harleman,  1966)  and  eq.3.31  for  a  turbidity 

current  with  u  -  0.5  cm  and  an  interface  thickness  n  ■  5  cm.  The 

calibration  results  are  given  in  Chapter  5. 

Model  verification  was  carried  out  by  checking  how  well  the 

calibrated  model  reproduced  the  phenomena  in  other  tests.  Here,  C0EU-5 

was  selected  as  a  verification  test.   Comparison  between  the  results 

using  the  calibrated  numerical  model  and  experimental  measurements  for 

many  aspects  of  test  C0EL^5  is  given  in  Chapter  5. 
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3.3  Analytical  Developments 
3.3.1  Mean  Sediment  Concentration  Below  Interface 

In  order  to  examine  the  longitudinal  variation  of  mean  sediment 
concentration  below  interface  in  a  stationary  sediment  wedge,  an 
attempt  at  solving  the  steady  state  convective-diffusion  equation 
analytically  was  made  as  follows. 

The  physical  properties  at  each  layer  are  assumed  to  be  uniform. 
That  is,  the  upper  layer  is  assumed  to  be  sediments-free,  while  the 
suspension  concentration  in  the  lower  layer  is  assumed  to  be  uniform. 
Also,  the  velocity  profiles  in  both  layers  are  assumed  to  be  uniform. 
A  schematic  figure  is  shown  in  Fig. 3. 2.  In  the  lower  layer,  upward 
vertical  velocity  of  the  fluid  is  small  compared  with  the  settling 
velocity  of  suspended  sediment.  Mass  diffusion  in  both  x  and  z 
directions  is  assumed  to  be  negligible  as  compared  with  mass  convection 
by  flow.  Under  these  assumptions  the  steady^state  convective-diffusion 
equation  can  be  expressed  as 


3C      3C 

u w3  —  =  0  (  3.19  ) 

9x      9z 


Integrating  eq.3.^9  from  bottom  to  interface,  n,  then  using  the 
boundary  conditions  at  bottom  (C  ■  C^,  at  z  =  0)  and  interface  (C  ■  0, 
at  z  =  n),  an  equation  for  mean  concentration  gradient  is  obtained 


3Cb       w3 


Cb  (  3.50 


3x 


where  Cb  is  the  mean  suspension  concentration  in  the  lower  layer. 
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Figure  3.2.  Schematic  of  Theoretical  Two-Layered  Model. 
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Furthermore,  eq.3.50  is  Integrated  frog  *«0  O.e.  the  entrance  or 
the  close-enrt  channel)  to  x.  Since  both  w3  and  u  are  functions  of  ■/., 
the  integration  of  the  terra  on  the  right  hand  side  of  eq.3.50  with  dx 
cannot  be  carried  out  unless  explicit  forms  of  the  functions,  wa(x/ 
and  u(x),  are  known.  According  to  the  experimental  observations  from 
the  present  study,  both  sediment  settling  velocity  and  flow  velocity 
deer-eased  exponentially  with  distance  from  the  entrance  of  the  side 
channel.  Therefore,  ws  and  u  can  be  expressed  as  w3  ■  w3-| -expOa,;: 
and  u  ■  u, •exp(-82x),  respectively,  where  wsi  and  ui  are  the  sediment 
settling  velocity  and  the  flow  velocity  at  x  =  0,  and  Si ,  62  are 
constants.  The  thickness  of  the  lower  layer,  ij,  is  assumed  unchanged 
because  it  decreases  very  slowly  with  x.  Accordingly,  the  solution  for 
the  integration  of  eq.3.50  can  be  expressed  as 


Cb         ws1        (62  -  Bi)x 

in  —  -  -  r [  e         »  %   J   (  3.51  ) 

Cb1     uin(B2  "  61 ) 


where  Cbi  is  the  mean  sediment  concentration  at  x=0.  The  term 
exp[(82  uBi)x]  can  be  expanded  using  Taylor  series,  i.e.  1  + (82*61  )x+ 

[(B2"Bi )x]2/2+ ,  and  its  second  and  higher  order  terms  can  be 

omitted  if  the  exponent  (B2~Bi)x  is  much  less  than  unity  (i.e.  the 
local  relative  settling  velocity  w3/u  is  very  close  to  the  ratio  of 
wai  .'u<  at  the  entrance).  Also,  it  is  assumed  that  n  •  costant-H  and  u-| 
■  constant 'U^,  to  relate  these  parameters  (n  and  u-| )  to  the  total  water- 
depth  H  and  the  densimetric  velocity  u^,  respectively.  Finally,  the 
dimenslonless  form  of  the  relationship  for  the  mean  concentration  below 
interface  can  be  obtained  as  follows: 
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Cb             ws1    x 
—  -  exp  r  »   B( )( )  !  (  5.52  ) 

Cbl  uA    H 


where  t   is  a  constant  and  oust  De  determined  empirically.  The  investi- 
gation of  lock  exchange  flow  (O'Brien  &  Cherno,  1 93^;  *ih,  1965)  showed 
that  the  interface  at  the  gate  and  the  initial  front  speed  are  equal  to 
one  half  of  the  total  water  depth  and  densimetric  velocity,  respective"- 
ly.  Therefore,  the  B  value  in  eq.3.52  could  be  expected  to  be  larger 
than  1,  since  ui  is  smaller  than  the  initial  front  speed  and  the 
interface  at  any  location  in  the  closed^end  channel  is  lower  than  that 
at  the  entrance. 
3.3.2  Flow  Velocity 

The  unsteady  state  fluid  motion  in  the  closed-end  channel  is 
impossible  to  obtain  by  solving  analytically  the  full  set  of  equation 
of  motion  due  to  the  nonlinear  convection  term  involved  (eq.3.21).  As 
previously  mentioned,  there  are  three  distinct  phases  found  la  the 
front  motion  of  a  gravity  current,  and  the  third  phase  is  the  so»*called 
viscous  self-similar  phase,  In  which  the  front  motion  is  governed  by 
gravity  and  viscous  forces.  Huppert  (1982)  solved  the  flow  velocity 
within  the  gravity  current  under  this  phase  in  a  deep  water  channel. 
In  the  present  study,  a  gravityt-viscous  force  balance  is  also  consi- 
dered for  solving  the  flow  velocity  in  a  sediment  wedge  arrested 
in  a  shallow  water  channel. 

The  channel  may  be  assumed  two-dimensional,  and  the  Lengtfc  of  the 
wedge  will  be  considered  to  be  much  greater  than  the  still  water  depth 
H.  The  layer  above  the  interface  will  be  assumed  to  be  sediment-free: 
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all  the  sediment  being  confined  to  the  wedge.  Thus,  momentum  diffusion 
in  the  flow  is  more  important  in  the  vertical  z^direction  than  in  the 
horizontal  x<Mirection.   In  addition,  the  closed  end  cf  the  olosedfend 
channel  will  result  in  a  water  surface  slope  towards  the  entrance, 
consequently,  the  hydrostatic  head  induced  by  the  elevated  water 
surface  will  balance  the  excess  pressure  due  to  the  density  gradient  of 
sedimentWaden  water  in  the  channel  as  well  as  the  shear  stress 
generated  due  to  the  vertical  velocity  gradient. 

Under  these  conditions,  the  steady  state  equations  of  motion  in  x 
and  z  directions  can  be  expressed  as  follows: 

1  3p      32u 

0  ■  " *  ez  —  (  3.53  ) 

p  3x      3z2 


1  3p 

(  3.54  ) 

p  3z 


The  horizontal  pressure  gradient  term  in  eq.3.53  can  be  substituted 
by  the  result  obtained  by  taking  the  horizontal  derivative  of  the 
hydrostatic  pressure  from  eq.3.51.  In  the  upper  layer,  the  horizontal 
pressure  gradient  is  related  only  to  free  surface  slope,  but  in  the 
lower  layer,  this  gradient  is  related  to  both  the  free  surface  slope 
and  the  horizontal  density  gradient.   Therefore,  eq.3.53  can  be  split 
into  two  equations  corresponding  to  both  layers,  which  are  thus 


Y  u     3h 


e, . 


3z2     3x 


n  <  z  <  h         (  3.55  ) 
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32u     3h  pw    Pn  -  p,j    3n   g  ,n  3p 

ez  -  g ♦  ■ g  —  t  -  J  —  dz   0  <  8  <  n   (3.56  ) 

3z2     3x  p       p,,      3x   p  '  z   3x 

where  h(x)  is  the  elevation  of  free  surface,  and  pn  is  the  density  of 
sedimenf-water  mixture  at  the  interface.  The  second  term  on  the  right 
hand  side  in  eq.3.56  can  be  neglected,  since  the  fluid  density  at  the 
interface,  p_,  is  very  close  to  the  water  density,  pw.  Here,  the 
sediment  concentration,  C.  is  defined  as  the  sediment  mass  per  unit 
volume  of  water^sediment  mixture.  Therefore  the  density  of  the 
water-»sediment  mixture  can  be  calculated  from  p  =  pw+C  ( 1  *>1 /Gs ) ,  where 
Gs  is  the  specific  gravity  of  sediment.  The  derivation  of  the  rela^ 
tionship  between  p  and  C  is  given  in  Appendix  B.1.  Then  the  integra- 
tion of  3p/3x  in  eq.3.56  can  be  carried  out  by  using  the  expression  of 
eq.3,52.   In  order  to  solve  eqs.3.55  and  3.56,  the  following  conditions 
are  utilized;  1)  no  shear  stress  at  the  water  surface  (i.e.  3u/3z  =  0, 
at  z  ■  h),  2)  nc-slip  boundary  condition  at  the  bed  (i.e.  u  =  0,  at 
z  =  0),  3)  the  horizontal  velocity  is  continuous  at  the  interface,  M) 
the  slope  of  vertical  velocity  profile  is  continuous  at  the  interface, 
and  5)  the  volumetric  rate  of  inflow  in  the  lower  layer  Is  equal  to  the 
rate  of  outflow  in  the  upper  layer  at  every  location  in  the  closed^end 
channel  (i.e.  J  udz  --  0  ).  Details  of  solving  equations  3.55  and  3.56 
are  given  in  Appendix  B.2. 

The  final  solution  of  flow  velocity  in  upper  and  lower  layers  are 
obtained  as  follows: 


gs0H2  1 

[  0.5  c2  "  C  '  ]    e  <  t.    <   1     (  3.57  ) 

ez  3(1  -  0.255) 
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[ ♦  (0.5 )  C^ 

ez     3530  ■  0.255*  52(1  I*   0.255) 

1 
+  ("1  * )  c  ]        0  <  c  <  F,        (  3.58  ) 

4.(1  *  0.255) 


where  the  water  surface  slope  s0  =  3h/3x,  the  dimensionless  vertical 
coordinate  t,  -   z/h,  and  the  dimensionless  interface  elevation  5  «  n/h. 
Eqs.3.57  and  3.58,  in  which  s0  and  5  are  the  function  of  x,  are  valid 
everywhere  in  the  closed^end  channel  except  at  the  entrance  where  the 
water  body  is  connected  to  open  water,  and  the  flow  condition  there 
is  critical  as  found  by  Schijf  and  Schonfeld  (1953).  Furthermore, 
normalizing  eqs.3.57  and  3.58  by  the  local  surface  outflow  velocity, 
us,  these  equations  become 

u    (1.5  -  0.3755)  S2  "  (3  *  0.755)  C  +  1 

5  <  C  <  I    I  3.59  ) 

us  "0.5  +  0.3755 

u    U/5)3+("3+1.552--0.37553)U/5)2+(3'-35+0.7552)U/O 

= 0  <  ;  <  5 

us  -0.5  +0.3755 

(  3.60  ) 


Note  that  us  is  always  negative  because  the  flow  is  towards  the 

entrance. 

3.3.3  Sediment  Flux 

Since  analytical  velocity  solutions  of  eqs.3.57  and  3.58  are  not 
applicable  to  the  entrance  of  the  closedf"end  channel,  an  attempt  to 
develop  a  practical  method  for  determining  the  sediment  flux  at  the 
entrance  of  the  channel  is  made  in  this  section. 
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According  to  present  and  previous  investigations,  the  Initial 
front  ap^ed,  Uf| ,  i.e.,  the  particle  velocity  at  the  entrance  immediat*1 
ely  after  the  gate  is  removed,  is  approximately  equal  to  one  half  of 
the  de.'isimetric  velocity  u^.  When  the  front  of  the  turbidity  current 
proceeds  downstream  along  the  closed-end  channel,  the  flow  velocity  at 
the  entrance,  u-|  ,  decreases  and  rapidly  reaches  a  steady  state.  The 
particle  velocity  variation  with  the  elapsed  time  at  several  locations 
in  the  closed-end  channel  is  shown  in  Fig. 5. 5  in  Chapter  b,   which 
presents  the  simulation  results  obtained  by  the  numerical  model.  From 
the  results,  it  is  observed  that  the  period  between  the  maximum 
velocity,  i.e.  during  the  front  passing,  and  the  steady  state  velocity 
is  much  shorter  compared  with  the  overall  duration  of  the  experiment, 
and  the  velocity  at  steady  state  is  proportional  to  the  maximum 
velocity.  Therefore,  it  is  reasonable  to  assume  that  a  constant 
particle  velocity  at  the  entrance  of  the  closed-end  channel  is  esta* 
blished  since  the  removal  of  the  gate  and  u-|  ■  constant-Ufi  .  There- 
fore, u-|  can  be  expressed  as  u-|  =  a  u&,  since  Ufi  =  0.5  u^,  where  a  is 
a  constant  and  the  densimetric  velocity  u^  »  [ApigH/pw]1'^.  The 
sediment  influx  through  the  lower  layer  at  the  entrance  per  unit  area 

is 

gH  1     1/2  _3/2 

s  .  urc,  -  a  I   — (1  -  — )  ]    C,         (  3.61  ) 

Pw  0* 

where  Cj  is  the  depth-averaged  concentration  at  the  entrance  of  the 

closed"end  channel.  Note  that  the  sediment  flux  into  the  closed-end 

_3/2 
channel  is  proportional  to  C^ 

Prediction  of  sediment  flux  computed  by  eq.3.61  is  compared  with 

laboratory  experimental  results  in  Chapter  5. 
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3.3.1  Water  Surface  Rise 

While  a  stationary  sediment  wedge  is  established  In  the  closed-end 
channel,  the  water  surface  in  the  region  beyond  the  toe  of  the  sediment 
wedge  must  be  raised  in  order  to  balance  the  excess  force  induced  by 
density  difference  at  the  entrance.  In  this  section,  an  attempt  is 
made  at  solving  for  the  maximum  water  surface  rise,  AHraax,  based  on 
momentmum  conservation.  A  schematic  diagram  appears  in  Fig. 3. 3. 

According  to  the  equation  of  momentum  conservation,  the  following 
expression  can  be  deduced 

F1  *  P1LQ2"2  ■  F2  +  P1uQl"1  <  3.62  ) 

in  which: 

F-|  ,  F2  ■  hydrostatic  forces  acting  on  the  cross  sections  at 
the  entrance  and  wedge  toe,  respectively 

P1U»  P1L  =  depth-averaged  densities  in  upper  and  lower  layers 
at  the  entrance,  respectively 

Ql  ,  Q2  -  discharges  in  upper  and  lower  layers  at  the  entrance, 
respectively  and 

u-| ,  U2  ■  depth-averaged  velocities  in  upper  and  lower  layers 
at  the  entrance,  respectively. 

However,  Q-j  ■  Q2  must  be  satisfied  due  to  the  volumetric  conserva*- 

tion  of  water  at  steady  state.  Accordingly,  u-|  is  approximately  equal 

to  U2,  since  the  interface  elevation  at  the  entrance  is  about  one  half 

of  the  water  depth.  Therefore,  eq.3.62  reduces  to 

F2  -  F1  "  (PIL  "   7luOQ2"2  (  3.63  ) 

in  which: 

H   H 
Pi  »  *  1  (  1  p,(z)gdz)  dz 
0   0 
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Figure  3.?.  Schematic  of  Stationary  Sediment  Wedge 
and  Water  Surface  Rise. 


,, 


B  2 

r2  -  -  (H  «  tHmax)   pwg,  arid 
2 

B  =  channel  width. 

In  eq.3.63,  the  term  on  the  right  hand  side  is  assumed  to  Be 
negligible  because  (p-\i  -■   puj)  is  small.   Also,  for  F2,  the  term  with 
the  order  of  (4Hmax^'  ls  neglected,  since  tHmax   is  very  small. 
Consequently,  the  rise  of  water  surface  at  the  region  beyond  the  toe  of 
the  stationary  sediment  wedge,  AHmax,  is  given  as 

1    1    H    H  H2 

AHmax  "  -  I  —  J   M  Pi(z)dz  ]  dz  «  —  }       (  3.6U  ) 
H   pw  0    z  2 

where  p^tz)  is  the  density  profile  at  the  entrance  of  the  closed^end 

channel. 

3-3.5  Tide-Induced  to  Turbidity  Currenflnduced  Deposition  Ratio 

In  many  coastal  regions,  sedimentation  occurring  in  closed^end 
channels  have  resulted  from  the  combined  effect  of  turbidity  currents 
and  tidal  motions,   it  is  of  concern  from  an  engineering  viewpoint  to 
distinguish  the  contribution  of  tide'-induced  deposition  from  total 
deposition.  In  this  section,  an  attempt  to  estabilish  a  relationship 
between  tide-induced  deposition  ana  turbidity  current-induced  deposi- 
tion is  m.^de.  A  relationship  for  sediment  flux  into  the  closed*-end 
channel  due  to  turoidxty  currents  has  been  derived  as  eq.3.6l.  The 
mass  of  sediment  deposition  over  a  tidal  period  can  be  computed  aj 


«0'SB  (H/2)  T  (  3.65 

■ 


where  Mn  Is  the  mass  of  deposition  Indiwsd  by  turbidity  currents,  and 
is  the  tidal  period. 

Regarding  tide-Induced  deposition,  a  simple  way  to  compute  r.he 
deposited  mass  over  a  tidal  cycle  's  using  the  concept  of  tidal  pri3m. 
This  assumes  that  suspended  sediment  moves  Into  the  channel  during 
flood  tide,  then  parts  of  it  settles  at  slack  water,  during  ebb  tide 
the  remaining  suspended  sediment  flows  out  of  the  channel.  According* 
ly,  the  deposition  mass  induced  by  tidal  motion  over  a  tidal  cycle,  Mj, 
is  obtained  as 

MT  =  20,  a0  L,  B  C,  (  3.66  ) 

where  a0  is  the  tidal  amplitude  and  0  is  the  portion  of  sediment 
depositing  Inside  the  channel  at  slack  water.  Q  ■  1  if  no  sediment 
outflow  during  ebb  tide,  and  a  <  1  if  some  sediment  flows  out  of  the 
channel.   Consequently,  the  ratio  of  tide-induced  to  turbidity  current 
■-induced  deposited  mass  (for  a  «  1  case)  can  be  obtained,  by  dividing 
eq.3.66  by  eq.3.65,  as 

My         Api  -1/2  a0   L 

_  ,  11. n2  ( )    (_)  (_)  (  3.6T  ) 

MD         Pw       H    L' 

where  Api  is  the  excess  density  at  the  channel  entrance,  and  L'  is  the 
tidal  wave  length  in  shallow  water-,  L'«(gH)1//2T.   In  eq.3.65,  three 
basic  dimensionless  parameters  are  included,  i.e.  Api/pw,  the  relative 
excess  density  at  the  channel  entrance,  a0/H,  the  relative  tidal 
amplitude,  and  L/L1 ,  the  channel  length  to  tidal  length  ratio. 


CHAPTER  1 
EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


Two  sets  of  experiments  were  conducted  in  the  present  study.  The 
first  was  performed  in  a  flume  system  at  the  U.  S.  Army  Corps  of 
Engineers  Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi. 
These  experiments  served  as  the  preliminary  study  for  an  understanding 
of  turbidity  current  and  sedimentation  in  a  closed-end  channel.  The 
experimental  results  have  been  presented  by  Lin  and  Mehta  (1986). 
During  the  tests,  considerable  temperature  gradients  in  the  flume  were 
observed.  The  second  set  of  experiments  was  conducted  in  a  smaller 
size  of  flume  system  in  a  room  with  temperature  control  at  the  Coastal 
Engineering  Laboratory  (COEL)  of  the  University  of  Florida.  Observa- 
tions in  detail  on  the  behavior  of  turbidity  current  and  associated 
sedimentation  were  made. 

1.1  Experiments  at  WES 
1.1.1  Experimental  Sef-up 

The  laboratory  tests  at  WES  were  performed  In  a  T-shaped  flume 
(Fig. 1.1).  A  9.1  m   long,  0.16  m  deep,  and  0.23  ■  wide  closed'-end 
channel  was  orthogonally  connected  to  a  100  m  long,  0.16  m  deep, 
and  0.23  m  wide  main  channel  at  the  approximate  midpoint  of  the  main 
channel.  Fig. 1.1  shows  the  experimental  set-up  schematically.  Both 
main  and  closed^end  channels  were  horizontal  and  made  of  6  mm  thick 
lucite.  The  walls  of  the  main  channel  were  artificially  roughened 
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using  the  lucite  elements.   Previous  measurements  In  the  main  channel 
showed  that  the  main  channel  had  a  Manning's  roughness  coefficient  n  ■ 
0.022  with  the  artificially  roughened  walls  (Dixit  et  al.  1982), 
Other  components  of  the  experimental  sef-up  comprised: 

1)  A  6.1  m  square  and  1.5  m  deep  basin  at  each  end  of  the  main 
channel.  The  upstream  basin  including  an  overflow  weir  structure 
served  as  the  headbay  and  the  basin  at  the  downstream  end  as  the 
tailbay. 

2)  Fluid  injection  apparatus  consisting  of  a  water  tank  of  111  liter 
capacity  and  a  slurry  tank  of  1500  liter  capacity.  As  shown  in 
Fig. 1.1,  a  quick  action  valve  controlled  the  water  supply  between 
the  water  tank  and  the  slurry  tank.  Four  sets  of  15  cm  propellers 
driven  by  a  915  rpm  electric  motor  and  a  mixing  pump  were  provided 
in  the  slurry  tank.  One  end  of  a  2.5  cm  diameter  hose  was  connect- 
ed with  the  outlet  of  a  rotameter,  and  the  other  was  placed  in 

the  main  channel  6.0  m  upstream  from  the  entrance  center  of  the 
Closed-end  channel. 

3)  Venturi  meter  attached  to  the  water  supply  pipe  towards  the  headbay 
In  order  to  monitor  the  flow  rate  from  the  supply  sump  to  the 
headbay. 

1)  V-notch  weir  to  measure  discharge  from  the  main  channel. 
1.1.2  Auxiliary  Equipment 

1 )  Point  gauges  for-  the  measurement  of  water  surface  and  bed  eleva- 
tions at  the  main  and  closed">end  channels. 

2)  Digital  electric  thermometers,  shown  in  Fig. 1.2,  to  monitor  the 
temperature  variation  of  fluid  in  the  main  and  closed«end  chan- 
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nels.  The  thermometers  oould  read  up  to  the  second  decima.1  place 
in  degrees  Fahrenheit. 
3)  A  rotating  cup  current  meter  for  measuring  the  distribution  of   the 
horizontal  velocity  over  the  vertical  in  the  main  channel.  The 
frequency  of  wheel  revolution  could  be  converted  to  the  flow 
velocity  from  a  calibration  curve. 
1)  The  suspension  sampling  apparatus  consisted  of  several  sampling 
devices  and  a  vacuum  pump.  Each  sampling  device  included  four  taps 
aligned  vertically,  as  shown  in  Fig.!). 3;  the  inner  diameter  of  each 
tap  was  3  mm.  One  end  of  each  tap  was  connected  by  a  long, 
flexible  plastic  tube  to  a  set  of  50  ml  plastic  jars  located  in  a 
closed  rectangular  plastic  box.  In  all,  30  Jars  were  placed  in  the 
closed  box  for  collecting  the  suspension  samples  from  different 
locations  and  elevations  during  each  test.  The  outlet  of  the 
closed  box  was  connected  to  a  vacuum  pump. 
1.1.3  Test  Materials 

Five  kinds  of  sediment  were  used  in  the  laboratory  experiments  at 
WES:  1)  kaollnite,  2)  silica  flour,  3)  flyash  I,  15  flyash  II,  and  5) 
Vicksburg  loess.  The  kaolinite  was  available  commercially  (Feldspar 
Coorporation,  Edgar,  Florida).  The  size  distribution  of  the  kaolinite 
is  given  in  Fig. 4.1.  The  median  particle  diameter  was  1.2  ym.  The 
cation  exchange  capacity  was  found  to  be  approximately  6.0  milliequiva- 
lents  per  100  grams  (Parchure,  1983).  Its  specific  gravity  was  2.57 
(Lott,  1986). 

The  particle  size  distribution  of  silica  flour  obtained  by  a 
Sedigraph  Particle  Size  Analyzer,  which  u3es  an  x-ray  beam  to  detect 
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Figure  1.2.  A  Digital  Electric  Thermometer  at  WES. 


Figure  1.3.  Suspension  Sampling  Apparatus  Deployed  in  the  WES  Flume. 
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7? 
the  concentration  decrease  due  to  settling  in  a  sample  cell,  is  shown 
in  Fig. 1.1.  The  median  particle  diameter  was  7  um.  The  specific 
gravity  of  the  sediment  was  2.65.  No  organic  matter  was  found. 

The  particle  size  distribution  of  flyash  I  was  obtained  using  the 
Sedigraph  Analyzer  and  is  shown  in  Fig. 1.1.  The  median  particle 
diameter  was  11  urn.  The  specific  gravity  of  the  sediment  was  2.15. 
Loss  on  ignition  was  found  to  be  3.7  percent. 

The  particle  size  distribution  of  flyash  II  was  also  obtained  by 
the  Sedigraph  Analyzer  and  is  shown  in  Fig. 1.1.  The  median  diameter 
was  10  vim.  The  specific  gravity  of  the  sediment  was  2.37.  Loss  on 
ignition  was  found  to  be  1.0  percent. 

The  particle  size  distribution  of  the  Vicksburg  loess  was  obtained 
by  Hydrometer  analysis  and  is  shown  in  Fig. 1.1.  The  median  particle 
diameter  of  the  sediment  was  18  um.  X-ray  diffraction  analysis  was 
conducted  at  WES  in  order  to  identify  the  clay  and  nonclay  minerals  of 
the  loess.   It  was  found  to  consist  of  five  clay  minerals,  which  are 
chlorite,  illite,  kaolinite,  montmorillonite,  and  vermiculite,  as  well 
as  three  major  nonclay  minerals,  which  are  quartz,  feldspar,  and 
dolomite.  The  cation  exchange  capacity  was  found  to  be  approximately 
19  milliequivalents  per  100  grams.  The  specific  gravity  was  2.68,  and 
the  organic  content  was  found  to  be  0.6  percent. 

As  mentioned  in  Chapter  3,  the  settling  velocities  of  raw  fine 
sediments  used  in  both  WES  and  C0EL  flume  tests  were  determined  by 
settling  column  analysis.  Two  different  sizes  of  settling  column  were 
used,  as  shown  in  Fig. 1.5;  one  was  1.15  m  high,  10  cm  square,  (design 
nated  as  column  no.1),  and  the  other  was  1.9  m  high,  10  cm  in  diameter, 
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Figure  U.S.  Two  Different  Size  Settling  Columns  with  a  Mixing  Pump. 
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(designated  as  column  no. 2).   A  mixing  pump  was  used  to  help  fully  mix 

the  water  and  fine  sediment  by  passing  air  bubbles  into  the  column 

prior  to  the  beginning  of  a  test.  Inasmuch  as  the  settling  velocity 

is  dependent  upon  sediment  concentration  for  cohesive  sediment,  several 

tests  were  conducted  with  different  initial  concentrations  for  each 

sediment.  The  relationship  of  median  settling  velocity  to  the  initial 

sediment  concentration  for  each  sediment  is  shown  in  Fig. 1.6.  A 

detailed  description  of  the  column  test  procedure  is  given  by  Lott 

(1986). 

The  chemical  composition  of  water  used  in  the  tests  at  the  WES 

flume  is  given  in  Table  1.1  (Dixit  et  al.  1982). 

Table  1.1  Chemical  Composition  of  a  water  sample  from 
the  100  m   flume  (after  Dixit  et  al.  1982) 


Ions 


PPm 


Ca  1.9 

Mg++  0.7 

Na+  13.1 

K+  6.6 

Mn  0.0 

+  +  + 

Fe  0.0 

Cl"  18.3 

Electrical 
Conductivity  0.19  mmho/cm 

PH  7.8 
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Figure  1.6.  Median  Settling  Velocity  vs.  Concentration. 
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1)  The  flow  rate  of  fresh  water  to  the  headbay  was  adjusted  until  a 
calibrated  reading  of  the  venturi  meter  was  attained.  Thereafter, 
the  flow  was  kept  running  for  about  1.5  hours  in  order  to  establish 
a  steady  flow  with  a  desired  water  depth  and  velocity  in  the  main 
channel.  Within  the  period  of  establishing  the  steady  flow,  the 
gatf  located  at  the  entrance  was  kept  open. 

2)  The  appropriate  amounts  of  test  sediment  and  water  were  mixed  in 
the  slurry  tank  to  obtain  a  desired  slurry  concentration.  The 
mixing  propellers  and  the  pump  started  operating  before  the 
sediment  was  added  to  the  tank.  For  easy  observation  of  the 
movement  of  turbidity  current  in  the  closed-end  channel,  Rhodamine 
dye  (red)  was  added  to  the  slurry  tank. 

3)  Fresh  water  from  the  water  tank  was  pumped  into  the  main  channel 
through  the  injection  hose  at  a  flow  rate  equal  to  the  slurry 
injection  rate  required.  The  injection  flow  rate  was  selected  to 
obtain  the  desired  suspension  concentration  in  the  flow  of  the  main 
channel. 

t)  Until  the  flow  in  the  main  channel  became  stable,  the  gate  of  the 
closed^end  channel  was  closed,  a  vertical  velocity  profile  at  the 
main  channel,  1.5  m  dowmstream  from  the  entrance  center  of  the 
closed^end  channel,  was  obtained. 

5)  A  quick  action  valve  was  switched  to  change  the  Injection  from 
fresh  water  to  slurry,  then  the  sediment  slurry  injection  to  the 
main  channel  was  kept  at  the  constant  flow  rate  throughout  the 
test. 
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6)  With  a  constant  flow  rate  of  slurry  la  the  main  channel,  the  gate 
of  the  closed'-end  channel  was  removed  allowing  the  sediment-laden 
water  to  Intrude  Into  the  closedftend  channe]  along  the  lower 
layer.  The  Interface  between  the  fresh  water  and  the  sediment 
front  was  traced  several  times  during  a  test  on  transparencies, 
which  were  attached  on  the  side  wall  with  a  scale  reference. 

7)  Thirteen  point  gauges  deployed  in  the  main  and  closed**end  channels 
are  shown  in  Fig. 4. 7.  The  distance  scale  of  the  point  gauge 
location  is  the  same  as  the  length  scale  of  the  closed^end  cha- 
nnel. Water  surface  and  bed  elevations  were  measured  using  point 
gauges  in  the  main  and  closedr-end  channels  at  three  different  times 
during  a  test,  (i.e.  during  a  constant  flow  rate  of  fresh  water  in 
the  main  channel  with  the  entrance  gate  closed,  a  constant  flow 
rate  of  fresh  water  in  the  main  channel  with  the  gate  open,  and  a 
constant  flow  rate  of  sedimentKladen  water  in  the  main  channel  with 
the  gate  open) . 

8)  Three  probes  of  digitalized  electric  thermometers  were  mounted 
on  point  gauge  supports  in  the  main  and  closed^end  channels  as 
shown  in  Fig. 1. 7.  Vertical  temperature  profiles  were  taken 
periodically  at  all  three  locations. 

9)  A  dye  injection  technique  and  floating  tracer  particles  were  used 
to  measure  vertical  and  lateral  velocity  profiles  in  the  closed'-end 
channel. 

10)  Seven  suspension  sampling  devices  were  deployed  in  the  main  and 

closed'-end  channels  as  shown  in  Fig.  1.7.  The  distance  scale  of  the 
sampling  device  location  was  the  same  as  the  length  scale  of  the 
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Figure  1.7.  Locations  of  Measurements  In  WES  Tests  (Plan  View). 
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closed>-end  channel.  Suspension  sampler  of  the  water*sediment 
mixture  were  withdrawn  from  all  taps  using  a  vacuum  pump,  towards 
the  end  of  a  teat. 

11)  After  the  suspension  samples  were  withdrawn,  the  gate  of  the 
closed*end  channel  was  closed,  but  a  slow  release  of  water  out  of 
the  channel  was  allowed  without  disturbing  the  deposited  sediment 
in  the  closed»-end  channel. 

12)  In  four  out  of  thirteen  tests,  the  deposits  were  collected  at 
various  subsections  from  the  bed  of  the  closed«end  channel  at 
approximate  21  hours  after  the  test  was  completed. 

A  summary  of  test  conditions  for  the  thirteen  tests  at  WES  is 
given  in  Table  1.2. 

1.2  Experiments  at  COEL 
1.2.1  Experimental  Set-up 

A  total  of  fourteen  tests  was  performed  at  COEL  in  a  r-shaped 
flume,  which  consisted  of  a  olosed>-end  channel  orthogonally  connected 
to  a  main  channel  at  2.83  m  from  the  upstream  end  of  the  main  channel. 
The  closedhend  channel  was  11.7  m  long,  0.1  m  wide,  0.19  m  deep,  and 
was  constructed  of  plexiglass  of  1.3  cm  thickness.  The  main  channel 
was  1.26  m  long,  0.16  a  wide,  0.19  m  deep,  and  was  also  made  of 
plexiglass.  A  plexiglass  gate  was  located  at  the  junction  of  the 
closed»end  channel  with  the  main  channel.  In  a  total  of  the  fourteen 
tests,  twelve  tests  (Nos.  1-8  &  11-12  &  13M1)  may  be  considered  steady 
flow  and  two  (Nos.  9  8.  10)  were  te3ts  with  tidal  effect.  The  hydraulic 
system  components  are   shown  in  Fig. 1.8. 
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Figure  14.8,   Schematic  of  Hydraulic  System  at  COEL   (Plan  View). 
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'['tie  components  of  the  experimental  set"up  for  the  steady  flow 
tests  and  the  tests  with  tidal  effect  comprised  the  following: 

1)  ft  1 .8  m  square  and  0.3  m   in  deep  basin  used  as  the  headtay.  An 
overflow  weir  structure  was  located  on  the  upstream  side  wall  of 
the  headbay  as  shown  in  Fig. it. 8.  The  elevation  of  the  gate  oould 
be  fixed,  or  varied  continuously  to  simulate  a  sinusoidal  tide 
using  a  variables-speed  electric  motor. 

2)  A  fresh  water  diffuser  lying  on  the  bottom  of  the  headbay.  The 
diffuser  was  made  of  a  section  of  perforated  PVC  pipe.  Fresh 
water  entered  the  headbay  through  either  side  of  the  diffuser  so 
that  water-  surface  disturbances  could  be  reduced  to  the  minimum. 

3)  A  rectangular  basin,  1.83  by  0.89  m,  used  as  the  tailbay.  A 
plexiglass  tail  gate  constructed  at  the  junction  of  the  main 
channel  with  the  tailbay  could  be  adjusted  to  obtain  the  desired 
water  depth  and  flow  velocity  in  the  main  channel. 

()  A  mixing  tank  consisting  of  a  primary  tank  and  a  secondary  tank. 
Each  tank  was  1.21  m  square,  1.82  m  in  deep  and  had  an  individually 
controlled  air  bubbler  apparatus.  The  air  bubbler  consisted  of 
steel  piping  perforated  with  holes  so  that  rising  air  bubbles 
induced  turbulence  and  mixing. 

5)  A  slurry  injection  diffuser  located  at  the  head  of  the  main 

channel.  A  slurry  pump  with  variable  speed  control  was  connected 
with  the  outlet  of  the  primary  tank  and  with  the  slurry  injection 
diffuser.  The  slurry  pump  was  used  to  establish  a  desired  flow 
rate  from  the  primary  tank  to  the  main  channel  through  the  slurry 
injection  diffuser. 
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6)  A  float  switch/level  control  pump  system  consisting  of  a  float 
switch  in  the  primary  tank  and  a  level  control  pump.  In  order  to 
obtain  a  constant  flow  rat:;  of  slurry  out  of  the  primary  tank,  a 
relatively  constant  water  surface  elevation  in  the  primary  tank 
had  to  be  maintained.  When  the  water  surface  elevation  dropped 
below  a  preset  level,  the  level  control  pump  automatically  turned 
on  and  began  pumping  the  slurry  from  the  secondary  tank  into  the 
primary  tank. 
Detailed  descriptions  of  the  use  of  the  control  valves  in  the 
piping  system,  flow  direction  in  the  pipeline,  and  the  design  of  the 
slurry  injection  diffuser  are  given  by  Lott  (1986). 
4.2.2  Auxiliary  Equipment 

1)  The  suspension  sampling  apparatus  was  deployed  at  eight  locations 
as  shown  in  Fig. 1. 9.  At  each  location,  seven  taps,  each  of  3  mm 
diameter,  were  vertically  and  equally  spaced  along  the  side  wall 
over  a  depth  of  10.5  cm.  Each  tap  was  connected  to  a  short  length 
of  a  plastic  tube.  The  plastic  tubes  were  upheld  by  a  tube  holder 
when  not  used. 

2)  Two  in-line  flowmeters  were  installed  in  the  fresh  water  pipeline 
and  slurry  Injection  pipeline,  respectively.   A  device  connected  to 
each  flowmeter  converted  the  rotation  frequency  of  the  flowmeter 
impeller  into  an  electrical  signal.  The  output  voltage  read  with  a 
volt  meter  corresponded  to  the  flow  rate  in  the  pipe. 

3)  An  electric  temperature  probe/meter  was  used  to  obtain  the  water 
temperature  at  any  desired  depth  and  location. 

1)  Two  point  gauges  were  used  for  measuring  the  water  surface  and  bed 
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Figure  ij.9.  Locations  of  Measurements  and  Sampling  Apparatus. 
Measurement  Locations  (a)  Plan. 
Suspension  Sampling  Apparatus  (b)  Plan,  (c)  Elevation. 


elevations  at  any  desired  location.  The  point  gauge  oould  read 
with  a  precision  of  0.1  mm. 

5)  A  13  mm  diameter  electromagnetic  current  meter,  mounted  on  a   point 
gauge  support  shown  in  Fig. 4. 10,  was  used  to  measure  vertical 
velocity  profiles  in  the  main  channel. 

6)  Sony  video  equipment  was  used  to  measure  front  speed,  observe  front 
shape,  and  measure  vertical  velocity  profiles  in  the  closed-end 
channel.  In  order  to  have  a  scale  reference  in  the  picture 
whenever  the  video  camera  viewed  through  the  closed^end  channel,  a 
plastic  tape  indicating  the  longitudinal  distance,  transparencies 
with  a  grid  reference,  and  transparent  metric  rulers  were  deployed 
along  the  channel  and  on  the  side  wall  of  the  channel.  The  video 
equipment  mainly  consisted  of  a  video  camera,  a  tape  recorder,  a 
time/date  generator,  and  a  TV  monitor.  A  dye  injection  technique 
in  conjunction  with  the  video  equipment  was  used  to  obtain  the 
vertical  velocity  profiles  at  the  desired  locations  in  the  closed* 
end  channel. 

7)  A  vacuum  pump  was  used  to  collect  the  deposited  sediment  from 
various  bed  subsections  in  the  main  and  closed*end  channels. 
1.2.3  Test  Materials 

Three  kinds  of  sediments  were  used  in  the  tests  at  COEL:  1) 
kaollnite,  2)  flyash  III,  and  3)  Cedar  Key  mud.   The  kaolinite  was  the 
same  as  that  used  at  WES. 

The  particle  size  distribution  of  flyash  III  was  obtained  using 
pipette  analysis  and  is  shown  in  Fig. 1.1.  The  median  particle  diameter 
was  13.7  (jm.  The  specific  gravity  of  the  sediment  was  2.63. 


Figure  4.10,  Electromagnetic  Current  Meter. 
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The  particle  size  distribution  of  Cedar  Key  mud  was  obtained  cud  * 
pipette  analysis  and  is  shown  in  Fig.".').  The  dark  soft  mud  was 
collected  from  Cedar  Key  on  the  Gulf  Cc?-st  of  Florida,  The  median 
particle  diameter  was  2.Z   um,  and  the  specific  gravity  of  the  mud  was 
2.36. 

A  settling  column  analysis  was  alsc  conducted  for  flyash  III  and 
Cedar  Key  mud,  and  the  relationships  of  median  settling  velocity  to  the 
initial  sediment  concentration  are  shown  in  Fig. 1.6. 

Tap  water  was  used  as  the  fluid  in  these  experiments,  the  chemical 
composition  of  this  water  is  shown  in  Table  4.3  (Dixit,  1982). 
1.2.1  Experimental  Procedure 

The  procedure  for  the  experiments  at  COEL  was  somewhat  different 
from  that  at  WES  due  to  the  different  laboratory  set*up,  equipment  and 
measurement  needs.  The  procedure  at  COEL  was  as  follows. 

Table  1.3 
Chemical  Composition  of  the  Tap  Water  (after  Dixit,  1982) 


CI  26  ppm 

NO3  0.07  ppm 

Fe  0.5  ppm 

K  1 . I  ppm 

Ca  25  ppm 

Kg  16  ppm 

Na  10  ppm 

Total  Salts  278  ppm 

PH  8.5  ppm 
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Steady  Flow  tests.  The  experimental  procedure  described  here  is 
for  the  steady  flow  tests. 

1)  At  the  beginning,  a  steady  flow  with  the  desired  flow  rate  and 
suspension  concentration  in  the  main  channel  was  established.  The 
operation  involved  adjustments  of  the  tail  gate  at  the  main  channel 
and  overflow  weir  at  headbay,  as  well  as  the  controls  of  the  water 
surface  elevation  in  the  channels  and  slurry  concentration  in  the 
primary  mixing  tank  in  coordination  with  the  in-line  flowmeters  and 
voltmeter.  A  detailed  description  of  this  operation  is  given  by 
Lott  (1986). 

2)  After  the  steady  flow  was  established  in  the  main  channel,  the  gate 
of  the  closed*end  channel  was  closed.  Subsequently,  the  main 
channel  velocity  profile  was  measured  using  the  electromagnetic 
current  meter.  The  first  set  of  suspension  samples  was  withdrawn 
at  tap  locations  A  and  B.  The  temperature  of  the  flow  in  the  main 
channel  and  of  the  water  in  the  closed^end  channel  (with  the  gate 
closed)  were  recorded.  All  measurement  locations  can  be  seen  In 
Fig. 1.9. 

3)  Then  as  the  gate  at  the  entrance  was  removed,  the  slurry  in  the 
main  channel  was  allowed  to  move  into  the  closed-end  channel.  As 
the  front  of  turbidity  current  moved  down  the  channel,  the  video 
camera  was  moved  alongside  to  continously  record  the  travelling 
time,  position,  and  shape  of  the  front  together  with  the  nearby 
scale  references. 

1)  The  suspension  samples  were  withdrawn  from  the  taps  immediateiy  as 
the  front  passed  each  sampling  location,  so  that  the  concentration 


92 
of  the  front  could  be  determined. 

5)  At  several  times  during  the  teat,  the  Interface  between  the  clear 
water  and  the  turbid  water  was  tracea  on  the  transparencies, 
which  were  attached  on  the  side  wall  of  the  closed'-end  channel. 

6)  During  the  test,  two  types  of  velocity  measurements  in  the  closed" 
end  channel  were  made.  One  was  surface  velocity,  which  was 
obtained  by  measuring  the  required  time  for  a  dye  parcel  (methyl 
blue)  to  move  over  a  known  distance  of  5  or  10  cm.  At  each 
location,  the  measurement  was  repeated  several  times  in  order  to 
obtain  a  mean  value.  The  other  measurement  was  the  vertical 
velocity  profile,  which  was  obtained  by  Injecting  a  vertical  streak 
of  dye  and  recording  its  horizontal  displacement  either  by  the 
video  camera  or  by  tracing  it  on  the  transparency. 

7)  A  set  of  suspension  samples  was  withdrawn  from  all  taps  In  the 
layer  of  the  watert-sediment  mixture  in  the  closed^end  channel  and 
the  main  channel  at  locations  A  and  B,  when  the  secondary  mixing 
tank  was  almost  empty. 

8)  Vertical  temperature  profiles  were  recorded  periodically  in  the 
main  and  closed>-end  channels  in  order  to  monitor  the  temperature 
variation  with  time  and  space  throughout  the  test. 

9)  If  a  test,  was  of  a  comparatively  long  duration,  the  nearly  empty 
secondary  tank  had  to  be  refilled  and  well  mixed  with  the  appro- 
priate amounts  of  tap  water  and  the  sediment.  During  the  refill 
ling,  a  measurable  drop  in  the  water  level  in  the  primary  tank  was 
observed. 

10)  Just  before  finishing  a  test,  a  set  of  suspension  samples  was 
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collected  from  all  taps  In  the  $l«sed*ttnd  channel  and  the  main 
channel  at  locations  ?.   «nd  B. 

1 1 )  At  tne  end  of  a  test,  the  gate  of  the  ciosed*-end  channel  and  the 
tall  gate  in  the  main  channel  were  closed.  The  water  in  the  main 
channel  and  the  cioserlftend  war  allowed  to  drain  slowly  through  the 
sampling  taps  and  the  drain  hole  on  the  bottom  of  the  headbay 
without  disturbing  the  bottom  deposits  in  the  main  and  closed-end 
channels. 

12)  About  21  hours  later,  the  main  and  closed-end  channels  were  divided 
into  several  sections  and  the  bottom  deposit  in  each  was  pumped 
into  a  plastic  jar  using  a  vacuum  pump. 

Tests  with  tidal  effect.  The  experimental  procedure  for  steady 
flow  tests  was  also  valid  for  the  tests  with  tidal  effect.  In  addi- 
tion, the  water  surface  elevation  at  the  main  channel  was  changed  by 
adjusting  the  elevation  of  the  overflow  weir  to  simulate  the  tidal 
effect  on  the  movement  of  turbidity  current  and  the  associated  sedimen* 
tation  in  the  closed-end  channel.  In  COEL-9  with  "flood  tide,"  the 
water  surface  elevation  in  main  channel  was  increased  twice  (each 
increase  by  a  step),  then  decreased  twice  (each  decrease  by  a  step),  to 
return  back  to  the  original  surface  level.  The  magnitude  of  each 
surface  elevation  change  and  the  corresponding  time  are  given  in 
Fig. 1.11.  In  COEL^IO  with  "ebb  tide,"  the  water  surface  elevation  in 
the  main  channel  was  dropped  twice  (each  drop  by  a  step),  then  increas- 
ed twice  (each  increase  by  a  step),  and  is  shown  in  Fig. Il.1t,  A 
summary  of  test  conditions  for  the  tests  performed  at  COEL  is  given  in 
Table  1.4. 
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CHAPTER  5 
RESULTS  AND  DISCUSSION 


5.1  Typical  Results 

'Typical'  refers  to  a  representative  format  of  the  results,  in 
that  results  given  in  this  section  generally  represent  directly 
measured  experimental  parameters  and  supplementary  results  generated 
using  the  numerical  model  and/or  computed  by  analytically  derived 
formulas.  These  results  are  from  an  individual  test  (C0EL*5),  which 
was  also  used  to  verify  the  calibrated  numerical  model.  A  complete 
set  of  typical  results  is  presented  to  Illustrate  the  fundamental 
mechanics  of  a  turbidity  front,  a  stationary  sediment  wedge,  and 
associated  sedimentation  in  a  closed^end  channel.  In  order  to  empha^ 
size  the  magnitude  of  these  physical  quantities  of  concern,  a  dimen<- 
sional  format  with  arithmetic  axes  is  chosen  for  the  plots. 
5.1.1  Main  Channel 

Vertical  velocity  profile.  The  purpose  of  establishing  a  steady 
flow  In  the  main  channel  is  to  keep  the  suspension  concentration 
unchanged  throughout  a  test.  The  locations  for  measuring  the  vertical 
velocity  profile  and  the  procedure  were  presented  in  Chapter  M, 
Fig. 5.1  exhibits  a  vertical  velocity  profile  measured  in  the  main 
channel  in  test  COEL^-5.  The  horizontal  velocity  in  the  main  channel 
is  denoted  as  u0>  and  the  depth-mean  velocity,  u0,  i3  assumed  to  be 
identical  to  the  horizontal  velocity  at  the  elevation  z  -  H/e,  where  e 
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is  the  base  of  natural  logarithm,  base'i  oil  ihe  logarithmic  velocity 
profile  assumption.  The  ranges  of  v^,   wert   6,0-21.0  cm/sec  in  WES 
teats,  I  nd  7.3«19.8  om/sec  in  COEL  tests,  respectively.  The  individual 
depth-mean  velocity  for  each  test  Is  given  in  Table  1.2  and  1.1. 

Vertical  sediment  concentration  profile.  The  vertical  distribu* 
tion  of  suspended  sediment  concentration  in  the  main  channel  near  the 
entrance  of  the  closed-end  channel  is  the  factor  controlling  the 
turbidity  current  and  sedimentation  in  the  closedfrend  channel.  The 
suspension  concentration,  C0,  given  in  Fig. 5. 2  was  taken  from  an 
arithmetic  mean  of  concentrations  at  locations  A  and  B  over  three 
withdrawal  times.  Fig. 5. 2  shows  an  approximately  uniform  vertical 
concentration  profile  in  test  C0EL*5.  The  depth*mean  concentrations  in 
the  main  channel,  denoted  as  C0,  for  all  tests  performed  at  WES  and 
COEL  are  given  in  Tables  1.2  and  1.1. 
5.1.2  Closed-End  Channel 

Time  variation  of  front  position.  As  mentioned  in  Chapter  2, 
three  distinct  phases  (i.e.  initial  adjustment  phase,  inertlal  self*si*> 
milar  phase,  and  viscous  self "similar  phase)  have  been  found  in  the 
front  propagation  of  salinity-induced  gravity  current  along  a  channel. 
A  plot  of  front  position,  xf,  with  the  elapsed  time,  t,  for  test  C0EL-5 
is  given  in  Fig. 5, 3.  Three  distinct  slopes  on  the  curve  are  observed. 
Among  Una,  two  slopes  at  r.he  early  stage  (slopes  of  1  and  2/3)  are 
Identical  to  the  characteristics  of  the  initial  adjustment  phase 
and  inertia!  self-similar  phase,  respectively,  of  sanility-induced 
gravity  currents  (Fig. 2. 3).  During  the  phase  of  the  initial  adjust- 
ment, the  turbidity  front  propagated  with  a  constant  speed,  which  was 
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Figure  5.1.  Main  Channel  Horizontal  Velocity  vs.  Elevation. 
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Figure  5.2.  Main  Channel  Concentration  vs.  Elevation. 
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motivated  by  the  initial  density  gradient  i.etween  the  main  and  closed" 
end  channels.  During  the  Inertial  self-.'.imilar  phase,  in  which  the 

force  balance  is  dominated  by  inertia  and  gravity  forces,  the  turbidity 

•1/3 

front  proceeded  at  a  speed  proportional  to  t    ,  where  t  is  the 

elapsed  time.  The  third  slope  observed  on  the  curve  in  Fig. 5.3  has  a 
value  of  0.29,  which  is  slightly  higher  than  the  slope  of  1/5  (the 
characteristic  value  of  the  viscous  self-similar  phase  for  salinity 
current).  It  suggests  that  during  this  period,  the  front  movement  was 
not  completely  controlled  by  gravity  and  viscous  forces,  possibly  due 
to  an  extra  force  introduced  by  a  longitudinal  temperature  gradient, 
which  was  observed  during  the  test. 

Distance  variation  of  front  speed.  The  variation  with  distance  of 
front  velocity,  Uf,  obtained  from  measured  travel  interval  and  travel 
time,  is  plotted  in  Fig.5.1).  The  data  show  that  the  front  velocity 
rapidly  decreased  as  the  front  propagated  farther  downstream  along  the 
channel.  The  deceleration  of  front  velocity  is  mainly  induced  by  a 

drag  force  per  unit  width,  which  is  against  the  front  movement  and 

2 
equal  to  O.5C0pwUfh2,  where  Cn  is  drag  coefficient  and  hg  is  the  front 

head  height.  Also,  the  motive  force  of  front  motion  per  unit  width, 

2 
0.5Ap'gh2,  decreased  with  distance  from  the  entrance  of  the  channel. 

The  reduction  of  the  motive  force  was  mainly  due  to  the  decrease  of 

front  density  by  sediment  deposition  and  interfacial  mixing.  There*- 

fore,  it  would  be  expected  that  the  rate  of  the  decrease  of  front 

speed  with  distance  of  a  turbidity  current  will  be  faster  than  that  of 

a  salinity-induced  gravity  current,  in  which  the  front  density  decrease 

es  due  only  to  interfacial  mixing. 
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Time  variation  of  hu.'-i/.ui'tdl  velocity   The  flow  velocities  in 
the  closed"Send  channel  generated  in  laboratory  tests  at  WEi  and  COEL 
were  less  than  1.0  cm/sec.  For  such  small  density  currents,  dye 
injection  technique  is  often  used  by  many  investigators  (e.g.  Sturm 
4  Kennedy,  1980;  Brocard  &  Harleman,  1980)  to  measure  the  horizontal 
velocity  profile  over  vertical.   Unlike  the  density  currents  induced  by 
salinity  and  temperature  gradients,  visibility  in  the  turbidity  current 
is  low.  Thus,  for  the  velocity  profile  measurement  in  the  present 
study,  a  vertical  streak  of  dye  was  injected  at  a  location  approximate" 
ly  1  cm  away  from  the  side  wall,  and  the  subsequent  times  of  dye 
displacement  were  recorded.  These  measured  magnitudes  of  horizontal 
velocity  were  smaller  than  the  flow  velocities  occurring  at  the  center 
of  the  channel  due  to  side  wall  resistance  (Lin  &  Mehta,  1986).  In 
addition,  due  to  the  lack  of  several  succesive  measurements  of  the 
vertical  velocity  profile,  the  temporal  and  spatial  variations  of 
horizontal  velocity  in  the  channel  are  demonstrated  by  simulation 
results  obtained  using  the  numerical  model. 

The  variation  of  the  maximum  horizontal  velocity  in  the  lower 
layer,  at  z  ■  2.5  cm,  with  the  elapsed  time  for  various  locations  in 
the  channel  is  shown  in  Fig. 5. 5.  The  velocity  was  obtained  using  the 
calibrated  numerical  model  with  the  given  conditions  of  test  C0EL>-5. 
In  Fig. 5. 5,  curves  of  horizontal  velocity  vs.  time  at  x  =  0.0  a   and  2.0 
m  show  that  the  horizontal  velocity  at  2.5  cm  elevation  increases  from 
zero  to  a  maximum  value,  then  decreases  and  eventually  reaches  <j   steady 
value.  It  was  observed  experimentally  that  the  maximum  value  (i.e.  the 
peak  of  the  curve)  occurred  at  the  moment  when  the  turbidity  front 
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passed  the  observation  point.  Note  that  the  time  required  for  the 
establishment  of  a  steady  horizontal  velocity  after  the  passing  of  the 
front  at  a  fixed  location  is  relatively  short. 

Time  variation  of  horizontal  velocity  profile.  Fig. 5. 6  exhibits 
the  vertical  velocity  profile,  generated  by  using  the  calibrated 
numerical  model,  at  x  =  2.0  m  for  four  different  elapsed  times.  In  the 
figure,  positive  velocity  Is  identified  as  being  towards  the  closed 
end.  At  times  60  minutes  and  220  minutes  after  gate  opening,  the 
occurrence  of  Identical  velocity  profile  suggests  that  a  steady 
horizontal  velocity  at  x  -  2.0  m  has  been  established.  Also,  a 
decrease  of  the  magnitude  of  horizontal  velocity  and  an  Increase  of  the 
Interface  elevation  of  the  velocity  profile  with  increasing  time  are 
observed. 

Distance  variation  of  horizontal  velocity  profile.  The  distance 
variation  of  the  steady  state  horizontal  velocity  profile  is  shown  in 
Fig. 5. 7.  These  profiles  were  produced  by  the  calibrated  numerical 
model  with  the  test  conditions  of  C0EL<*5.  It  can  been  shown  that  the 
discharge  In  the  layer  below  interface  Is  equal  to  the  discharge  in  the 
layer  above  interface  at  all  three  locations,  which  satisfies  the 
condition  of  volumetric  conservation  of  water  at  steady  state.   In 
addition,  from  the  figure  a  decrease  in  the  magnitude  of  the  horizontal 
velocity  and  of  the  interface  elevation  with  increasing  x  are  observ- 
ed. A  contour  plot  of  the  steady  state  flow  regime  in  the  channel  Is 
given  in  Fig. 5. 8. 

A  comparison  of  the  vertical  velocity  profile  among  experimental, 
numerical,  and/or  analytical  results  was  made  for  test  C0EL*6  (not 
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Figure  5.6.  Horizontal  Velocity  vs.  Elevation  at  Four  different  times 
After  Gate  Opening. 
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sufficient  measured  data  off  Horizontal  velocity  in  COEL-S  test).  These 
rerults  are  presented  in  9  dimensionlese  format  and  discussed  in 
Appendix  B.3. 

Time  variation  of  suspension  concentration.  As  mentioned  in 
section  3 • 2 . 4 ,  test  C0E1M  was  selected  as  a  calibration  test  to 
calibrate  five  unknown  coefficients  (i.e.,  horizontal  and  vertical 
momentum  and  mass  diffusion  coefficients  ex,  ez,  ex,  ez,  and  Manning's 
n),  in  the  numerical  model.  After  comparing  the  numerical  results  with 
the  corresponding  measurements  in  COEL-I,  these  calibrated  values  for 
the  bestf-f itting  comparisons  were  found  to  be  the  following:  ex  =  £x 
■  5.0  x  10^3  m2/sec,  ez  »  ez  «  10**  m2/sec,  and  n  =  0.03.  Also  as 
noted,  the  values  of  Ax  ■  0.5  m,  Az  ■  1.0  cm,  and  At  -  0.15  sec,  were 
adopted  in  the  numerical  simulation.  Fig. 5. 9  shows  the  calibration 
result  of  the  time  variation  of  suspension  concentration  at  z  ■  0.7  cm 
for  various  locations  in  the  closed^end  channel.  In  each  test,  the 
suspension  concentration  profile  at  each  sampling  location  was  general'* 
ly  measured  at  three  different  times,  and  the  first  of  the  three 
measurements  was  taken  during  the  passage  of  the  turbidity  front.  The 
symbols  given  in  Fig. 5. 9  represent  the  measured  data  at  different  times 
and  various  locations,  while  the  lines  represent  the  calculations  by 
using  the  numerical  model  with  the  aforementioned  coefficients. 

The  calibrated  numerical  model  was  verified  using  the  measured 
data  of  C0EL-5,  and  the  verification  result  is  shown  in  Fig. 5. 10.   In 
Fig. 5. 10,  the  variation  of  bottom  suspension  concentration  (at  z  ■  0.7 
cm)  with  elapsed  time  at  various  locations  in  the  channel  is  given. 
The  symbols  in  Fig. 5. 10  are  measurements  at  different  times  and  various 
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locations,  and  the  lines  indicate  predicted  results  using  the  calibraf- 
ed  numerical  model  at  the  corresponding  locations.  A  reasonable 
agreement  between  the  measured  data  and  numerical  results  is  observed. 
It  is  also  noted  that  the  suspension  concentration  at  a  particular 
location  during  the  front  passage  was  lower  than  that  at  the  steady 
state. 

Distance  variation  of  concentration  profile.  Fig. 5. 11  shows 
the  calibration  result  of  the  distance  variation  of  the  vertical 
concentration  profile  in  the  channel  for  test  COEL-M.  The  vertical 
concentration  profiles  measured  at  86  minutes  and  130  minutes  after- 
gate  opening  at  various  sampling  locations  as  well  as  profiles  generat*- 
ed  by  the  numerical  model  (with  the  calibrated  coefficients)  at  the 
corresponding  times  and  locations  are  shown  in  Fig.5.11. 

Fig. 5. 12  shows  model  verification  of  the  distance  variation  of 
vertical  concentration  profile  in  the  channel  for  COEL"5.  The  vertical 
concentration  profiles  measured  at  67  minutes  and  190  minutes  after 
gate  opening  at  various  sampling  locations,  as  well  as  the  correspond- 
ing numerical  results,  are  shown  In  Fig. 5. 12.  A  good  agreement  between 
the  measured  concentration  profile  and  numerical  prediction  is  observed 
for  the  locations  of  x  =  0.1,  2.1,  and  1.6  m.  But,  at  locations  x  = 
7.9  and  11.0  m,  the  numerical  results  are  observed  to  give  slightly 
lower  values  than  those  of  measured  data.  This  discrepancy  may  be 
induced  by  using  relatively  high  settling  velocity  in  the  numerical 
model  for  the  location  near  the  end  of  the  channel,  where  the  sediment 
settling  velocity  could  be  smaller  due  to  the  sediment  sorting  effect 
of  a  nonuniform  flow  regime,  as  well  as  measurement  errors  involved 
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Figure  5.11.  Concentration  at  Steady  State  vs.  Elevation  at 
Five  Locations.  (Calibration  Results) 


113 


12 


10 


67min 

I90min 

i 

a 

■      C 

0.1m 

~       & 

•     D 

2.1m 

o 

•     E 

4.6m 

o 

•     F 

7.9m 

♦     G 

II.  Om 

V 

COEL- 5 


Numerical 
'  Results 


200  400  600 

CONCENTRATION,  C  (mq//) 


800 


Figure  5.12.  Concentration  at  Steady  State  vs.  Elevation  for 
Five  Locations.  (Verification  Results) 
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in  the  gravimetric  procedure  for  low  sediment  concentration.  The 
concentration  profiles  at  times  67  m  nutw  and  190  minutes  obtainec 
by  using  the  numerical  model  are  identical  at  various  locations, 
suggesting  that  the  steady  state  concentration  profile  has  been 
established  67  minutes  after  the  gate  opened.  Yet,  the  measured 
concentration  profiles  at  the  time  67  minutes  and  190  minutes  are  not 
identical,  as  shown  in  Fig.5.12.  This  is  mainly  due  to  a  slightly 
higher  boundary  concentration  at  main  channel  introduced  during  the 
period  between  67  minutes  and  190  minutes  after  gate  opening  and  due  as 
well  to  the  analysis  errors  of  the  collected  liquid  samples. 

Distance  variation  of  concentration  at  front  head.  As  mentioned 
before,  the  decrease  of  the  sediment  concentration  at  the  turbidity 
front  with  distance  is  induced  by  sediment  deposition  and  to  some 
extent  by  interfacial  mixing.  Fig.5.13  shows  the  distance  variation  of 
sediment  concentration  at  front  head  in  C0EL-5.  The  data  point  in  the 
figure  represents  the  suspension  concentration  at  z  =  0.7  cm  (i.e.  the 
elevation  of  the  lowest  tap  at  each  sampling  location),  during  the 
passage  of  the  turbidity  front.  Note  that  no  front  concentration  was 
taken  from  section  C  (the  sampling  location  at  x  =  0.1  a   of  the 
channel),  because,  there,  the  turbidity  intrusion  from  the  main  channel 
was  driven  by  shear-induced  circulation  (within  the  gyre  zone  as  shown 
in  Fig. 1.1).  According  to  Keulegar.'s  ('957)  investigation  on  the  lock 
exchange  flows  of  saline  water,  front  dilution  is  more  pronounced  in 
the  lock  flows  having  relatively  small  Reynolds  number.  Nevertheless, 
there  are  no  adequate  experimental  data  from  his  investigation  that  can 
be  used  to  compare  with  the  result  shown  in  Fig.5.13.   Since  the  front 
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dilution  of  salinity-induced  density  current,  is  only  affected  by 
interfaciai  mixing,  it  should  be  expected  that  the  decrease  of  front 
density  with  distance  at  the  turbidity  front  will  be  faster  than  that 
of  a  salinity  front  under  the  same  initial  conditions.  Fig. 5. 13 
exhibits  a  rapid  decrease  of  suspension  concentration  at  the  turbidity 
front  with  distance. 

Distance  variation  of  mean  concentration  below  interface.  Under  a 
steady  state,  longitudinal  variations  of  mean  concentration  below  the 
Interface,  C^,  Is  given  in  Fig. 5. 11  and  Fig. 5. 15,  where  the  interface 
is  defined  as  the  zerof-veloclty  elevation.  Fig. 5. 11  shows  the  calibra^ 
tlon  result,  including  measured  data  and  numerical  results  for  COEL-t. 
The  measured  mean  concentration  below  the  interface  at  the  sampling 
location  was  obtained  by  averaging  the  mean  concentration  below  the 
interface  at  86  minutes  and  130  minutes  after  gate  opening.  Fig. 5. 15 
shows  the  verification  result  by  comparing  the  measured  data  and  the 
corresponding  numerical  result  for  C0ELf5,  and  In  addition,  includes 
analytical  result  obtained  by  using  eq.3.52.   It  can  be  seen  from 
Fig. 5. 15  that  the  measured  mean  concentration  below  the  interface 
(obtained  by  averaging  the  concentrations  at  times  67  minutes  and  190 
minutes)  exponentially  decreases  with  distance  from  the  entrance.  The 
decrease  of  mean  concentration  below  the  interface  along  the  channel 
resulted  from  the  deposition  of  suspended  sediment  and  the  interfacial 
mixing  of  tne  sediment  through  the  interface.  Also,  the  rapidly 
decreasing  concentration  suggests  that  the  assumption  of  constant 
sediment  concentration  along  the  channel,  made  by  McDowell  (19711  in 
his  theoretical  study  on  the  settling  solidsMnduced  currents,  may  be 
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inappropriate  for  solids  having  high  settling  velocity.  The  lengitudi" 
nai  distribution  of  mean  concentration  below  the  interface  obtained  by 
the  numerical  model  and  the  analytical  approach  are  also  shown  in 
Fig. 5. 15.  The  calibrated  coefficient  6  -  7.C  in  eq.3.52  indicates  that 
the  stationary  sediment  wedge  had  a  velocity  at  entrance  u-|  -  0.57  Ufi , 
which  is  close  to  the  experimental  result,  u-|  =  0.70  Uf|  ,  found  by  Gole 
et  al.  (1973).  Comparison  shows  a  good  agreement  between  the  numerical 
and  theoretical  results  with  the  measured  data.  In  addition,  a  plot  of 
concentration  contour,  obtained  from  the  results  of  numerical  Simula^ 
tlon,  for  test  C0EL>-5  at  steady  state  is  shown  In  Fig.5.16. 

Longitudinal  variation  of  water  surface.  The  longitudinal  varia** 
tion  of  the  water  surface  in  the  channel  discussed  in  this  section  is 
represented  by  the  elevation  difference,  AH(x),  between  local  water 
depth  and  the  water  depth  at  the  entrance.  Since  the  magnitude  of  AH 
was  on  the  order  of  0.01  mm  in  all  tests,  the  measurement  of  AH  using  a 
point  gauge  was  out  of  the  question.  Instead,  the  longitudinal 
variations  of  AH  in  the  channel  for  C0ELU5  were  simulated  using  the 
numerical  model,  as  shown  In  Fig. 5.17.  The  longitudinal  variations  of 
AH  at  5  minutes  and  15  minutes  may  be  referred  to  as  the  unsteady 
cases,  and  at  times  90  minutes  and  180  minutes  (identical  curve)  may  be 
considered  as  the  steady  state  case.   During  unsteady  state,  the 
propagation  of  small  surface  waves,  induced  by  the  intrusion  of  dense 
fluid  into  the  clear  water,  can  be  seen  in  the  figure.  At  steady 
state,  the  water  surface  rises  from  the  entrance  towards  the  closed 
end,  and  reaches  a  horizontal  level  at  the  location  beyond  the  toe  of  a 
stationary  sediment  wedge. 
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The  maximum  rise  of  water  3urface,  AHroax,  can  be  computed  accord* 
ing  to  the  analytical  derivation  (eq.3.6';)  with  given  data  on  H,  pH, 
and  p-,u).   in  test  COEL*i.  H-9.3  cm,  pw  -  0.99683b  g/cm3,  and  pi(s), 
the  density  profile  at  entrance,  was  measured.  Substituting  these 
known  vaJues  Into  eq.3.61,  AHmax  is  found  to  be  0.024  mm,  which  is 
close  to  the  numerical  result  of  0.018  mm.  The  good  agreement  between 
the  theoretical  and  numerical  results  of  AHmax  confirms  that  a  station'- 
ary  sediment  wedge  can  be  established  in  a  closed«end  channel  in  the 
steady  state  if  provided  the  entrance  boundary  conditions  remain 
unchanged.   Detailed  calculation  of  AHmax  using  eq.3.64  for  COEL^S 
is  given  in  Appendix  C. 

Distance  variation  of  deposition  rate.  At  about  2t  hours  after 
the  end  of  each  test,  the  bed  deposit  collected  from  each  subsection 
(of  selected  length  depending  upon  the  total  length  over  which  measure 
able  sediment  mass  was  deposited)  was  analyzed  for  the  deposited  mass 
and  particle  size  distribution.  Mean  sediment  deposition  rate  was 
computed  by  dividing  the  deposited  mass  by  the  subsection  area  and 
effective  duration.  Fig. 5. 18  shows  the  calibration  result  of  longitu- 
dinal distribution  of  mean  sediment  deposition  rate.  The  dashed  line 
in  Fig. 5. 18  represents  the  distance  variation  of  mean  sediment  deposi- 
tion rate,  which  was  obtained  from  test  COELM.  The  solid  line  in 
Fig.5.l8  shows  the  corresponding  deposition  rate  distribution  using  the 
numerical  model. 

The  verification  result  of  the  numerical  model  is  shown  in 
Fig. 5. 19.  In  this  figure,  the  dashed  line  shows  the  distance  variation 
of  mean  sediment  deposition  rate  of  C0EL-5.  Typically,  there  was 
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relatively  very  high  deposition  in  the  zone  Inside  the  entrance, 
followed  by  an  approximately  e/ponentialL-type  decrease  of  deposition 
proceeding  farther  along  the  channel.  The  heavy  deposition  in  the  zone 
near  the  entrance  occurred  when  suspended  particles,  which  were 
aggregated  under  a  relatively  strong  turbulent  flow  in  the  main 
channel,  moved  into  the  relatively  quiescent  closed"-end  channel  (most 
of  them  being  heavier  than  the  current  could  carry),  and  deposited  on 
the  bed  almost  immediately.  The  solid  line  in  Fig. 5. 19  shows  the 
prediction  of  mean  sediment  deposition  rate  of  COEL/'S  using  the 
numerical  model.  A  reasonable  agreement  between  the  numerical  simula- 
tion with  the  experimental  results  is  observed. 

Dispersed  particle  size  distribution  of  deposit.  Fig. 5. 20  shows 
the  distance  variation  of  the  dispersed  particle  size  distribution  of 
the  deposit.  The  particle  size  distribution  was  obtained  by  hydrometer 
analysis  for  each  bed  deposit  sample.  The  dispersed  particle  size 
shown  in  the  figure  actually  represents  the  primary  particle  size, 
i.e.,  particles  from  which  floes  are  formed  during  aggregation.  There^- 
fore,  the  particle  sizes  in  the  figure  should  not  be  interpreted  as  the 
size  of  floes  depositing  from  suspension  during  the  test. 

In  Fig. 5. 20,  it  can  be  seen  that  the  median  particle  diameter  of 
deposit  decreased  with  distance  from  the  entrance.  Note  that  the 
deposit  sample  in  the  main  channel  had  the  maximum  median  particle 
diameter,  dm  ■  3.9  urn.  Three  particle  size  distribution  curves  shown 

in  the  figure  are  more  or  Jess  parallel  to  each  other.  A  modified 

1  /? 
sorting  coefficient  ,  which  may  be  defined  as  (dg5/dm)    for  the 

present  study,  can  be  used  as  an  index  representing  the  degree  of 


12C. 


c 

o 

m 

~< 

O 

IS 

o 

a 

01 

I 

h 

P 

, ^ 

jj 

E 

M 

a. 

03 

(0 

CO 

z 

0) 

o 

tr 

%4 

m   o 

O 

-Q  5 


M 


-8 


< 

CD 


(%)    1H9I3M    A9    JQNIJ  !N33U3d 


■  .. 


127 
uniformity  of  particle  size  distributiv"  of  the  deposit.  The  modified 
sorting  coefficients  for  all  three  cur.''  towards  the  smallest  particle 
size  3hown  In  Fig. 5. 20  are  2.4,  2.2  and  1.9,  respectively.  This  trend 
of  relatively  constant  sorting  coefficients  was  found  to  he  generally 
true  for  all  the  tests.  The  median  dispersed  particle  size  of  the 
deposit  and  the  modified  sorting  coefficient  of  each  particle  size 
distribution  for  all  COEL  tests  are  given  in  Appendix  D. 

5.2  Characteristic  Test  Parameters 
Tables  5.1  and  5.2  summarize  the  characteristic  parameters  for  all 
of  the  tests  performed  at  WES  and  COEL,  respectively.  The  parameters 
included  in  these  tables  are  either  characteristic  In  the  sense  that 
they  are  representative  factors  in  a  physical  sense  of  a  test,  or 
characteristic  in  the  sense  that  they  are  frequently  used  for  computing 
or  normalizing  the  variables  of  interest  in  the  subsequent  sections  of 
this  chapter. 

5.3  Characteristic  Results 
The  results  of  steady  flow  tests  represented  in  this  section  are 
considered  to  be  characteristic  because  they  generally  include  the  data 
from  WES  and  COEL  te3ts.  Most  study  results  are  presented  in  a 
dimensionless  or  normalized  format,  while  some  results,  in  which  the 
characteristics  of  significant  physical  factors  are  examined,  are 
presented  in  a  dimensional  format.   The  results  of  tidal  effect  tests 
conducted  at  COEL  are  included  and  discussed  in  the  subsequent  section. 
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5.3-1  Front  Behavior 

Front  nose  height,  to  head  height  ratio.   According  to  dimensional 
analysis  (eq.3.2),  the  ratio  of  front  nose  height  to  head  height, 
hj/hj,  may  be  dependent  upon  the  local  head  Reynolds  number,  Ufhj/v, 
and  the  fraction  depth,  h2/H.  However,  it  has  been  pointed  out  by 
Lawson  (1971)  that  h^/fig  is  independent  of  the  fraction  depth  h2/H. 
In  order  to  compare  with  previous  investigation  results,  the  effect  of 
the  local  head  Reynolds  number,  Ren  =  Ufhj/v,  on  h^/h2  is  examined  and 
shown  In  Fig. 5. 21,  even  though  both  examining  dimensionless  parameters 
Ufh2/v  and  h-|/h2  include  the  same  parameter  h2.  The  data  points  of 
present  study  are  obtained  from  all  COEL  tests.  A  trend  of  decreasing 
h-|/h2  with  increasing  Ren  is  evident.  The  result  indicates  that  with 
higher  relative  viscous  force  In  the  flow  (i.e.  low  Re^,  value),  the 
ratio  of  hj/112  will  be  higher,  as  expected.   In  Fig. 5. 21,  in  addition 
to  the  results  from  the  present  study,  the  data  deduced  from  Braucher 
(1950),  Keulegan  (1958),  Wood  (1965),  Lawson  (1971),  and  Simpson  (1972) 
are  included.  The  atmospheric  results  obtained  by  Lawson  (1971)  show 
considerable  scatter,  and  only  the  upper  and  lower  bounds  of  the  data 
points  are  shown  in  the  figure.  The  consistent  trend  given  by  all  of 
the  results  suggests  that  the  relationship  between  h-|/h2  and  Ren  shown 
in  this  figure  is  valid  for  all  types  of  gravity  currents. 

Front  head  height  to  neck  height  ratio.  According  to  dimensional 
analysis  (eq.3.5),  the  ratio  offff  front  head  height  to  neck  height, 
(12/113,  may  be  dependent  upon  the  local  neck  Reynolds  number,  Ren  ■ 
Ufh^/v,  the  fraction  depth,  hj/H,  and  the  local  relative  settling 
velocity,  ws/uf.  Fig. 5. 22  represents  the  effects  of  the  local  neck 
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Reynolds  number,  Ren.  and  the  fraction  depth,  t^/H,  on  the  ratio  of 
i12/'h3 ■  Tne  data  wider  consideration  here  were  obtained  from  all  COEL 
teats  over  the  range3  of  1en   «  20*1,000  and  the  fraction  depth,  ho/H  - 
0.19*0.39.  A  relatively  constant  value  of  h2/ti3,  with  an  arithmetic 
mean  of  1.3,  Is  obtained  over  the  entire  Ren  range.  Keulegan  (1958) 
also  found  a  relatively  constant  value  of  h2/h3,  with  an  arithmetic 
mean  2.08,  over  ranges  of  Ren  -  3,000«30,000  and  I13/H  -  0.1»0.21. 
Based  on  experimental  observations,  Simpson  and  Britter  (1979)  found 
that  the  ratio  of  112/113  decreases  with  increasing  fraction  depth  I13/H. 
Therefore,  the  different  constant  values  of  h2/ti3,  as  shown  In  Figure 
5.22,  obtained  in  the  present  study  and  by  Keulegan  (1958),  should  have 
resulted  from  different  fraction  depth  in  both  studies.  For  a  constant 
113/11,  the  ratio  of  113/113  is  nearly  independent  of  the  local  neck 
Reynolds  number  Ren.  Finally,  the  relatively  constant  ratio  of  h2/h3 
In  the  present  study  suggests  h2/h3  Is  independent  of  the  relative 
settling  velocity,  ws/uf. 

Dlmenslonless  initial  front  3peed.  According  to  the  results  of 
dimensional  analysis  (eq.3.8),  the  ratio  of  initial  front  speed  to 
densimetric  velocity,  Ufj/u^,  may  be  dependent  upon  the  denslmetric 
Reynolds  number,  ReA  -  u4H/\>,  and  the  wldth»to«depth  ratio,  B/H, 
Fig. 5. 23  shows  the  effect  of  the  denslmetric  Reynolds  number,  Re6  - 
uAH/v,  and  the  width^tocdepth  ratio,  B/H,  on  the  ratio  of  Ufi/uA.  Tne 
data  points  from  the  tests  of  WES  and  COEL,  are  distinguished  by  using 
different  symbols.  The  width»to*depth  ratio,  B/H,  of  WES  tests  varied 
from  2.5  to  1.6,  and  the  B/H  of  COEL  test  was  approximately  1.0.  There 
exists  a  relatively  great  scatter  of  data  points,  which  Is  believed  to 
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Figure  5.23.  Initial  Front  Speed  to  Densiraetrio  Velocity  Ratio 
vs.  Denaimetrie  Reynolds  Number. 
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be  due  to  the  presence  of  sneor*Hnduced  gyre  at  the  entrance.  The  gyre 
length  in  x  direction  was  lr.  the  range  of  0.2*0.1  ■  in  WES  tests  and 
0.1(40.2  m  in  COEL  tests.  The  measured  front  speed  at  x  -  1.0  ra  for  WES 
tests,  and  x  -  0.75  m  for  COEL  teats  was  taken  as  the  initial  front 
speed,  which  was  defined  as  the  front  velocity  induced  due  to  the 
Initial  density  gradient.  The  initial  front  speed  of  each  test  is 
listed  in  Tables  5.1  and  5.2. 

From  Fig. 5. 23,  no  measurable  difference  between  the  dimensionless 
initial  front  speeds  obtained  from  WES  and  from  COEL  tests  is  observe 
ed.  This  result  may  suggest  that  the  ratio  of  ufl/uA  is  independent  of 
the  width*to»depth  ratio,  B/H.  The  arithmetic  mean  value  of  0.13  for 
Uf|/u^  Is  obtained  from  all  data  points  and  is  shown  as  a  dashed  line. 
Due  to  the  relatively  large  data  scatter,  It  is  difficult  to  conclude 
whether  or  not  a  clear  trend  with  respect  to  ReA  or  B/H  exists.   In 
previous  investigations  on  the  lock  exchange  flow,  the  ratio  of  uf1/uA 
was  found  to  be  Independent  of  the  densimetric  Reynolds  number  and  the 
width*tO"depth  ratio,  and  equal  to  0.5  from  theoretical  development  by 
Yih  (1965).  Keulegan  (1957)  obtained  an  approximately  constant 
experimental  value  of  0.16  for  ufl/uA  over  the  range  of  Re&  -  100* 
20,000  and  B/H  -  0.11*1.0.  Comparing  the  ratio,  Uf^/u^,  of  0.13 
obtained  from  the  present  study  with  keulegan 's  0.16,  the  agreement, 
notwithstanding  the  presence  of  the  gyre  at  the  entrance  of  the 
closedtiend  channel,  suggests  that  the  front  motion  beyond  x  -  1.0  m  in 
the  WES  channel  and  x  -  0.75  m  in  the  COEL  channel  towards  the  closed 
end  of  the  channels  was  mainly  controlled  by  the  turbidity  gradient 
with  the  unchanged  concentration  source  in  the  main  channel. 
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Distance  variation  of  normalized  front  speed.   According  to  the 
result  of  dimensional  anal>sis  !eq.3.10),  the  dlmensionless  front 
speed,  ur/uA,  was  found  to  be  dependent  upon  the  front  position,  x/H, 
the  wldth*tortdepth  ratio,  B/H,  the  denslmetrlc  Reynolds  number,  ReA  - 
uAH/v,  and  relative  settling  velocity,  w3/uA.  Keulegan  (1957,  1958) 
examined  the  effects  of  ReA  and  B/H  on  the  distance  variation  of 
dlmenslonless  front  speed  In  lock  exchange  flow.  He  proposed  the 
following  relationship  between  Uf/uA  and  x/H 

uA   uA     x 

—  -  —  +  T-  (  5.1  ) 

U{       Uf|    H 

where  the  parameter  y   is  a  function  of  ReA  and  B/H.  The  relationship 
between  uf/uA  and  x/H  given  in  eq.5.1  for  the  salinlty»lnduced  density 
current  was  found  not  to  be  valid  for  turbidity  current  due  to  the 
Influence  of  the  local  relative  settling  velocity,  w3/uA.  According  to 
the  results  of  a  large-scale  flume  study,  Barr  (1967)  suggested  that 
front  movement  of  gravity  current  is  Independent  of  ReA  if  ReA  >  1,000. 

In  the  present  study,  a  normalized  front  speed,  Uf/Tjfi ,  was 
adopted  insteady  of  uf/uA,  while  noting  that  the  other  influential 
dimensionles3  parameters  should  be  the  same.  Here,  uf1  is  the  average 
front  speed  over  a  distance  0.1A2.0  m  from  the  entrance  of  the  channel 
for  COEL  tests,  (witnln  >:  .  0.9*0.4  m,  the  front  speed  was  influenced 
by  shear*induced  circulation).  Also,  because  the  settling  velocity  of 
suspension  Is  a  distancenvarying  parameter,  the  settling  velocity  of 
suspension  at  the  entrance  of  the  closed-end  channel,  w3, ,  may  be  used 
as  the  reference  parameter  lr.  the  present  and  subsequent  discussions. 
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Fig. 5. 21  shows  the  longitudinal  variation  Of  diroer.sicnle  s  front  3peed, 
uf /uf1 >  "ith  reference  to  the  values  of  w8l 'u^   ami  He4.   It  can  be  3een 
that  the  Re^  value.?  of  all  five  tests  shown  in  the  figure  are  greater 
that,  1,000,  Accordingly,  the  distance  variation  of  front  speed  shown 
in  Fig. 5. 21  is  Independent  of  Re^  cased  on  3arr's  finding.  The  curves 
in  the  figure  represent  smoothed  trend  lines,  which  are  obtained  from 
the  corresponding  dimensionles3  plot  for  each  test. 

Dimensional  analysis  indicates  the  relative  settling  velocity, 
ws1/uA,  is  another  significant  parameter  Influencing  the  front  speed  of 
turbidity  current  (eq.3.10).  However,  the  parameter  w3l/u4,  a  refer" 
ence  value  at  the  entrance,  cannot  uniquely  describe  the  local  relative 
settling  velocity,  w3/uit  for  different  types  of  sediments  (cohesion* 
less  and  cohesive  sediments),  due  to  different  compositions  of  the 
settling  particles  of  coheslonless  and  cohesive  sediments.  Therefore, 
at  present,  the  examination  of  the  variation  of  curve  3lope  with  the 
influencing  parameter,  wal/uA,  must  be  considered  separately  for 
cohesive  kaollnite  (C0ELB2,  3,  and  5)  and  cohesionless  flyash  III 
(COEL.^6  &  7).  For  each  sediment,  the  slope  of  the  curve  in  Fig. 5. 21 
becomes  steeper  as  the  relative  settling  velocity  increases,  as  could 
be  expected.  The  break  in  slope  observed  in  the  curves  of  C0EL"2,  3 
and  6,  is  believed  to  be  due  to  the  phase  shift  from  lnertlal  selfSsl* 
milar  phase  to  viscous  self "similar  phase. 

Local  denslmetrlc  Froude  number.  According  to  the  dimensional 

analysis  (eq.3.11),  the  local  denslmetrlc  Froude  number,  FrA  =  uf/ 

1/2 
[(Ap'/pw)gh3]   ,  could  be  expected  to  be  a  function  of  local  neck 


Reynolds  number,  Ren  •  ufi^/u,  internal  stratification  of  ooncentrati 
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at  the  front,  fraction  depth,  IWH,  and  width-1?;  Cdepth  ratio,  B/H.  For 
al.l  COEL  tests,  B/H  was  kept  approximately  at  1.0.  Also,  the  fraction 
depth,  ho/H,  of  data  selected  for  examining  the  trend  of  Fr^  was  in  the 
range  between  0.19^0.39,  which  has  but  little  effect  on  Fr^  according 
to  Simpson  and  Brltter  (1979).  Therefore,  an  examination  of  the  effect 
of  Ren  only  on  Fr^  may  be  made. 

Fig. 5. 25  shows  the  effect  of  local  neck  Reynolds  number,  Ren  = 
Ufhj/v,  on  local  densimetric  Froude  number,  FrA.  The  data  points  shown 
In  the  figure  are  from  all  COEL  tests.  Although  there  is  considerable 
data  scatter,  a  trend  of  increasing  FrA  with  increasing  Ren  In  the 
range  of  Ren  -  25^320  is  evident.  Keulegan  (1958)  investigated  similar 
variables  for  experiments  with  a  saline  water  front  in  lock  exchange 
flow  and  found  that  FrA  -  0.22  Ren  ,  for  Ren  <  100  and  Fr4  -  1.05  for 
Ren  >  400,  which  is  plotted  as  a  solid  line  In  the  figure.  The  local 
densimetric  Froude  number  in  the  present  study  is  slightly  lower  than 
that  obtained  by  Keulegan  (1957)  for  the  same  Ren  value.  This  may  be 
due  either  to  the  concentration  stratification  effect  of  the  front,  or 
to  the  use  of  sediment  concentration  at  the  lowest  tap  (at  z  -  0.7  cm) 
of  the  sampling  location  to  represent  the  mean  concentration  in  the 
front.  Either  of  these  effects  will  result  In  the  smaller  FrA  values. 
In  fact,  the  ratio  of  Uf|/uA  noted  in  Fig. 5. 23  is  the  densimetric 
Froude  number  at  the  entrance  and  is  a  special  value  of  Frir  as 
dlcussed  above. 

Dimenslonless  front  position.  Based  on  dimensional  analysis 
(eq.3.16),  the  dimenslonless  front  position,  xf/H,  was  found  to  be 
dependent  upon  the  wldth»to»»depth  ratio,  B/H,  relative  settling 
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velocity,  x's/iiA,  densimetrl-.  r!e»:iolds  number,  Re^  -  uAH7v,  Mid  d-mf-i  • 
sionless  elapsed  time,  uAt/H.   In  Barr's  '1967)  large-scale  flu*?  study 
of  look  exchange  flow,  he  found  that  the  front  movement  of  gravity 
current  was  strongly  dependent  on  Re&  when  ReA  <  1,000,  and  was  almost 
Independent  of  Re&  when  ReA  >  1,000.  For  all  C0EL  tests,  Re4  of 
each  test  was  greater  than  1  ,000  (which  can  be  seen  In  Table  5.2). 
Thus  the  front  movement  for  each  test  would  be  expected  to  be  independ* 
ent  of  ReA.  Also,  all  C0EL  tests  had  approximately  the  same  wldth^to" 
depth  ratio,  B/H.  Fig. 5. 26  shows  the  variation  of  xf/H  with  uAt/H  for 
the  steady  C0EL  tests  with  given  reference  parameter,  wsi/uA.  For 
tests  using  kaollnlte,  (C0ELH2,  3,  5,  and  11),  although  there  exists 
considerable  data  scatter,  a  trend  curve,  which  comprises  three 
distinct  slopes  of  1,  2/3,  and  1/5,  can  be  discerned  in  the  figure. 
However,  for  flyash  III  tests  (COEL*"  6,  7,  and  8),  there  are  only  two 
slopes  (2/3  and  1/5)  in  the  trend  curve  of  each  test.  The  scatter  of 
data  points  from  kaolinite  tests  may  be  due  to  different  ws1/uA,  but, 
as  evident,  the  dependence  of  Xf/H  versus  u^t/H  relationship  within 
inertial  self«-similar  and  viscous  self«slmllar  phases  on  ws-|/uA  seems 
to  be  slight.  Also,  it  is  observed  that  the  duration  of  front  propagat- 
ion in  the  initial  adjustment  phase  3eems  longer  for  the  front  with 
smaller  ws1/uA  than  for  that  with  higher  w^/u^.  For  flyash  III  tests, 
no  Initial  adjustment  phase  (in  which  the  front  propagates  at  a 
constant  speed)  was  round  due  to  the  high  w3i/uA  values.  Also,  the 
duration  of  front  propagation  m  the  inertial  self ''similar  phase  (slope 
«  2/3)  was  longer  for  the  front  with  smaller  ws-|/uA  than  for  that  with 
greater  wg1/uA. 
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Figure  5.26.   Dlraenslonless  Front  Position  vs.   Dlmenslonless 
Elapsed  Time. 
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Distance  variation  of  normalized  excess  front  density,  fig. 5. 27 
shows  tne  distance  variation  of  the  normalized  excess  front  density, 

Ap'/Apd,  where  Ap'  and  Ap<j  are  the  excess  front  densities  at  z~  0-7  cm 
at  each  sampling  location  and  x-2.1  in  ,  respectively.  It  can  oe  seen 
that  Ap'/Apd  of  each  test  decreased  approximately  exponentially  with 
distance  from  the  location  of  x/H  -  20.  It  has  been  pointed  out  by 
Keulegan  (1957)  that  the  front  dilution  of  salinity-induced  density 
current  through  lnterfaclal  mixing  la  faster  In  lock  flows  having 
relatively  small  Reynolds  numbers.   Due  to  the  decrease  of  front 
density  of  turbidity  current  with  distance  depending  on  both  sediment 
deposition  and  lnterfaclal  mixing,  the  change  of  the  curve  slopes  shown 
in  Fig. 5.27  may  be  examined  with  respect  to  two  dimenslonless  parame^ 
ters,  ws-|/ua  and  u^H/v.  The  ranges  of  ws-|/u^  and  u^H/v  of  these  tests 
shown  in  the  figure  were  0.0023«0.01 72  and  1 ,260«2,630,  respectively. 
It  13  observed  that  the  change  of  curve  slope  is  consistent  with  the 
change  of  w^/u^,  which  suggests  that  the  dependence  of  the  decrease  of 
excess  front  density  with  distance  on  the  sediment  settling  is  greater 
than  on  lnterfaclal  mixing  under  the  specific  ranges  of  w^/u^  and 
u^H/v.  Also,  the  slopes  of  the  curves  in  Fig. 5. 27  become  steeper  with 
increasing  w3l/u^,  as  one  would  expect. 

Distance  variation  of  mean  concentration  below  Interface.  Figures 
5.28(a)  and  (b)  3how  the  longitudinal  distribution  of  dimenslonless 
mean  concentration  below  the  interface  for  kaolinite  and  flyash, 
respectively.  Test  number,  the  dimenslonless  relative  settling 
velocity,  w3i/uA,  and  sediment  in  each  test  are  given  in  the  figure. 
In  general,  It  can  be  seen  that  the  dimenslonless  mean  concentration 
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below  Interface,  %/C^   ,  decreases  exponentially  with  distance  from  the 
entrance.  A  bestifittlng  straight  line  for  one  test  or  a  group  of 
tests  Is  given  in  both  figures.   It  Is  observed  that  the  slopes  of 
these  straight  lines  become  3teeper  with  Increasing  Wyj/u^,  regardless 
of  the  type  of  sediment.  Furthermore,  the  numerical  value  given  for 
each  straight  line  is  the  value  of  Wgi/uj  computed  by  using  eq.3.52 
with  B  -  7.0  (the  value  used  for  the  analytic  prediction  in  Fig. 5. 15). 
Comparing  the  analytically  computed  w^/u^  value  with  experimental 
ws1^uA  value,  a  good  agreement  Is  obtained  for  smaller  experimental 
ratios  of  Wg-j/u^.  However,  when  Wg^/u^  becomes  greater,  the  disore* 
pancy  between  the  analytic  and  experimental  results  becomes  larger.  It 
can  be  Interpreted  that  with  high  ratio  of  wsl/ua,  the  local  relative 
settling  velocity,  w3/u,  will  become  much  different  from  the  ratio  of 
relative  settling  velocity  at  the  entrance,  ws-|/ui,  where  u  and  u-|  are 
the  local  and  entrance  particle  velocities  of  the  fluid,  respectively. 
Consequently,  the  term  of  exp[(B2"Bi )x]  in  eq.3.51  cannot  be  reduced  to 
be  [1+(B2"8i )x],  as  was  done,  due  to  the  significant  truncation  error 
introduced  by  omitting  the  higher  order  expansion  terms. 
5.3.2  Deposition  Related  Characteristics 

Median  column  settling  velocity.  The  results  of  tests  for 
obtaining  median  settling  velocities,  w3ra,  of  raw  sediments  in  settling 
columns  are  discussed  In  this  section.  Flg.t.6  shows  that  the  settling 
velocities  of  flyashes  I,  II,  III,  and  silica  flcur  were  not  entirely 
independent  of  the  initial  column  concentration,  C,  but  slightly 
increased  with  increasing  concentration.   A  similar  trend  was  found  for 


1 48 
Vieksburg  loess.  The  mecliaa  settling  velocity  of  Cedar  Key  mud  was 
found  to  be  0.001  cm/sec  for  C  ■  323*373  mg/1. 

Settling  column  tests  w  re  performed  for  kaollnite  with  the 
initial  suspension  concentration  varing  from  M13  mg/1  to  10,027  mg/1. 
Two  distinct  relationships  for  the  median  settling  velocity  with 
initial  column  concentration  were  found.   In  the  range  of  C  -  '443--1 ,500 
mg/1,  wsal  of  kaoiinite  Increased  with  increasing  C,  according  to 

wara  -  0.0001  C0-79  (  5.2  ) 

where  wsm  is  in  cm/sec,  and  C  is  in  mg/1.  With  reference  to  eq.2.t, 
eq.5.2  indicates  k-|  -  0.0001,  and  n^  -  0.79.  In  the  range  of  C  = 
1,500^10,027  mg/1,  wsm  remained  unchanged,  equal  to  0.033  cm/sec. 

Longitudinal  variation  of  suspended  sediment  settling  velocity. 
Flg3.5.29(a)  and  (b)  show  the  distance  distribution  of  dimensionless 
settling  velocity  of  suspended  sediment,  ws/wsi .  Because  of  the 
complexity  involved  in  directly  measuring  the  settling  velocity  of 
settling  floes  (for  cohesive  sediment),  an  indirect  method  was  adopted 
to  compute  the  settling  velocity  using  the  mass  balance  equation, 
i.e.  eq.2.6.  For  all  WES  and  C0EL  experimental  tests,  the  ratio  tij/t0<j 
was  on  the  order  of  0.01  due  to  the  small  velocity  in  the  closed"end 
channel.  Therefore,  the  probability  of  a  particle  sticking  to  the  bed, 
1r,T>/Tcd>  in  l-ns  "lass  balance  equation  was  assumed  to  be  unity.  Also, 
the  deposition  rate  per  unit  bed  area,  dm/dt,  can  be  computed  using  the 
deposited  sass  divided  by  the  effective  duration  (I.e.  the  period  fr»m 
the  passage  of  the  front  at  the  location  of  Interest  until  the  end  of 
the  test),  and  collection  area.   In  addition,  the  spatially  interpo" 
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Figure  5.29. (a)  Dimensionless  Settling  Velocity  vs.  Dlmensionless 
Distance  From  Entrance  for  Cohesive  Sediments. 
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Figure  5.29. (b)  Dimensionless  Settling  Velocity  vs.  Dimensionless 
Distance  From  Entrance  for  Cones ionless  Sediments. 
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lated  mean  concentration  below  interface,  Cbl  was  used  for  C  in  the 
mass  balance  equation.  The  settling  velocity  of  suspended  sediment  was 
calculated  under  these  assumptions. 

Due  to  the  distinct  characteristics  of  cohesionless  and  cones Ive 
sediments,  examination  of  the  distance  variation  of  dlmenslonless 
settling  velocity  for  cohesive  sediment,  which  included  kaolinite  and 
Vicksburg  loess,  may  be  considered  separately  from  cohesionless 
sediment,  which  included  flyash  I,  II,  and  III.  As  seen  in  Figs. "5. 29 
(a)  and  (b),  in  general,  w3/wsi  decreased  with  increasing  x/H.  In 
addition,  the  slope  of  the  curve  may  be  related  to  w^/u^,  which  was 
found  to  be  a  significant  parameter  for  all  of  the  depositlon«related 
characteristics.  From  Figs. 5. 29(a)  and  (b),  it  Is  observed  that  for 
both  cohesive  and  cohesionless  sediments,  the  slope  of  the  lines  became 
steeper  with  increasing  wa1/u^. 

Longitudinal  variation  of  dispersed  particle  size  of  deposit. 
Fig. 5.30  shows  the  distance  variation  of  dlmenslonless  median  dispersed 
particle  size  of  the  deposit,  Am/<im-\  ,  where  dm-|  is  the  median  dispersed 
particle  size  of  the  deposit  at  the  entrance.  As  mentioned  in  section 
5.1.3,  the  dispersed  particle  size  represents  the  primary  particle 
size  from  which  floes  are  formed.  There  is  a  trend  for  dm/dmi  to 
decrease  with  distance  from  the  entrance. 

The  dependence  of  the  rate  of  the  decrease  of  dn/dm  with  distance 
on  w3i/u£  may  be  examined  using  the  results  of  COEL  test.  According  to 
the  plot,  the  rate  of  decrease  of  dm/dra.  with  distance  is  approximately 
the  same  for  all  of  four  kaolinite  tests,  but  is  different  from  two 
other  tests  with  flyash  III.  In  order  to  explain  why  dm/dffl  decreased 
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Figure  5.30.  Diraensloniess  Median  Dispersed  Particle  Size  of 
Deposit  vs.  Diraensloniess  Distance  From  Entrance. 
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with  distance  at  the  same  rate  for  ail  four  kaollnite  tests,  one 
significant  feature  of  cohesive  sediment  roust  be  noted;  In  the  case  of 
a  cohe3ive  sediment  such  as  kaollnite,  the  settling  units  or  floes  are 
composed  of  primary  particles.  Oftentimes,  faster  settling  floes  are 
composed  of  finer  particles  (Dixit  et  al.,  1982).  On  the  other  hand 
weakly  aggregated  floes  composed  of  larger  primary  particles  tend  to 
have  lower  settling  velocities  and  travel  further  than  floes  composed 
of  finer  particles.  There  is  therefore  no  unique  correlation  between 
the  settling  velocity  of  a  floe  and  the  primary  constitutive  particles. 

For  cohesionless  sediment,  such  as  flyash  III,  no  significant 
particle  aggregation  is  involved.  Accordingly,  the  trend  of  the  slope 
of  the  curves  in  Fig. 5. 30  becoming  steeper  with  increasing  Wg-j/u^ 
could  be  expected. 

Variations  of  settling  velocity  and  dispersed  particle  size  with 
main  channel  concentration.  Fig. 5. 31  3hows  the  variations  of  sediment 
settling  velocity,  w3,  and  median  dispersed  particle  size  of  deposit, 
dm  (mean  values  for  all  deposits  located  in  the  region  x  »  0. 0*3-0  m   in 
the  closed^end  channel),  with  the  corresponding  mean  main  channel 
concentration,  C0.  The  purpose  of  this  figure  is  to  examine  the 
influence  of  main  channel  concentration  on  settling  velocity  and 
dispersed  particle  size  at  the  area  near  the  entrance.   It  is  believed 
that  most  of  the  deposited  floc3  (for  cohesive  sediment)  in  the 
entrance  area  were  formed  in  the  main  channel  under  turbulent  conditi- 
ons, and  deposited  immediately  when  they  moved  into  the  closed^end 
channel,  where  the  flow  velocity  was  much  lower. 


15"4 


MEDIAN    DISPERSED  PARTICLE  SIZE  OF  DEPOSIT.d^m) 

lO  Q 


U"> 


O 

C\J 


1                  1 
(V  ro  <q-  up 

1              1 

JJi^-J 

o 

UJ  UJ  W  UJ 

o  o  oo 

"<n 

uouo 

1 

•  «  •    ♦ 

o    <    a    0 

- 

/ 

/ 
/ 

— 

/ 
/ 
/ 
/ 

y 

c** 

f  X 

^^^                in 

- 

J.-' 

/* 

•    ' 

/ 

/ 

!                           1 

1                           1 

lO 

to 


■q  | 

_r  ■ 

UJ    E 


< 

5 


z 
o 

5 


z 

a 

z 
o 
o 

z 

< 

UJ 


in 

9 

CO 

<\J 

o 

O 

O 

o 

CJ 

o 

a 

o 

o 
d 


o 
cd 
o 
o 


pas/wo)   sM'xil30"I3A   9NI1113S   1N3WKJ3S 


155 

It  can  be  seen  from  Fig. 5- 31  that  ws  increased  rapid.lv  with 
increasing  C0  as  C0  >  I ,180  mg/1.  while  dffi  iecreased  with  Increasing 
C0.  Since  one  would  expect  that  greater  aggregation  of  cohesive 
sediment  would  occur  under  conditions  of  higher  concentration  and 
turbulent  flow.  It  is  not  surprising  that  Fig. 5. 31  shows  that  the  floes 
formed  under  higher  suspension  concentration  possessed  higher  settling 
velocities  and  were  composed  of  smaller  primary  particles,  as  was  also 
noted  previously. 

Critical  deposition  shear  stress.  Table  5.3  gives  values  of  the 

critical  deposition  shear  stress,  tc(j,  of  kaolinite.   tc(j  is  defined  as 

the  shear  stress  below  which  all  initially  suspended  sediment  eventual" 

ly  deposits.  The  purpose  of  computing  tCC]  of  kaolinite  in  the  present 

study  was  to  explore  the  deaositional  property  of  kaolinite,  and  to 

provide  an  input  value  required  for  the  application  of  the  mathematical 

model.  The  determination  of  tgcj  was  based  on  the  mass  balance  equation 

with  measured  deposition  rate,  dm/dt,  mean  concentration,  C0,  settling 

velocity,  w3,  and  mean  flow  velocity,  u0,  in  the  main  channel.  Here, 

the  settling  velocities  of  suspended  sediments  in  main  channel  are 

assumed  to  be  the  same  as  the  ones  obtained  by  settling  column  analy« 

sis,  even  though  the  flow  conditions  were  quite  different  in  the  two 

cases.  The  bed  shear  stress,  iD,  was  calculated  using  the  formula,  i^ 

_2 
»  f0puu0/8,  where  the  friction  factor,  f0,  was  determined  using  a  Moody 

diagram  (Dally  S  Karleman.  1S66J  based  on  the  Reynolds  number,  4Hu0/v, 

where  R  13  the  hydraulic  radius.  The  mean  tc(j  of  kaolinite  from  test; 

COEL-3,  5,  9,  13.  and  1 1  is  equal  to  0.13  N/m2,  which  compares  well 

with  0.15  N/m2  obtained  by  Mehta  (1973). 


I9€ 

Table  5.3  Critical  Dfposltion  Shear  Stress  of  Kaollnlte 


Test  Number 

(COEL)  3        5        9       13       I* 


dm/dt 
(g/em2->sec)x 

105 

2.8 

1.1 

1.9 

1.1 

1.7 

Co 
(mg/1) 

2,218 

908 

1,926 

1,133 

1,692 

wsm 

(cm/sec) 

0.032 

0.022 

0.032 

0.031 

0.032 

(N/m2) 

0.019 

0.072 

0.075 

0.110 

0.097 

Tb/l:cd 

0.61 

0.15 

0.69 

0.68 

0.69 

Tcd0 
(N/m2) 

0.08 

0.16 

0.11 

0.16 

0.11 

Distance  variation  of  flocculatlon  factor.  Fig. 5. 32  shows  the 
longitudinal  distribution  of  the  flocculatlon  factor,  F,  defined  as  the 
ratio  of  the  median  settling  velocity  of  the  floes,  w3,  to  the  settling 
velocity  of  the  primary  particle  sizes,  w3d,  from  which  the  floes  are 
formed.  The  data  represented  in  the  figure  include  all  steady  COEL 
tests  for  kaolinite  except  C0EL-11,  in  which  deposition-*related  data 
were  not  obtained.  It  can  be  seen  from  this  figure  that  three  distinct 
trends  of  flocculation  factor  varying  with  distance  appear  and  correal 
pond  to  different  suspension  concentrations  in  the  main  channel.  The 
reason  for  Increasing  flocculation  factor  (in  COEL'2  &  3)  between  x  - 
0.0*1.5  m  is  that  most  suspended  floes  from  the  main  channel  were 
deposited  at  x  «  1.5*2.0  m  by  the  convectlve  transport  of  the  shear' 
induced  gyre  at  the  entrance.   It  was  previously  mentioned  that 
deposits  of  cohesive  sediment  can  consist  of  two  parts:  individual  or 
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weakly  aggregated  particles  and  relatively  large  floes  formed  by  finer 
primary  particles.  The  peak  of  the  COELr-2  curve  indicates  that  the 
particles  with  a  relatively  high  degree  of  aggregation,  formed  In  main 
channel  under  high  suspension  concentration,  dominated  the  deposit  at 
x  -  1.5"2.0  m.  Beyond  this  location,  F  could  be  expected  to  decrease 
with  x  because  most  large  floes  had  already  settled  out.  The  relative* 
ly  constant  F  values  between  x  -  1.5«5.2  m  observed  in  the  curve  of 
COEL1-3  suggest  that  at  moderate  main  channel  concentrations,  approxima-* 
tely  the  same  amounts  of  individual  particles  and  floes  had  deposited 
in  the  closed^end  channel.  At  low  concentrations  in  the  main  channel 
(COEL*t  &  5),  small  F  values  at  x  -  1.5  m  indicate  that  the  larger 
(weakly  aggregated)  individual  particles  dominated  the  deposit,  and 
thereafter,  F  increased  gradually  with  increasing  x  due  to  the  settling 
of  the  floes  with  a  lower  degree  of  aggregation,  which  were  formed  in 
the  main  channel  under  low  suspension  concentration. 

Flocculatlon  factor  variation  with  dispersed  particle  size. 
Fig. 5. 33  shows  the  flocculatlon  factor  variation  with  median  dispersed 
particle  size,  which  was  obtained  from  both  the  WES  and  COEL  tests.  It 
can  be  seen  from  the  figure  that,  in  general,  F  increased  with  decreafc 
sing  median  dispersed  particle  size,  dm.  Note  that  for  kaolinite 
tests,  higher  suspension  concentration  in  the  main  channel,  C0,  seems 
to  have  generated  floes  with  larger  flocculatlon  factor  formed  by  the 
same  dm.  The  flocculatlon  factors  of  flyash  I,  flyash  II,  and  flyash 
III  are  observed  to  be  very  low,  indicating  very  weak  cohesion.  Also, 
the  low  flocculation  factor  for  Vicksburg  loess  was  found  under  low 
suspension  concentration.   There  13  a  reasonable  agreement  with  the 
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Figure  5.33.  Flocoulatlon  Factor  vs.  Median  Dispersed  Particle  Size. 
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results  obtained  by  Mlgniot  ("'968)  and  Dixit  (1982),  notwithstanding 
the  3omewhat  different  test,  conditions  among  the  three  studies. 

Longitudinal  variation  of  deposition  rate.  According  to  dimen* 
sional  analysis  (eq.3.19),  the  normalized  deposition  rate,  6/<5-|,  was 
found  to  be  dependent  upon  the  dlmensionless  distance,  x/H,  the 
widttWto^depth  ratio,  B/H,  the  relative  settling  velocity,  Wg/u^,  and 
the  density  difference  at  the  entrance,  Ap^/pw.  Due  to  Insufficient 
data  for  examining  the  effects  of  all  of  these  dlmensionless  parameters 
on  6/S-\,   only  the  effect  of  the  major  Influencing  parameter,  w3-|/U£,  on 
6/ 61  is  examined  in  this  section  by  assuming  that  the  dependence  of 
S/«1  on  B/H  and  on  Ap-|/pw  are  relatively  small  under  the  specific  test 
conditions. 

Figs.S.S'Ua)  and  (b)  show  the  longitudinal  variation  «/«i  for 
cohesive  and  cohesionless  sediments,  respectively.  In  each  test,  it  is 
observed  that  6/61  decreased  with  increasing  x/H  from  the  entrance. 
Once  again,  the  relative  settling  velocity  parameter,  w^/u^,  should  be 
important  for  examining  the  slopes  of  the  curves  (I.e.  rate  of  decrease 
of  dlmensionless  deposition  rate  in  the  x  direction).  Due  to  the 
deposited  samples  being  collected  at  different  locations  near  the 
entrance  for  each  test,  it  Is  somewhat  difficult  to  examine  the  slopes 
of  the  curves  on  a  comparative  basis,  since  these  start  from  different 
locations.  Yet,  in  general,  the  trend  of  the  slope  of  the  curves 
becoming  steeper  with  increasing  Wg-j/u^  can  be  easily  seen  from  both 
figures.  Note  that  the  slope  trends  of  the  curves  for  S/6i  are  similar 
to  those  of  curves  shown  in  Figs. 5. 29(a)  and  (b)  for  the  dimensionless 
settling  velocity. 
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Figure  S.S'l.Ca)  Dimensionless  Deposition  Rate  vs.  Dimensionless 
Distance  From  Entrance  for  Cohesive  Sediments. 
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Relationship  between  sediment  flux  and  suspension  concentration  at 
entrance.  Fig. 5. 35  presents  the  results  of  an  attempt  to  provide  a 
relationship  of  some  practical  value  for  predicting  sedimentation  in 
closed«end  channels  of  the  type  examined  in  the  present  study  under 
similar  conditions.  The  data  points  shown  in  the  figure  are  obtained 
from  all  st  eady  tests  at  WES  and  COEL  except  C0EL*11,  in  which  no 
deposited  sediment  samples  were  obtained.  The  mean  sediment  flux,  S, 
was  computed  by  dividing  the  total  deposit  mass  by  test  duration  and 
the  area  of  the  lower  half  cross  section  at  entrance.  The  analytical 
derivation  of  eq.3.61  provides  a  way  to  predict  the  sediment  flux 
through  the  lower  layer  into  the  channel,  given  the  boundary  concentrate 
tion,  Ci ,  at  the  entrance,  water  depth  H,  and  specific  gravity  of 
sediment  G3.  Arithmetic  mean  values  of  H  and  Gs  were  computed  for 
12  tests  Included  in  Fig. 5. 35  and  are  equal  to  8.8  cm  and  2.55, 
respectively.  Also,  an  experimental  result  that  the  surface  outflow 
velocity  at  the  entrance  of  the  closedtJend  channel  at  steady  state  is 
about  30J  lower  than  the  initial  front  speed,  found  by  Gole  et  al. 
(1973),  is  used  to  obtain  a  value  (a  -  0.35)  in  eq.3.61.  Accord* 
ingly,  eq.3.61  can  be  reduced  to 

1  .5 

S  -  0.C15  C-,  (  5.3  ) 

where  the  sediment  flux,  S,  i3  in  g/m2«minute  and  C-|  is  the  depth* 
averaged  concentration  at  the  entrance  in  mg/1.  The  straight  line 
representing  eq.5.3  is  plotted  against  corresponding  data  points 
on  logrtlog  axes,  and  shows  a  good  agreement  between  prediction  and  the 
experimental  results. 
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Equation  3.6!  should  b*  applicable  not  only  for  prototype  closed* 
end  channels,  but  also  for  prototype  opan*ended  ones,  as  lung  as  tne 
flow  field  in  the  channel  is  quiescent  and  the  sediment  transport 
mechanism  is  dominated  by  a  turoidity  gradient.  Nonetheless,  for 
basins  with  short  lengths,  such  as  I  marina,  all  suspended  sediment 
transported  by  turbidity  current  may  not  deposit  before  it  reaches  the 
closed  end  of  the  basin.  As  a  consequence,  a  turbidity  current 
reflected  from  the  closed  end  will  occur.  Such  reflected  turbidity 
fronts  propagating  towards  the  entrance  were  observed  during  the  tests 
performed  at  WES.  For  that  particular  case,  the  use  of  eq.3.6l  to 
predict  the  sediment  flux  into  a  basin  needs  further  study. 

5.1  Results  of  Tidal  Effect  Tests 
In  coastal  regions,  tidal  niotion  is  an  important  factor  in  the 
transport  of  fluid  and  sediment.  The  study  of  tidal  effect  on  turbidi- 
ty current  and  associated  sediment  transport  noted  in  this  section  is 
preliminary,  as  it  is  based  on  only  two  tests  (COEL-9  &  10)  and  one 
reference  test  (C0EL*11),  i.e.  test  without  imposed  water  surface 
variation  in  the  main  channel.  The  effects  of  water  surface  variation 
on  front  speed  of  turbidity  current,  phases  of  front  propagation, 
and  mean  sediment  concentration  below  the  interface  are  described.  An 
example  of  the  contribution  of  tldefinuuced  deposition  in  a  closed*end 
channel,  relative  to  contribution  from  turbidity  currenfinduoed 
deposition,  is  given 
5.1.1  Longitudinal  Variation  of  Front  Speed 

Figure  5.36  shows  the  distance  variation  of  the  dimensionless 
front  speed  in  the  two  tidal  effect  tests  and  in  the  reference  test. 
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Arrows  shown  along  the  horizontal  axis  indicate  the  direction  of  water 
surface  variation  Imposed  In  the  headbay  (upward  and  downward  arrows  of 
dashed  line  indicate,  respectively,  rise  and  fall  of  water  surface  in 
test  G0Eu"9;  while  the  corresponding  arrows  of  the  solid  line  indicate, 
respectively,  rise  and  fall  of  water  surface  in  COEL*10),  and  the 
approximate  front  position  at  the  time  of  water  surface  change  in  the 
headbay.  The  Imposed  variation  of  water  surface  with  the  elapsed  time 
during  the  tests  can  be  seen  in  Flg.1.11. 

From  the  curves  of  C0EL»9  and  10,  the  response  on  front  speed  to 
the  first  pair  of  changes  in  water  level  can  be  seen.  It  appears  that 
at  the  distance  of  x/H  ■  85«95,  the  effect  of  an  increase  of  water 
level  in  C0EL*9  was  a  lessening  of  the  slope  of  the  curve;  however,  no 
effect  on  the  curve  slope  of  C0EL*10  by  decreasing  water  level  is 
observed.  In  the  range  of  x/H  -  115*125,  no  effect  on  the  curve  slope 
of  C0ELH9  by  increasing  water  level  is  observed;  however  the  effect  of 
a  decrease  of  water  level  in  COELMO  was  a  steepening  of  the  curve 
slope.  A  comparison  between  the  curves  of  C0ELS9  and  11  reveals  that 
an  increase  of  water  level  in  the  main  channel  seems  to  have  prevented 
the  front  speed  from  decreasing.  However,  due  to  a  lower  main  channel 
concentration  in  COEL*10,  its  dimensionless  front  speed  was  higher  than 
that  of  reference  te3t  C0EL*11  before  any  change  of  water  level  was 
imposed.  Therefore,  it  is  probably  unjustified  to  make  any  comparison 
of  the  effect  of  Imposed  water  surface  variation  for  this  case,  with 
the  reference  test. 
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5. 4. 2  Front  t  Ml  tier.   Var'J  a  t 1  on  wj  th  Elapggj  Time 

Aa  shown,  in  Fig. 5. 26,  three  characteristic  phases  (i.e.  initial 
adjustment  phase,  inertial  self^similar  phase,  and  viscous  3elf-*similar 
phase)  were  found  in  the  front  propagation  in  the  steady  COEL  tests. 
The  effect  of  water  surface  variation  on  the  phase  change  during  front 
propagation  may  be  examined.  Fig. 5. 57  shows  the  variation  of  dimen* 
slonless  fru.it  position,  Xf/H,  with  dimensionless  elapsed  time,  u^t/H, 
for  two  tidal  effect  tests  and  for  the  reference  test.  Again,  the 
arrows  shown  at  the  horizontal  axis  Indicate  the  direction  of  water 
surface  variation  imposed  in  the  headbay  and  the  corresponding  dimens 
sionless  eiapsed  time  at  the  time  of  change  of  water  level  in  the 
headbay. 

It  can  be  seen  from  Fig. 5. 37  that  the  front  motion  in  all  three 
tests  started  from  the  initial  adjustment  phase  (slope  -1),  and 
entered  the  inertial  self*similar  phase  (slope  ■  2/3)  at  approximately 
the  same  time  of  u^t/'H  »  57.  Subsequently,  the  curve  slope  of  refer* 
ence  test  COELM  1  started  decreasing  from  the  slope  of  2/3  at  u^t/H  - 
400,  and  eventually  became  1/5  (the  slope  for  viscous  selfsimilar 
phase)  at  u4t/H  »  850.  However,  for  both  the  tidal  effect  tests,  in 
which  the  water  surface  variation  imposed  in  the  headbay  started  at 
about  t.Al/H  *  330,  the  curve  slope  of  CO  ELMO  ("ebb  tide")  did  not 
decrease  from  the  Inertial  seif*sirailar  phase  until  u4t/K  ■  580. 
However,  the  front  in  C0EL«9  ("flood  tide")  still  propagated  within  the 
inertial  self  •-similar  phase  when  it  reached  the  closed  end  (at  u*t/K  » 
800).  Based  on  these  observations,  a  preliminary  conclusion  can  be 
drawn  that  water  surface  variation  in  main  channel,  regardless  of 
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increasing  or  decreasing  water  level,  will  extend  the  duration  of 
front  propagation  within  the  inertlal  aelfrisimllar  phase.  The  reason 
may  be  that  the  oscillation  in  the  flow,  generated  by  water  surface 
change,  provides  extra  inertial  force  to  the  propagating  front. 
5.1.3  Longitudinal  Variation  of  Mean  Concentration  below  Interface 

Fig. 5.38  shows  the  distance  variation  of  dimen3ionless  mean 
concentration  below  interface  for  the  two  tidal  effect  tests  and  the 
reference  test.  It  can  be  seen  that  the  imposed  water  surface  varla* 
tions  in  both  tests  had  no  significant  effect  on  the  longitudinal 
distribution  of  mean  concentration  below  the  Interface  at  steady  state. 
Possibly,  there  were  also  no  significant  effects  on  the  amount  of 
sediment  deposited  in  the  channel. 

Note  that  the  changes  of  water  surface  in  C0ELr9  and  10  are  not  a 
realistic  simulation  of  astronomical  tide.  Astronomical  tide  simulaR 
tion  requires  a  periodic  variation  of  water  surface,  with  tidal 
amplitude  and  period  scaled  appropriately  to  represent  similarity  with 
prototype  tide.  To  evaluate  the  tidal  effect  on  the  mechanics  of 
front  propagation  and  associated  sediment  deposition,  further  investig* 
ation  is  required. 

5.1.1  Tlde-ilnduced  to  Turbidity  Current*Induced  Deposition  Ratio 
In  section  3.3.5,  the  relationship  between  tideflinduced  and 
turbidity  eurrenttinduced  deposition  in  closed«end  channels  has  been 
considered  (eq.3.67).  According  to  this  relationship,  the  ratio  of 
tide-'induced  to  turbidity  current-Hnduced  deposition  ma3S,  MT/MD,  can 
be  obtained  if  the  excess  density  at  the  channel  entrance,  Api/pw, 
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17? 
relative  tidal  amplitude,  a0/K,  and  channel  length  to  tidal  wave  length 
ratio,  !.,/L' ,  are  known. 

In  order  to  examine  the  deposition  contributed  by  tidal  motions  as 
percentage  of  that  by  turbidity  current  in  prototype  scale  channels  and 
tidal  conditions,  one  typical  example  is  sufficient.  Consider  a  I ,000 
m  long  and  5   m  deep  olosed*end  channel,  under  the  influence  of  tide 
having  a  0.5   m  amplitude  and  a  12  hours  period,  with  a  relative  density 
gradient  at  the  entrance,  Ap-|/pw  ■  0.001.  The  local  tidal  wave  length 
can  be  computed  as  L'  =  CgH)1/2T  =  302,100  m.  Hence  the  tlde*induced 
to  turbidity  current>*induced  deposition  mass  ratio  is 


MT  Api  --1/2  a0   L 

—  -  11.1(2  ( )    (— )  (— )  =  0.12    C  M  ) 

MD  pw      H    L' 


The  result  indicates  the  deposition  mass  induced  by  tidal  motions  i3 
approximately  12?  of  that  induced  by  turbidity  currents.  It  can  be 
seen  from  the  formula  that  the  contribution  of  tide*induced  deposition 
becomes  important  when  the  area  of  water  suface  in  the  channel  and 
tidal  range  are  large,  or,  the  channel  is  3hallow  and  the  density- 
gradient  is  small. 


CHAPTER  6 
CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  Summary  of  Investigation 

In  order  to  explore  the  distinction  between  low  Reynolds  number 
front  behavior  of  a  turbidity  current  and  a  non^settling  gravity 
current  (e.g.  salinity  or  teraperature»induced  density  current),  two 
sets  of  laboratory  experiments  were  conducted  in  specially  designed 
flume  systems.  Characteristic  results  are  described.  These  results 
have  been  discussed  on  the  basis  of  significant  dimensionless  paramef- 
ters  obtained  from  dimensional  analyses.  Discussed  physical  aspects  of 
front  behavior  include  front  nose  height  to  head  height  ratio,  front 
head  height  to  neck  height  ratio,  dimensionless  initial  front  speed, 
dimensionless  instantaneous  front  speed,  local  densimetrlc  Froude 
number,  dimensionless  front  position,  dimensionless  excess  front 
density,  and,  most  important,  relative  settling  velocity. 

Typical  results  of  the  flow  regime  in  the  channel,  obtained  by 
simulation  using  a  specially  developed  and  calibrated  numerical  model, 
include  time  variation  of  horizontal  velocity,  and  time  and  distance 
variations  of  vertical  velocity  profile.  An  attempt  at  deriving 
analytical  solutions  for  the  horizontal  velocity  profile  at  steady 
state  was  made.   Also,  tyoical  results  related  to  suspension  concentrat- 
ion obtained  from  both  experimental  measurement  and  numerical  model 
simulation  are  given.  These  Include  time  variation  of  suspension 


173 
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concentration,  and  distance  variation  of  the  vertical  concentration 
profile.  The  longitudinal  distribution  of  mean  concentration  below  the 
interface  at  steady  state  was  examined  using  experimental,  numerical, 
and  analytical  results, 

Deposltionfrelated  characteristic  results,  which  were  mainly 
obtained  from  experimental  measurements,  include  the  longitudinal 
distributions  of  settling  velocity,  median  dispersed  particle  diameter, 
flocculatlon  factor,  and  deposition  rate.  Formulations  for  estimating 
turbidity  driven  sediment  flux  into  the  channel,  and  the  ratio  of 
tldeMnduced  to  turbidity  currentninduced  deposition  have  been  dis* 
cussed. 

Finally,  the  effect  of  water  surface  variations  in  the  main 
channel  on  Instantaneous  front  speed,  duration  of  each  characteristic 
phase  in  front  propagation,  and  mean  concentration  below  the  Interface 
were  examined. 

6.2  Conclusions 

The  groupings  of  diraensionless  parameters  derived  from  dimensional 
analysis  for  front  nose  height,  front  head  height,  initial  front  speed. 
Instantaneous  front  speed,  front  position,  and  deposition  rate  in 
section  3.1  have  provided  an  appropriate  framework  for  discussing  and 
presenting  experimental  results.  Also,  some  results  of  dimensionless 
analysis  have  been  confirmed  by  comparing  present  experimental  results 
with  those  obtained  by  other  investigators. 

In  accordance  with  calibration  and  verification  results  of  the 
developed  numerical  model,  this  two-dimensional,  explicit,  f lnite*dlffrf 
erence  model  is  shown  to  have  a  good  capability  of  simulating  flow 
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regime,  concentration  distribution,  and  sediment  deposition  in  a 
closedflend  channel.  In  the  model,  the  effect  of  cohesive  sediment 
aggregation,  buoyancy  effect  due  to  fluid  stratification,  and  tidal 
effect  have  been  taken  into  account. 

Experimental  facilities,  specially  designed  for  keeping  flow 
velocity  and  suspension  concentration  In  the  main  channel  unchanged 
throughout  a  test,  effectively  established  a  stable  flow  regime  and 
kept  concentration  constant  for  finer  sediments  (e.g.  kaolinlte).  For 
sedlment3  having  high  settling  velocities  (e.g.  flyash  I),  mixers  in 
the  slurry  mixing  tanks  were  not  capable  of  preventing  the  material 
from  settling  in  the  tank.  It  should  be  noted  that  turbulent  flows  in 
the  main  channel  affect  the  behavior  of  turbidity  current  in  the 
closedftend  channel  by  changing  the  vertical  concentration  profile  at 
the  main  channel,  the  degree  of  aggregation  of  suspended  particles,  and 
the  length  of  mixing  gyre  at  the  entrance  of  the  closed^end  channel. 

Due  to  the  difficulty  of  measuring  the  vertical  velocity  profiles 
in  the  slow  moving  and  turbid  waters,  the  temporal  and  spatial  variaf* 
tions  of  horizontal  velocity  In  the  channel  were  examined  by  using 
numerical  simulation  results.  It  was  observed  that  the  horizontal 
velocity  at  2.5  cm  elevation  at  a  fixed  location  increased  from  zero  to 
the  maximum,  then  decreased  and  reached  a  steady  value  as  time  proceeds 
ed.  The  maximum  velocity  occurred  at  the  moment  when  the  front  passed 
the  point  of  interest.  A  steady  velocity  was  attained  even  though  the 
front  continued  to  propagate  farther  downstream  in  the  channel.  It  was 
also  found  that  the  Interface  elevation  of  the  velocity  profile  at  a 
fixed  location  Increased  with  time,  and  that  the  velocity  magnitude  and 
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interface  elevation  decreased  with  distance  from  the  channel  entrance. 
At  steady  state,  the  magnitudes  of  flow  discharge  In  the  layers  aboye 
and  below  the  interface  were  equal,  satisfying  volumetric  conservation 
of  the  fluid. 

The  suspension  concentration  at  a  fixed  position  increased  rapidly 
from  zero  to  a  certain  value  (when  the  front  was  passing  the  location 
of  interest),  then  continued  increasing  and  eventually  reached  a  steady 
value.  It  appears  that  the  propagating  front  has  little  effect  on 
concentration  distribution  upstream.  Measured  vertical  concentration 
profiles  revealed  that  water  in  the  layer  above  the  interface  is  not 
completely  free  of  sediment.  This  suggests  that  a  part  of  the  finer 
particles  has  been  transported  upward  by  diffusion  and  convection. 
However,  the  characteristics  of  the  vertical  concentration  profile, 
including  relatively  small  concentrations  in  the  upper  layer,  the 
elevation  of  maximum  concentration  gradient  located  at  approximately 
zero"velocity  elevation  and  relatively  uniform  concentration  in  the 
lower  layer  all  suggest  that  a  two-layered  model  with  clear  upper  layer 
Is  adequate  for  theoretical  development.  The  shape  of  the  vertical 
concentration  profile  was  similar  at  various  locations  in  the  channel, 
but  the  magnitude  decreased  with  distance  from  the  entrance. 

In  view  of  the  good  agreement  for  concentration  distributions 
between  experimental  results  and  numerical  simulation,  the  settling 
velocity  computed  from  the  mass  balance  equation,  and  used  as  input 
data  in  the  numerical  simulation,  can  be  considered  to  appropriately 
represent  the  settling  velocity  of  suspended  particles  in  the  channel 
under  specific  flow  conditions. 
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Small  propagating  surface  waves,  which  result  from  the  elevateo 
water  surface  dur  ing  the  intrusion  of  dense  fluid,  wee  observed  in 
numerical  simulation  results.  When  a  stationary  sediment  wedge  is 
established,  the  water  surface  rises  from  the  entrance  towards  the 
wedge  toe,  then  continues  horizontally  to  the  end  of  the  channel. 
Momentum  conservation  principle  over  the  channel  reach  from  the 
entrance  to  the  wedge  toe  was  demonstrated  by  the  good  agreement  of 
AHmax  (the  maximum  surface  rise),  obtained  by  an  analytic  approach  and 
by  numerical  simulation. 

Front  nose  height  to  head  height  ratio,  h-|/h2,  was  found  to  be 
strongly  influenced  by  bottom  friction,  and  was  found  to  decrease  with 
increasing  local  head  Reynolds  number,  Ren  =  Ufhj/v.  The  consistent 
trend  of  h-j/hj  vs.  Ren  obtained  from  various  investigations  suggests 
that  the  relationship  between  hi/hg  and  Ren  shown  in  Fig. 5. 21  is  valid 
for  all  types  of  gravity  currents. 

An  approximately  constant  value  of  the  ratio  of  front  head  height 
to  neck  height,  hj/hj,  over  the  ranges  of  local  neck  Reynolds  number, 
Ren  =  urh3/v  -  20' 1 ,000,  and  fraction  depth,  h3/H  ■  0. 19*0. 39  was 
observed.  The  arithmetic  mean  of  all  data  points  for  this  ratio  is 
equal  tu  1.3.  The  difference  between  hj/hj  obtained  in  the  present 
study  and  that  of  Keulegan  (1958),  who  found  h2/h3  =  2.0,  resulted  from 
different  ranges  of  hj/H  in  the  two  studies.  This  observation  can  be 
supported  by  the  finding  of  Simpson  and  Brltter  (1979)  that  hp/t^ 
should  decrease  with  increasing  the  fraction  depth,  h3/H.   [f  should 
furthermore  be  noted  tnat  h2/h3  was  found  to  be  nearly  Independent  of 
the  relative  settling  velocity,  ws/uf. 
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Dimensionless  Initial  front  speeu,  Ufi/u^,  was  found  t.c  be 
independent  of  the  den.olmetric  Reynolds  number,  Re  ,  and  the  wldthi-to* 
depth  ratio,  B/H.  The  arithmetic  mean  of  Uf^/u^  was  0.1)3,  which  is 
close  to  0.16  obtained  In  Keulegan's  11957)  Investigation  of  lock 
exchange  flow  of  saline  water.  This  observation  Is  Indicative  of  the 
dynamic  similarity  between  turbidity  current  and  sallnityMriven 
current  during  the  Initial  stages  of  front  propagation. 

An  exponentlalRtype  decrease  of  the  dlmenslonless  front  speed 
Uf/Uf|  with  distance  observed  in  the  present  study  was  due  to  drag 
force  against  front  motion,  and  the  reduction  of  motive  force  by 
sediment  settling  and  interfaclal  mixing.   Also,  the  ratio  of  uf/Uf| 
was  found  to  be  strongly  dependent  upon  the  relative  settling  velocity, 
ws1/uA,  for  both  of  kaollnlte  and  flyash  III.  The  dimenslonless  front 
speed,  Uf /Uf 1 ,  decreased  faster  with  distance  as  the  relative  settling 
velocity,  Wgj/u^,  increased.  For  turbidity  current,  the  rate  of  the 
decrease  of  the  front  speed  with  distance  Is  faster  than  that  of 
salinity-induced  gravity  current,  in  which  the  front  density  decreases 
with  distance  only  by  Interfaclal  mixing. 

Local  densimetrlc  Froude  number,  FrA,  was  found  to  increase  with 
Increasing  local  neck  Reynolds  number,  Ren,  in  the  range  of  Ren  » 
25-320.  The  values  of  FrA  obtained  in  the  present  study  were  slightly 
lower  than  what  Keuiegan  0958)  found,  which  might  result  from  the 
effect  of  concentration  stratification  at  the  front  and  the  use  of 
bottom  concentration  for  representing  the  mean  concentration  at  the 
front. 
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Three  distinct  phases  (initial  adjustment  phase,  inertial  self"sl« 
milar  phase,  and  viscojs  selfsimilar  phase)  found  in  the  front 
propagation  of  sallnltyUnduced  gravity  current,  were  also  found  in 
the  front  propagation  of  turbidity  current.   Ail  three  characteristic 
curve  slopes  (slope  -  1,  ?/3,  and  1/5)  corresponding  to  tnree  phases 
were  observed  in  tests  with  kaollnite.  The  duration  of  front  propaga" 
tion  in  the  initial  adjustment  phase,  in  which  the  front  proceeds  at  a 
constant  speed,  lasted  longer  for  the  front  with  smaller  wsi/u^.  Yet, 
the  wsl/u^  had  little  effect  on  the  durations  of  propagation  in  the 
inertial  self»>similar  and  viscous  selfHsimilar  phases.  In  flyash  III 
tests,  only  two  characteristic  curve  slopes  (slope  ■  2/3,  and  l/5)f 
corresponding  to  the  last  two  phases,  were  observed  due  to  the  high 
relative  settling  velocity,  Wg^/u^.  In  addition,  the  duration  of  front 
propagation  in  the  inertial  selfslmllar  phase  (slope-2/3)  was  longer 
for  the  front  having  the  smaller  Wg^/u^. 

Keulegan  (1957)  has  pointed  out  that  the  front  dilution  of 
sallnltyMnduced  density  current  through  interfacial  mixing  is  faster 
in  lock  flows  at  relatively  small  Reynolds  numbers.  In  the  present 
study,  an  exponentlal^type  decrease  of  dimenslonless  excess  front 
density,  Ap'/Ap^,  with  distance  from  the  entrance  was  observed.   In 
this  case,  the  decrease  is  due  to  the  combined  effect  of  sediment 
settling  and  interfacial  mixing.  The  rate  of  the  decrease  was  found  to 
increase  with  increasing  relative  settling  velocity,  wsl/u^. 

An  exponentialrttype  decrease  with  distance  of  mean  concentration 
below  the  interface  at  steady  state  was  also  observed.   Accordingly, 
the  assumption  of  a  constant  mean  concentration  below  the  Interface 
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made  by  many  Investigators  for  the  purpose  of  theoretical  development 
is  inappropriate.  In  all  tests,  the  rate  of  the  decrease  of  dimension- 
less  mean  concentration  below  the  Interface,  C^/C^,  ,  with  distance 
became  greater  as  ws-|/ua  increased.  The  ability  of  the  numerical  model 
for  simulating  the  longitudinal  distribution  of  mean  concentration 
below  the  interface  was  good.  The  analytic  derivation,  eq.3.52,  for 
the  longitudinal  distribution  of  mean  concentration  below  the  interface 
was  found  to  be  only  good  for  relatively  small  w3-|/uA  values.  ■  Uso, 
the  discrepancy  between  analytic  and  experimental  results  became 
greater  with  increasing  values  of  w^/u^. 

The  variation  of  settling  velocity  with  concentration  for  fine" 
grained  sediments  used  was  investigated  by  using  settling  column 
tests.  It  was  found  that  the  median  settling  velocities,  w3m,  of 
flyash  I,  II,  III,  and  silica  flour  were  not  completely  Independent  of 
the  initial  column  concentration,  C,  but  increased  slightly  with 
increasing  concentration.  The  median  settling  velocity  of  Vicksburg 
loess  was  also  found  to  slowly  increase  with  increasing  C,  over  the 
range  of  C  =  t23"1,939  mg/1,  and  wam  of  Cedar  Key  mud  was  equal  to 
0.001  cm/sec  for  C  -  323C373  mg/1.  The  median  settling  velocity  of 
kaolinlte  was  found  to  be  concentration^dependent  over  the  range  of  C  ■ 
M3M«9M  mg/1,  and  could  be  expressed  as 


0.79 
=  0.0001  C  (  6.1  ) 


where  wsrn  is  in  cm/sec,  and  C  is  in  mg/1.   In  the  range  of  C  -  1,b00- 
10,027  mg/1,  wsm  of  kaolinlte  was  found  to  be  independent  of  C. 
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As  mentioned  above,  the  settling  velocity  of  the  deposit  compute" 
from  the  mass  balance  equation  ,»as  found  to  adequately  represent  the 
local  settling  velocity  of  suspended  sediment  In  tne  channel.  The 
settling  velocity  of  the  deposit  was  found  to  decrease  with  distance 
from  the  entrance  in  both  kaolinite  and  flyash  III  tests.  The  rate  of 
decrease  of  the  dimensionless  settling  velocity,  w3/w3i  ,  with  distance 
became  greater  as  the  relative  settling  velocity,  w3-|/u^,  increased. 

The  uniformity  of  particle  3ize  distribution  of  the  deposit  was 
examined  using  a  modified  sorting  coefficient,  (dgj/djj)'' 2.  This 
coefficient  was  found  to  be  approximately  the  same  everywhere  In  the 
channel  for  both  kaolinite  and  flyash  III  tests.  The  modified  sorting 
coefficient  in  all  steady  tests  using  kaolinite  ranged  from  1 .9  to  2.5, 
while  in  flyash  III  tests  it  ranged  from  1.3  to  1.6.  This  result 
suggests  that  the  particle  size  distribution  of  kaolinite  deposit  was 
more  uniform  than  that  of  flyash  III.  Again,  the  median  dispersed 
particle  size  of  the  deposit  was  found  to  decrease  with  distance  from 
the  entrance  because  of  sorting  effect.  Because  there  is  no  unique 
correlation  between  the  settling  velocity  of  a  floe  and  the  primary 
constitutive  particles  size,  as  mentioned  above,  the  dependence  of 
distance  variation  of  the  dimensionless  median  dispersed  particle  size, 
dm/dmi ,  on  w3i/u£  in  kaolinite  test.3  was  not  found  to  be  significant. 
In  flyash  III  tests,  the  rate  of  the  decrease  of  dm/dmi  with  distance 
became  greater  as  w3^/U£  Increased. 

Under  turbulent  flow  conditions  (Re  -  4Ru0/v  =  i 2, 965*27, 350)  in 
the  main  channel,  it  was  found  that  the  settling  velocity  of  kaolinite, 
w3,  increased  rapidiy  from  0.007''  cm/sec  at  suspension  concentration  C0 
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-  1.375  mg/1  to  ws  -  0.015  cm/sec  at  C0  -  3, 7!*6  mg'l.  On  the  other 
hand,  the  median  dispersed  particle  size,  dm,  from  which  the  floes  were 
formed,  lecreased  from  2.78  pm  at  C0  ■■  1,375  mg/1  to  1.39  um  at  co  " 
3,716  mg/1.  The  critical  deposition  shear  stress  of  kaolinlte,  t0(], 
which  is  defined  as  the  shear  stress  below  which  all  Initially  sus* 
pended  sediment,  eventually  deposits,  was  found  to  be  0.13  N/m2,  very 
close  to  0.15  N/m2  obtained  by  Mehta  (1973). 

A  trend  of  flocculation  factor,  F,  increasing  with  decreasing 
median  dispersed  particle  size,  dm,  was  observed,  similar  to  the 
findings  by  Migniot  (1968)  and  Dixit  et  al.  (1982).  Also,  that  primary 
particles  at  higher  concentrations  tend  to  form  relatively  large  floes 
was  observed.  Three  typical  trends  of  longitudinal  distribution  of  the 
flocculation  factor  with  respect  to  different  ranges  of  concentration 
were  identified.  The  flocculation  factor  Increased  from  entrance  to 
x  *  1.5  m  due  to  shear*induced  circulation,  and  at  x  =  1.5*2.0  m,  F  was 
large,  corresponding  to  a  high  main  channel  concentration.  At  loca^ 
tions  beyond  x  =  2.0  m,  F  rapidly  decreased  with  distance  in  the  test 
having  the  highest  main  channel  concentration,  C0  -  3,716  mg/1, 
remained  unchanged  with  distance  in  the  test  having  C0  =  2,218  mg/1, 
and  increased  slowly  with  distance  in  the  tests  having  C0  »  908*1  ,375 
mg/1. 

Typically,  there  was  significant  sediment  deposition  in  the  zone 
near  the  entrance,  followed  by  an  approximately  exponentialfttype 
decrease  of  deposition  farther  along  the  channel.  This  was  the  case 
because  most  particles  with  relatively  high  settling  velocities  kept 
suspended  in  the  main  channel  deposited  immediately  due  to  quiescent 
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water  in  the  channel.  The  rate  of  the  decrease  of  diraensionlcss 
deposition  rate,  S/6],   with  distance  consistently  increased  as  tne 
relative  settling  velocity,  ws-|/u4,  increased,  regardless  of  sediment 
type.  The  numerical  model  was  shown  to  have  good  capability  for 
simulating  the  deposition  rate  in  the  channel. 

The  relationship  derived  for  predicting  the  sediment  flux,  S,  into 
the  closed^end  channel  by  turbidity  currents  was 


gH      1     1/2  3/2 

S  =  a  [  —  (  1  »  —  )  ]   C:  (  6.2  ) 

Pw       Gs 


This  relationship  is  valid  not  only  for  the  closedaend  channels,  but 
also  for  the  open-ended  channels,  as  long  as  the  flow  condition  in  the 
channel  is  quiescent.  Here,  H  is  the  water  depth,  Gs  is  the  specific 
gravity  of  sediment,  and  C-|  is  the  depthimean  concentration  at  the 

entrance.  Note  that  the  sediment  flux  is  proportional  to  H1/2  and 

—3/2 

C]      ,  but  is  independent  of  sediment  settling  velocity. 

The  tide*induced  to  turbidity  current^induced  deposited  mass 

ratio  was  derived  as 


MT  Api  i1/2  a0   L 

—  -  11.42  ( )    (— )  (— )  (  6.3  ) 

MD  pw       H    L' 


where  Hj,  Mn  are  tlde-'lnduced  and  turbidity  currents-induced  deposited 
sediment  mass  over  a  tidal  cycle,  respectively,  a0  is  the  tidal 
amplitude,  and  L'  Is  the  tidal  wave  length  in  shallow  water.  The 
significance  of  tide^induced  deposition  becomes  obvious  when  the  area 
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of  water  surface  in  the  ehannei  la  large,  tidal  range  is  high,  the 
channel  la  shallow,  and  the  density  gradient  at  the  entrance  is  low. 

Based  or.  the  results  of  two  tidal  effect  tests,  a  preliminary 
conclusion  wa3  drawn  that  water  surface  variations  imposed  in  the  main 
channel,  regardless  of  whether  the  level  is  increasing  or  decreasing, 
will  extend  the  duration  of  front  propagation  within  the  inertlal 
selffcslmilar  phase,  in  which  the  front  proceeds  at  a  speed  proportional 
to  t*1/3.  Water  surface  variation  had  no  significant  effects  on  the 
longitudinal  distribution  of  mean  concentration  below  the  interface  at 
steady  state. 

6.3  Recommendations  for  Future  Study 
Although  the  two-dimensional,  explicit,  finite  diffence  numerical 
model  realistically  simulates  the  vertical  movement  of  the  free-3urface 
and  hence  the  propagation  of  surface  waves,  the  numerical  time  step  is 
severely  restricted  by  the  time  required  for  the  surface  wave  to 
propagate  a  distance  of  approximately  one  spatial  grid.  It  should 
however  be  noted  that  the  amplitude  of  the  surface  wave  generated  by 
the  bottom  density  current  is  much  smaller  than  the  total  water  depth, 
and  thus  is  not  physically  of  concern.  Using  a  rigid-lid  model  (Sheng, 
1986)  instead  of  a  free^surface  model  for  simulating  the  physics 
associated  with  density  currents  may  save  much  computer  time,  and  still 
achieve  the  same  results. 

la  numerical  simulation  there  must  be  temporally  and  spatially 
continuous  distributions  of  physical  quantities  to  keep  the  numerical 
scheme  stable;  therefore  the  suspension  concentration  at  the  sharp 
front  was  stretched  out  in  the  longitudinal  direction  during  Simula- 
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tion.   AcGcdinglj ,  this  mathematical  model  Is  not  applicable  to 
detailed  stud'.es  of  front  Shape  and  front  speed  of  gravity  currents. 
For  investigators  interested  in  the  front  behavior  of  such  currents,  a 
special  treatment  which  adopts  a  higher  order  difference  approximation 
of  a  Hermlte  interpolation  function  in  the  convective*'diffuslon 
formulation,  as  suggested  by  Fischer  (1976),  can  be  U3ed  to  Improve  the 
numerical  stability  and  reduce  numerical  errors  In  the  concentration 
prediction  at  the  sharp  front. 

Most  experiments  conducted  in  the  present  study  generated  low 
Reynolds  number  turbidity  currents  (ReA  -  321*3,830),  and  the  physics 
associated  with  such  currents  Is  Reynolds  numberfcdependent.  However, 
prototype  turbidity  currents  typically  have  much  higher  Reynolds 
numbers.  Thus,  more  realistic  prototype  phenomena  at  higher  Reynolds 
number  turbidity  currents,  in  which  the  physics  of  concern  is  Reynolds 
numberoindependent,  should  be  laboratory  tested. 

The  method  of  settling  column  analysis  is  a  good  alternative  to 
determine  the  settling  velocity  distribution  of  fine-grained  sedi'- 
ments.  A  problem  encountered  when  conducting  the  analysis  is  the 
significant  drop  of  water  surface  in  the  column  during  suspension 
sample  withdrawals.  A  minimum  of  30  ml  of  liquid  sample  is  needed  for 
determining  sample  concentration  by  the  gravimetric  procedure.  One 
possible  way  to  reduce  water  surface  reduction  In  the  settling  column, 
and  thus  improve  accuracy,  is  by  using  a  turbidimeter  for  which  only  a 
few  ml  of  sample  is  required  per  withdrawal. 

Three  physical  phases  were  identified  in  the  front  propagation 
of  gravity  currents  by  Huppert  (1982)  and  Rottman  and  Simpson  (1983). 
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In  each  phase,  front  propagatlor  is  dominated  by  specific  forces,  ana 
the  front  position  is  a   specific  function  of  the  elapsed  time.  But 
no  unique  functions  were  found  tc  describe  front  movement  from  the 
instant  of  gate  opening  to  when  it  reaches  the  closed  end.  In  future 
research,  if  one  can  examine  the  front  movement  of  gravity  currents  by 
looking  at  the  net  forces  acting  on  the  front,  i.e.  £Fj  -  mf(duf/dt), 
where  Fj  is  the  component  force  acting  on  the  front,  and  mf  is  the 
front  mass,  it  may  be  possible  to  find  a  unique  relation  between 
front  speed  and  elapsed  time. 


APPENDIX  A 

SPATIAL  GRID  SYSTEM,  FLOW  CHART, 

AND  STABILITY  CONDITION  FOR  NUMERICAL  MODELING 


A.I  Spatial  Grid  System  for  Numerical  Modeling 


(i-l)lti    Colunn    (i+l)th 
"i-l,N 


Fig.A.1.  Spatial  Grid  System  for  Numerical  Modeling 
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k.2   Flow  Cha-t  for  Numerical  Calculation 


Start 


/ 


Read  the 
|  Input  Data 


Set  Up 
Initial  Conditions 


Compute  Water  Surface 
elevation  for  the  Next  Time  Step 
(Equation  3.1)1  ) 


Compute  Density 
Gradient  and  Pressure  Gradient 

(Equations  3. 30  and  3. 31  ) 


Compute  Bed  Shear  Stress 
and  the  Probability  of  Deposition 
(Equation  3.3I  ) 
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Compute  Settling  Velocity  of 
Suspended  Sediments,  Richardson  Number, 
Momentum  and  Mass  Diffusion  Coefficients 
(B33ed  on  Equation  3-t5  and  Empirical  Formulas) 


Compute  Horizontal 

Velocity  for  the  Next  Time  Step 

(  Equation  3.t3  ) 


Compute  Sediment 

Concentration  for  the  Next  Time  Step 

(  Equation  3. It  ) 


Compute  Vertical 

Velocity  for  the  Next  Time  Step 

(  Equation  3.t2  ) 
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Print 
Water  Surface  Elevation, 
Horizontal  and  Vertical 
Velocities,  Sediment  Concentration 


<^   Stop  or  Not 


Reset  Variables 

and 
Advance  Time  Step 
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To 

A. 3  Derivation  of  Numerical 

Stability  Cond 

ition 

investigate 

the  linear   st 

.ability 

of  eq.3 

. ii7 ,   the  continuous 

functior 

1  f(x)   can  be  expressed  as  a  form 

of  complex  Fourier  series 

(Abbott, 

1979),    i.e. 

n 

f     - 
J 

n                 2itk 

I  F     exp(    i  J    Ax 

k     k                    X 

) 

(   A.I    ) 

where  * 

is  the  spatl 

al  interval  of  perio 

diclty. 

Using 

the  same   complex 

Fourier 

series  expression  for  subsequent 

continuous  functions 

n+1 
f 
J 

n+1 
=   I  F       exp( 
k     k 

2nk 
1  J 

X 

Ax    ) 

(   A. 2   ) 

n«1 
f 
J 

n«1 
-IF       exp( 
k     k 

2-nk 

1  J 

X 

Ax    ) 

(   A. 3   ) 

n*1 
f 
j  +  1 

n»1 

-If       exp[ 
k     k 

2irk 

i  ( 

X 

j+1)    Ax 

] 

(  M  ) 

n*1 
f 
j.41 

n*1 
-If       exp[ 
k     k 

2itk 

i  (J«1)   Ax 

X 

] 

(    A. 5    ) 

and  substituting  these  expressions   into 
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Next,   eq.A.2   is  divided  by  the   term  exp(  12irkjAx/Ji),   and  then  reduced  to 

M 
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(    A. 7    ) 

From  ea.A-3.    an  "amplification  factor"  Rk  for  each  K  can  be  obtained  as 

follows 

n+1 
2  Pjj  2AAt  2 

Rk  -  =  1  ♦  —  [21  Bin   (nkAx/>.)   ]  (  A. 8   ) 

F"*1  (Ax) 2 

then, 


8AAt  2  1/2 
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Since  3in^(  irkAx/  &)    is   always  positive  and  less   than  or  equal   to  one, 
the  requirement   that    |  R|<  |    <  1    is 
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Because  A  representing  momentum  or  ma3S  diffusion  coefficent  is 
positive,  AAt/(Ax)2  is  always  greater  than  zero.   Consequently, 


AAt    1 

<  -  (  A.!1  ) 

(Ax)2   H 


APPENDIX  B 
FLOW  REGIME  IN  A  CLOSED»END  CHANNEL 


B.1  Relationship  Between  Suspension  Concentration  and  Penalty 
In  the  present  study,  the  suspension  concentration,  C,  Is  defined 
as  the  sediment  mass  contained  in  a  unit  volume  of  water-*sediment 
mixture.  The  relationship  between  C  and  p  can  be  derived  as  follows: 


Ma 
p  =  ■  (  B.1  ) 

•a 


where  Ma  is  total  mass  of  the  water«sediment  mixture,  and  H0  Is  total 
volume  of  the  mixture.  However,  the  total  mass  Ma,  includes  the  mass 
of  fresh  water,  Mw,  and  the  mass  of  sediment,  Mg.  Then,  eq.B.1  becomes 


Mw  +  Ms     Pw(*o"»3)  +  Ps*s 
.  (  B.2  ) 


where  *3  is  the  volume  of  sediment  contained  in  the  mixture.  The  term 
on  the  right  hand  side,  p3Vs/V0,  i3  equivalent  to  the  suspension 
concentration,  C,  i.e.  p3«s/¥0  ■=  C.  Accordingly,  the  proportion  of 
volume  occupied  by  sediment  in  the  mixture  is  Vs/¥0  ■  C'p3  ■  C/pwG3. 
Consequently,  eq.B.J  can  be  expressed  as 


1 
C  (  I  n  —  )  i.  B.3  ) 


G 


3 
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or  the  'tensity  difference,    Ap,    hetween  the  water^sediment  mixture  and 
fresh  water  becomes 


1 

Ap   -    P   4    PW   -  C    (    1    *  —   )  (   B.4    ) 


B.2  Derivation  of  Steady  State  Horizontal  Velocity 
Equations  3.55  and  3.56  can  be  further  reduced  to 


32u     3h 

ez  -  g  —  n  <  z  <  h  (  B.5  ) 

3z2     3x 


32u     pw  3h        1   g  .n  3C 

-  g(— ) —  +  (1  -  — )— J  —  dz   0  <  z  <  n   (  B.6  ) 

3Z2     p   3x        Gs  p  z  3x 


The  integration  in  eq.B.6  can  be  carried  out  by  substituting  the 
solution  of  eq.3.52  into  the  suspension  concentration  C 


n  3C      n  3  wa1   x  w3l 

J  —  dz  =  J  —  {C,  expOB( )(-)]}dz  »  *g  (n  -  z)  C    (  B.7  ) 

z  3x       z  3x  u&   H  U4H 


Accordingly,  eq.B,6  becomes 

32u     pw  3h    Ap     ws1 

ez  -  g(— )—  +  g—  (-6  )(nftz)  (  B.8  ) 

3z2     p   3x    p      u^H 

where  Ap  «  p-pw.  Tlie  terms  on  the  right  hand  side  of  eq.B.5  and  eq.B.6 

are  not  functions  of  z,  hence  both  equations  can  be  integrated  twice 
with  respect  to  z  to  obtain 
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gli0   z2 

u  -    80   +    3i?    * n  <   z   <  h  (   B.9    I 

*z     2 

z2        65 
and  11  -    62   +    83  z  ♦    (6i)   *    1185) —  «  —  z3         0   <  z   <   n       (   B.10   ) 

?         6 

where  s0  Is  the  slope  of  free  water  surface  ■  3h/3x,  614  »  gpw3o''(ezf|) ■ 

65  -  SBgApWg^tpu^Hez),  and  all  of  8OI  »\,    62,  83  are  Integral  const* 

ants. 

Applying  the  boundary  condition  of  zero  shear  stress  at  water 

surface  (I.e.  3u/3x  -  0)  and  the  no*slip  boundary  condition  at  the  bed 

(I.e.  u  =  0)  to  eqs.B.9  and  B.10,  two  integral  constants  of  8i  and  82 

can  be  found  as  8-]  -  "gs0h/e2  and  82  -  0.  Consequently,  equations  B.9 

and  B.10  can  be  further  reduced  to 


gs0   z2 

u  -  (  —  «  hz  )  +  e0      n<z<h      (  B.11  ) 

ez    2 

z2   85 
and      u  =  83z  +  (8n  +  065)  —  »   —  z3   0  <  z  <  n      (  B.12  ) 

2    6 


Again,  two  other  boundary  conditions  at  the  Interface  are  utl« 
llzed.  These  are:  1)  the  horizontal  velocity  at  the  Interface  is 

continuous  (i.e.  ul  +  ■  ul   ),  and  2)  the  slope  of  vertical  velocity 

'  n     '  n- 

proflle  at  the  interface  is  continuous  (i.e.  3u/3x|  +  -  3u/3x|   ). 

1  n        '  n- 

Accordingly,  the  Integral  constants,  80  and  83,  In  eqs.B.11  and  B.12 
can  be  determined  as 


85       ga0       n 

—  n3+(  "81,)  —        (B.13) 

6         e,         2 
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Wo  K3o  95 

and        fo  -  1 h  -  ( — -  *  6i()  q  -  —  -t2     (  B.14  ) 

e,.      e.,  2 


Finally,  volumetric  conservation  requires  that  the  net  volumetric 
flow  at  any  cross  section  In  the  channel  must  be  equal  to  zero  (I.e. 


J  udz  -  'J).  This  Is  used  to  determine  6c.  Therefore  65  Is  found  as 
'0 

2gs0 

e5  ■  n (  B.15  ) 

h£3(i  M  0.25E.)ez 


where  5  is  the  dimen3ionjess  Interface  elevation  ■=  n/h. 

Substituting  e0,  63,  and  85  Into  eqs.  B.11  and  B.12,  the  solutions 
for  horizontal  velocity  in  upper  and  lower  layers,  respectively,  are 
obtained  as  follows: 

ks  H2  1 

u  =  [  0.5c2  -   C  + ]    5  <  C  <  1     (  B.16  ) 

ez  3(1  a   0.25C) 

gs0H2      jS  1  1 

u  = [ ♦  [0.5* ]r2  ♦  l"1+ J«  i 

nz       353(i"0.25£)        £2(1-0.25O  C(1--0.25E.) 

0  <  6  <  g   (  B.17  ) 
According  to  eq.B.16,  the  surface  velocity,  u3,  I.e.  r,  ■  1  ,  is 


8S0H  1 

u, i>0.5  * ]  (  B.I  8  1 

ez  3(1  *  0.25U 


us  is  used  to  normalize  eqs. B. 16  and  B.17.  Eqs. B.16  and  B.!7  tnen 
become 
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u     (1.5'0.375F.)  C2  1  (3KJ.75F.)  r  ♦  1 


t  <  C  <  1    (  ".19  ) 


»0.5  *  0.375  F, 
(C/t)3+C-3*1.5e2^0.37553)(c/5)2+(3«35+0.75C2)(5/C) 


0  <  c  <  F. 


u3  «0. 5  +  0.375  £ 

(  B.20  ) 

B.3  Comparison  of  Vertical  Velocity  Profiles 
A  comparison  of  vertical  velocity  profile  between  experimental, 
numerical,  and/or  analytical  results  is  made  in  this  section.   As 
mentioned  before,  the  vertical  velocity  profile  was  measured  by 
injecting  a  vertical  streak  of  dye  at  locations  approximately  1  cm 
away  from  the  side  wall  of  the  channel.  Numerical  predictions  of  the 
horizonal  velocity  were  obtained  from  model  simulation,  which  did  not 
Include  the  side  wall  resistance  effect.  Analytical  solutions  of 
steady  state  flow  regime  were  obtained  from  eqs.B.19  and  B.20. 
Figs. B. 1(a)  and  (b)  show  the  comparison  of  vertical  velocity  profiles 
between  measured  data  and  numerical  predictions  for  C0EL"6  test  unoer 
unsteady  conditions.  The  horizontal  velocity  is  normalized  by  |us|, 
where  u3  is  the  surface  outflow  velocity.  In  Fig. B. 1(a),  at  i  =  1.8  a 
and  t  •  40  minutes,  us  obtained  from  experimental  measurement  and 
numerical  simulation  are  equal  to  -0.36  cm/sec  and  *0.53  cm/see, 
respectively.   In  Fig.E.I(b),  at  x  -  1.0  m  and  C  =  40  minutes,  experi- 
mental and  numerical  us  are  equal  to  ^0.26  cm/sec  and  -0.25  cm/sec, 
respectively. 

Flgs.B.2(a;  and  (b)  show  the  comparison  of  vertical  velocity 
profiles  between  experimental,  numerical,  and  analytical  results 
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for  COEL-6  Lest  at  steady  state.   In  Fig. B. 2(a),  at  x  -  1.8  m,  us 
obtained  from  experimental  measurement  and  numerical  prediction  are 
equal  to  "0.18  cm/sec  and  1-0.53  cm/sec,  respectively.   us  of  analytical 
solution  can  not  be  computed  from  the  analytical  formulas  because  of 
the  involvement  of  the  water  surface  slope,  s0,  which  is  a  very  small 
value  and  not  easily  measurable.  In  Fig. B. 2(b),  at  x  ■  1.0  m,  the 
experimental  and  the  numerical  u3  are  equal  to  r»o.  1 1  cm/sec  and  -0.25 
em/sec,  respectively. 

It  is  observed  from  all  the  plots  mentioned  above,  that  the  shapes 
of  vertical  velocity  profiles  obtained  from  experimental,  numerical, 
and  analytical  results  are  fairly  similar.  However,  the  magnitude  of 
the  horizontal  velocity  obtained  from  laboratory  measurement  was 
smaller  than  that  from  numerical  prediction,  which  would  be  expected 
due  to  the  boundary  effect  of  the  side  wall. 
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APPENDIX  C 
COMPUTATION  Of  MAXIMUM  SURFACE  RISE 


When  a  stationary  sediment  wedge  Is  established,  the  maximum  rise 
of  water  surface  In  the  closed^end  channel  can  be  computed  using 
eq.3.62,  which  Is 

%        1  M        M  H2 

^max  -  -  {  —  J  C  j   Pi(z)  dz  *  —  }        (  3.M    ) 
H   py  0    z  2 

In  test  C0ELH5,  the  water  depth  at  entrance,  H-9.3  cm,  the  density 
of  water,  pw  «  0.996836  g/cm3,  and  the  density  profile  at  entrance, 
Pl(z),  was  measured.  Measured  values  of  p^  at  different  elevations 
at  the  entrance  of  the  closedftend  channel  are  given  in  Table  C.i, 
The  term  of  double  integrations  in  eq.3.62  can  be  expressed  as  a 
summation  according  to 


,H  M  8   8 

[J   p,(z)  dz]  dz  -  j.C  I   P1J  Azj  ]  AZl        (  C.1  ) 


f 

0   'z  1-1  J-I 


Computations  of  this  summation  are  given  in  Table  C.1. 
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Table  C.1  Computation  of  Maximum  Surface  Rise 


zl      Pll      Pti**J       IPIjAzj     (IpijAzi)4Zi 

J-i  j-i 

(era)        (kg/m-)  (kg/m2)  (kg/m2)  (kg/m) 


9.3         996.836 


8.5  997.116  7.977  7.977  0.032 


7.0         997.199  11.957  22.931  0.232 


5.5         997.213  11.958  37.892  0.156 


1.0         997.251  11.959  52.851  0.681 


2.5         997.251  11.959  67.810  0.905 


1.0         997.252  11.959  82.769  1.129 


0.0         997.252  9.970  92.710  0.878 


I    (    [pljAZj    )Azt    •   1.313 
i-1    j-i 


Therefore, 


1  1  (9.3)2 

^max  "  f x   1.313   x   10  N  }    =   0.021  mm 

9.3       0.996836  2 


APPENDIX  D 
CHARACTERISTICS  OF  CLOSED*END  CHANNEL  DEPOSITS 


Characteristics  of  the  deposits  of  closear-end  channel  at  COEL  are 
listed  in  the  following  table,  which  include  sample  location,  median 
dispersed  particle  size,  dm,  85?  finer  particle  size,  dgc,,   and  modified 
sorting  coefficient. 

Table  D.1  Characteristics  of  Closed^End  Channel  Deposits. 


Flume  Test 

Sample 

Location, 

dm 

d85 

Sorting 

Number 

Number 

X 

(i 

n) 

(um) 

Cum) 

Coefficient 

2 

1 

o.oou 

* 

0.50U 

4.12 

23.6 

2.1 

2 

0.00D 

H 

0.50D 

3.08 

15.5 

2.2 

3 

0.00 

M 

0.50 

1  .82 

9.6 

2.3 

1 

0.50 

* 

1.00 

1.10 

8.3 

2.1 

5 

1.00 

fi 

1.50 

1  .10 

6.8 

2.2 

6 

1.50 

ft 

2.00 

1.38 

5.6 

2.0 

7 

2.00 

H 

2.50 

1.18 

5.1 

2.1 

8 

2.50 

M 

3.50 

1.17 

5.3 

2.1 

9 

3.50 

M 

5.00 

1  .00* 

1.5 

2.1 

10 

5.00 

t» 

7.00 

— 

— 

— 

3 

1 

0.50U 

It 

0.50D 

1.32 

20.3 

2.2 

2 

0.00 

rt 

1.00 

2.50 

13.3 

2.3 

3 

1.00 

H 

2.00 

1.81 

8.0 

2.1 

4 

2.00 

p 

3.50 

1.H6 

6.1 

2.0 

5 

3.50 

- 

7.00 

0.96* 

1.6 

2.2 

1 

1 

0.75U 

H 

0.50D 

».m 

18.2 

2.1 

2 

0.00 

H 

3.00 

2.78 

10.7 

2.0 

3 

3.00 

" 

10.00 

0.91* 

1.3 

2.2 

*  indicates  extrapolated  data 

U  represents  a  location  at  the  upstream  of  the  main  channel  from 

the  center  of  the  clcsedflend  channel 
D  represents  a  location  at  the  downstream  of  the  main  channel  from 

the  center  of  the  closedrtend  channel 
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Table  D.1  —  continued 


Flume  Test 

Sample 

Location, 

dm 

d85 

Sorting 

Number 

Number 

X 

(m) 

(urn) 

(um) 

Coefficient 

5 

1 

o.oou 

« 

0.50U 

3.72 

23.3 

2.5 

2 

0.00 

4 

3.00 

2.78 

11.3 

2.2 

3 

3.00 

a 

10.00 

1.07* 

3.9 

1.9 

6 

1 

o.oou 

H 

0.50U 

22.90 

13.0* 

1.1 

2 

0.00 

I 

2.00 

12.10 

21.9 

1.1 

3 

2.00 

rt 

6.00 

6.01 

12.1 

1.1 

7 

1 

0.50U 

** 

0.50D 

21.55 

11.5* 

1.3 

2 

0.00 

* 

1.50 

13.11 

27.5 

1  .1 

3 

1.50 

n 

1.50 

6.21 

11.2 

1.3 

8 

1 

0.50U 

f 

0.50D 

29.17 

50.8* 

1.3 

2 

0.00 

f< 

2.00 

11.15 

31-5 

1.6 

9 

1 

O.OOU 

* 

0.50U 

2.89 

15.7 

2.3 

2 

0.00 

i 

3.00 

1.69 

9.1 

2.1 

3 

3.00 

t- 

10.00 

0.77* 

1.1 

2.1 

12 

1 

0.00 

it 

8.00 

9.71 

38.2 

2.0 

13 

1 

0.50U 

^ 

0.75U 

1.81 

11.2 

2.8 

2 

0.00 

» 

0.56 

1.39 

9.9 

2.7 

3 

0.63 

«■ 

2.63 

1.26 

8.7 

2.6 

1 

1.00 

^ 

7.00 

0.52* 

2.9 

2.1 

11 

1 

0.50U 

« 

0.75U 

1.83 

11.1 

2.8 

2 

0.00 

ft 

0.53 

1  .70 

13.2 

2.8 

3 

0.69 

H 

2.69 

1.71 

7.6 

2.1 

1 

1.00 

f 

7.00 

0.61* 

3.2 

2.3 

*  indicates  extrapolated  data 

U  represents  a  location  at  the  upstream  of  the  main  channel  from 

the  center  of  the  closedflend  channel 
D  represents  a  location  at  the  downstream  of  the  main  channel  from 

the  center  of  the  olosedfiend  channel 
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